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SECTION I 


SUMMARY 

A new simulation system, ROCETS, was designed and developed to allow cost-effective 
computer predictions of liquid rocket engine transient performance. The system allows a user to 
generate a simulation of any rocket engine configuration using component modules stored in a 
library thru high-level input commands. The system library currently contains 24 component 
modules, 57 sub-modules and maps, and 33 system routines and utilities. FORTRAN models from 
other sources can be operated in the system upon inclusion of interface information on comment 
cards. Operation of the simulation is simplified for the user by Run, Execution and Output 
Processors. The simulation system makes available steady-state trim balance, transient operation, 
and linear partial generation. The system utilizes a modern equation solver for efficient operation 
of the simulations. Transient integration methods include integral and differential forms for the 
trapezoidal, first order Gear, and second order Gear corrector equations. 

A detailed technology test bed engine (TTBE) model was generated to be used as the 
acceptance test of the simulation system. The general level of detail of the model was that 
reflected in the SSME DTM (Reference 2) . The model sucessfully obtained stady-state balance in 
main stage operation and simulated throttle transients including engine start and shutdown. A 
NASA fortran control model was obtained, ROCETS interface installed in comment cards, and 
operated with the TTBE model in closed-loop transient mode. 
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SECTION II 
INTRODUCTION 


The National Aeronautics and Space Administration (NASA) Facilities such as the George C. 
Marshall Space Flight Center (MSFC) require analysis and simulation of pump fed liquid rocket 
engine transient performance. The types of analysis and simulation include control design and 
analysis, design parametric studies, research and development, failure investigation, real-time 
simulation, feasibility studies, and software design, development, and testing. Therefore, multiple 
simulations representing different engine configurations with various levels of fidelity and transient 
response ranges are needed to support these studies. An analytical tool to meet these needs in a 
cost-effective manner is a digital computer simulation system. 


A computer simulation system named ROCket Engine Transient Simulation (ROCETS) was 
designed and developed under this program. An engine transient performance simulation normally 
consists of mathematical representations of the engine components interfaced together to describe 
the engine system performance. These component-by-component engine simulations (Figure 
2-1) require interfacing the component models together in a computer program, with appropriate 
program controls to interpret user commands, execute the program, and provide outputs to the 
user. All of this can be accomplished with in-line computer code that- is a free-standing 
simulation. However, a simulation system provides many benefits relative to individual, 
free-standing simulations. 
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Figure 2-1. A Rocket Performance Simulation Consist of 
Component-By-Component Models 

A simulation system (Figure 2-2) allows generation of simulations representing different engine 
configurations without expensive new computer code production and verification. The system acts 
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as a repository so that the same engineering methodology representing the components is utilized 
in different simulations to ensure prediction consistency. In addition, the system provides the 
latest modeling technology of verified numerical techniques and utilities; new advances placed in 
the system can easily be shared by all simulations operating in the system. A simulation system 
also provides a common operating base for all users to minimize required operational training after 
the initial start-up experience is obtained. 

■ Re-Use Of Developed/Verified Model Codes 

■ Repository For Methodology 

■ Advanced Modeling Technology & Techniques Easily Adaptable 

■ Reduces Required User Training 



Figure 2-2. Simulation Systems Are Effective Tools 

The ROCETS program to design and develop a simulation system consisted of nine (9) 
technical tasks: 

1 . Architecture 

2. System Requirements 

3. Component and Submodel Requirements 

4. Submodel Implementation 

5. Component Implementation 

6. Submodel Testing and Verification 
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7. Subsystem Testing and Verification 

8. TTBE Model Data Generation 

9. System Testing & Verification 

The Architecture definition determined there would be five major components of the ROCETS 
system: 

1 . Library System 

2. Executive programs (or Processors) 

3. Simulation Input and Output 

4. Documentation 

5. Maintenance Procedures 

The requirements were developed and documented in the System Requirements Specification 
(SRS) of P&W FR-20283, 25 November 1988 (Reference 3). The component and submodel 
implementation and testing/verification is contained in the System Design Specification (SDS) of 
P&W FR-20284, 25 July 1990 (Reference 4). The Technology Test Bed Engine (TTBE) Model 
description and system testing/verification are contained in this report. 
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SECTION III 

SYSTEM DESCRIPTION 


The Rocket Engine Transient Simulation (ROCETS) System was designed to use modular 
building blocks to represent engine components, and an architecture to interface these modules in 
any configuration desired by the user when generating an engine simulation. The architecture 
structure does not include any specific rocket engine configuration, and thus the flexibility exists to 
configure any rocket cycle of the future. The five components of the ROCETS system are: 

1. Library - A central source of all software code to allow multiple-users. 

2. Executive Programs - Software processors that conduct system functions. 

3. Simulation Input/Output - User inputs to configure a simulation, to execute the 
simulation, and to output the desired parameters. 

4. Documentation - System standards, engineering descriptions, user’s manual, 
programmers manual, and qualification test plants (contained in FR-20284). 

5. Maintenance Procedures - Instructions for system upkeep (contained in 
FR-20284) . 

3.1 SYSTEM OVERVIEW 

The ROCETS System has engineering models of all major engine components which are 
implemented as FORTRAN subroutines. These subroutines are called "modules". Standard 
engineering modules, once fully verified and documented, are put into a library so they can be 
accessed by all system users. A unique aspect of the ROCETS system is that engineering 
modules use comment cards to interface with the system. This allows ROCETS modules to be 
used outside the system as well as the ability to quickly adapt existing code to be used inside the 
system. 

Virtually any engine cycle can be represented by connecting the engineering modules in a 
desired order. While the modules could be connected by hand (i.e., an engineer building a main 
concatenating routine), this is time consuming, tedious, and error prone. The ROCETS system 
uses a Configuration Processor to accomplish this task (Figure 3-1). An engineer builds a 
configuration input file using high-level commands and the Configuration Processor generates an 
executable FORTRAN main program. The Configuration Processor also scans the execution order 
to identify algebraic loops required by the model. Algebraic loops are caused by variables which 
are used before being calculated, or variables which are outputs of more than one module. 

The ROCETS systems also uses a high-level command language to supply necessary inputs to 
run a particular model experiment. The input is read and interpreted by a Run Processor. The 
Run Processor initializes necessary inputs and sets appropriate flags to carry out the users 
instructions. The Run Processor allows the user to input schedules, set-up additional algebraic 
balances, and tailor a variety of integration options. 
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Figure 3-1 . ROCETS System Overview 

Execution control is provided by an Execution Processor. It controls looping, print, balancing, 
and linearization. Within the Execution Processor are calls to numerical utilities that provide 
steady-state balancing, transient integration, and linearization. It provides a centralized location for 
all numerical operations so that adding new features to the system in the future is simplified. 

Output is controlled by an Output Processor that allows the user to specify parameters to be 
printed and plotted. Plot information is supplied to a interface routine designed for a particular 
plotting software package. Therefore, implementing plotting on a different system requires only a 
change in the plot interface routine. Linearization output is not controlled by the Output Processor, 
but rather all necessary information is passed to a separate interface routine. This feature allows 
tailoring of the linearization output by changing only the interface routine. 

The ROCETS system has three run modes: Steady-state trim balance, Transient, and 
Linearization. The steady-state trim balance mode iterates dynamic states and algebraic balance 
variables until time derivatives and algebraic balance error terms are zero (within a specified 
tolerance). Transient integration normally integrates dynamic states using a predictor-corrector 
scheme with the corrector equations and algebraic balances closed simultaneously. The 
linearization mode linearizes about a steady-state or transient point and provides state-space 
model partials which can be used for other applications. Table 3-1 presents a summary of the 
ROCETS system significant features. 
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Table 3-1 . ROCETS System Significant Features 

Library 

1 . Storage for re-use of developed codes. 

2. Access for multi-users 

3. Repository for modeling methodology 

4. Allows adaptable future modeling technology 

Component Based 

1 . Component models are non-linear representation 

2. Generic component modules; unique characteristics in distinct sub-modules. 
Interface Structure 

1 . Component models use comment cards to Interface with system 

2. Existing models linked in any arrangement to simulate all engine cycles. 

3. Any FORTRAN model can easily be used in system simulations. 

Configurati on Processor 

1. Structured, English-like input 

2. Automatic scanning for required inputs and algebraic balances 

3. Generates FORTRAN main program 
Run Processor 

1. Structured, English-like input 

2. Schedule (curve) input available 

3. Definition of additional algebraic balances 

4. Activation for states and balances 

5. Three run modes: steady-state, transient, linearization 
Close-loop Integration With State-of-the Art Numerical Utilities 

1. Trapezoidal, first and second order Gear methods; others can be adapted. 

2. Ability to activate/deactivate states 

3. Ability to remove dynamic effects of states (force derivatives to zero) 

4. Advanced non-linear equation solver to close corrector and algebraic balances 
simultaneously. 

Automatic Linear Partial Generation 

1 . Repeatability and linearity checking 

2. Analytic handling of algebraic balances 

3. Analytic handling of states forced to steady-state. 

3.2 ENGINEERING MODULES 
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Engineering modules are stand-alone engineering representations of individual entities that are 
singular in purpose. The modular approach separates engineering modules, sub-modules, 
component data and generic data (properties) into the basic building blocks of the simulation. For 
example, a generic turbine module can be used multiple times in a single simulation simply by 
changing the component performance characteristics or map as well as being used in multiple 
simulations. This reduces the amount of code required while providing consistent methodology. 

The approach taken in modeling gives primary preference to engineering first principals 
followed by empirical correlations and transfer functions. However, modules of similar functions 
can be built with different modeling approaches and varying levels of complexity. The user then 
has the flexibility to select different approaches and level of detail used in a simulation. 

During the design phase of ROCETS, it was evident that the use of existing engineering 
representations would be desirable. To achieve this goal, it was decided to separate system 
functions from the engineering representations. This was accomplished by using call lists for 
communication to the engineering modules and keeping all system dependent code out of the 
individual modules. An additional benefit is that the modules can be operated as individual entities 
during design and verification. Modules only communicate to the ROCETS system through the 
subroutine call list. Commons are not used to communicate with the main or other modules. 
However, common blocks can be used in certain cases for communication between a module and 
a sub-module. 

Modules are interfaced to the ROCETS system using three blocks of comment cards at the 
beginning of the subroutine. These comment card blocks are called "interface cards" and are 
read by the Configuration Processor. The interface blocks relate call list names to system names, 
define the status of each variable for system operation, define the I/O status of each variable, and 
the FORTRAN variable type. Virtually any FORTRAN subroutine can easily be converted to the 
ROCETS system by adding the interface information on comment cards. However, the module 
history including author, dated revisions and internal code documentation should also be included. 

3.3 TRANSIENT MODELLING ASPECTS 

In general, the dynamics which are modelled in a rocket engine consist of volume dynamics, 
flow inertia, rotor speed integration, and thermal capacitance. Volume dynamics implement the 
laws of conservation of mass and energy using density and internal energy as dynamic states. 

Flow inertia dynamics implement conservation of momentum using flow rate as the dynamic state. 
Thermal nodes implement heat transfer laws and the energy equation applied to a metal mass 
using the metal temperature as the dynamic state. 

The baseline transient integration scheme is a predictor-corrector with the corrector equations 
closed by a modified Newton-Raphson iteration. Using a closed-loop integration offers advantages 
which are incorporated into the system. One item of particular usefulness is the capability of 
forcing states to their steady-state value during a transient. This is accomplished by using a 
steady-state error term (i.e., forcing the time derivative to zero) instead of closing the corrector 
equation for specified states. When this is done, the dynamic effects of the specified states are 
removed thereby allowing a variety of studies to be conducted. An obvious use of this feature is 
to obtain reduced order linear state-space models. However, it has also proven extremely 
valuable during model verification and validation. 

With the closed-loop integration, using density and internal energy as states causes numerical 
problems in liquid systems due to the extreme sensitivity to pressure to density and the difficulty in 
providing first guesses for internal energy. It would be considerably better to use pressure and 
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enthalpy as states but it is not possible to write appropriate differential equations. The solution to 
this problem is to make a change of iteration parameters. Instead of using density and internal 
energy as the iteration parameters to close the corrector equations, pressure and enthalpy are 
used. 

3.4 GLOBAL COMMUNICATION 

While engineering modules communicate through call lists, the system functions do not 
because ROCETS provides maximum flexibility by dynamically building system communication 
without using a predefined data structure. Therefore global commons are constructed which 
contain all the variables passed into or out of the engineering modules. These commons are used 
to communicate between the engineering modules, the interpretive reader, and the Execution 
Processor. Figure 3-2 depicts the communication flow. 

The global commons are divided by FORTRAN variable type: real, integer, character, double 
precision, complex, and logical. In addition to the 6 global variable commons, additional system 
commons are used to communicate information concerning states, derivatives, and additional 
balances as well as other necessary information. 



Figure 3-2. ROCETS Global Communication 
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3.5 CONFIGURATION PROCESSOR 


The goal when defining a simulation is to converts an abstract concept into a mathematical 
representation in a flexible, reliable, and convenient manner. Therefore it is desirable to automate 
the simulation creation to the extent possible, freeing the user from the tedious aspect of 
assembling a simulation. In ROCETS, a Configuration Processor is used to automatically create a 
simulation. (Figure 3-1). 

The configuration input consists of user commands defining a particular simulation in a simple, 
structured high-level format. The user defines the system to be modeled in the configuration input 
file by specifying component types, design characteristics, the relationship between various 
elements of the system, property packages to be used and what properties are to obtained, and 
definition of algebraic balances. (Note that algebraic balances can also be defined at run time). 

The processor performs two functions in generating a simulation: first it reads the 
configuration input, then it reads the interface definitions for the modules specified in the 
configuration. The processor needs the engineering module interfaces to determine the required 
variables, call lists, variable status (input, output, state, derivative, etc.) and variable type (real, 
integer, array, character, etc.). The processor cross references the configuration input and the 
module interface information to generate the specific variable names. These names are used to 
generate the appropriate call lists for the FORTRAN main program/module communication. This 
methodology is what allows the engineering modules to remain separate from the system code. 

The global commons are dynamically built for individual simulations during configuration. The 
commons consist of the variable names created from the module call lists along with required 
system variables. 

3.6 RUN PROCESSOR 

Run input consists of user commands to execute a configured simulation (Figure 3-1). The 
user input contains information required to define schedules, set inputs, define algebraic loops, 
specify output, and control execution. The input is in a high-level structured language. 

The ability to define and use schedules is quite powerful. Besides allowing schedules for time 
inputs, schedules can be set-up to define desired functional relationships in conjunction with 
algebraic loops. As an example, schedules can be defined to set a requested chamber pressure 
and mixture ratio and algebraic balances defined to vary valve areas until the requested values are 
obtained. 

Integration options can be tailored through run input to optimize model operation. The 
integration method, perturbation sizes, tolerance, convergence criteria, and activation can be set. 
The inputs are divided into "defaults” and "exceptions”. It is generally easier to set-up default 
information which is adequate for most states and then to override the defaults for specific states 
when necessary. Currently the system includes Euler, trapezoidal, first order Gear, and second 
order Gear integration schemes. However, other integration schemes can be easily added. 

The default for all states is to be active. However, it is often convenient to turn states off at 
various times. Three selections are possible for operations with each state: 

ON = the state is active 

OFF = the state is inactive and held constant 

STEADY-STATE = the state is always iterated to steady-state thereby removing the 

dynamic effect of the state 
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ROCETS provides the capability to define algebraic balances at run time. An independent 
variable can be varied until a dependent variable is equal to another dependent variable or until a 
dependent variable is equal to a value. The value may be an input or read from a schedule. This 
is especially useful when running operating lines or generating control schedules. 

Balance options can be tailored through run input similar to the integration options. The 
perturbation sizes, tolerance, convergence criteria, and activation can be set. Like the integration 
options, the inputs are divided into “defaults" and “exceptions". 

Linearization options can also be defined at run time. Included are values to be used for 
repeatability and linearity checking. Partials are generated by making a forward perturbation, 
backward perturbation, and repeating the forward perturbation. Repeatability is checked by 
comparing variable values on the two forward perturbation passes. If the percentage difference is 
more than the specified value a warning message will be written. Linearity is checked by 
comparing the forward and backward difference partials. If the percentage difference is more than 
the specified amount a warning message will be written. In addition to the check values, 
parameter names for linear model inputs and outputs are specified by the user. 

Linearization defaults are defined to establish the perturbation size for generating partials and 
exceptions to the defaults may also be specified. These functions are similar to default and 
exception declarations for states and balances. 

Simulation output options are handled through the Run Processor. This includes optional print 
during Jacobian evaluations and convergence attempts, options for debug output, and specifying 
simulation output. 

Three run modes are supported: Steady state, Transient, and Linearize. For steady-state, 
the number of consecutive steady-state points to be run is also specified. A system parameter 
POINT is available for reading schedules with the steady-state point number. This is useful for 
running steady-state operating lines or generating control schedules. For a transient, the time 
increment, print time, plot time, and termination time are specified. Time is available for reading 
schedules. The linearization mode perturbs each state and specified input to generate 
state-space model partial derivative matrices. 

3.7 MULTI-VARIABLE NEWTON-RAPHSON SOLVER 

ROCETS employs a state-of-the-art non-linear equation solver which is the heart of efficient 
system operation. It is a modified multi-variable Newton-Rahpson technique, which has been 
optimized to operate effectively with the large systems of equations encountered in the rocket 
modeling problems. The basic method operates on the matrix equation. 

AX = J-i AY 

Where AY is the amount the errors, or dependent variables, need to change to be zero and 
AX is the associated change in the independent variable. The solver Jacobian, J, is a matrix of 
partial derivatives generated with the model. The solver makes a number of passes equal to the 
number of iteration variables plus 1 through the model to generate the Jacobian. 

To improve the efficiency of the solver, the Jacobian is scaled using a modified version of the 
method given by McLaughlin (Reference 8) . The normalization factors from Jacobian scaling are 
then used in determining convergence as well as limiting allowed change of the independent 
variables. (This is necessary in non-linear systems to prevent excessive movement leading to 
exceeding map bounds, etc.) 
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Further enhancements include an algorithm following the Broyden method for updating the 
inverse Jacobian. Broyden' s method updates the inverse Jacobian without evaluating or inverting a 
new matrix, providing a large savings in number of total passes through the model. The matrix 
update is basically a secant-type method and is performed during convergence attempts. 

In steady-state operation the solver is used to drive all state time derivatives to zero while 
simultaneously driving algebraic balance error terms to zero (within a specified tolerance) . 
Transiently, the solver is used to provide simultaneous closed-loop integration of states and closing 
of algebraic loops. Closed-loop integration entails iteration on the simulation state variables (or 
state iteration variables) until they are equal to calculated values. This technique provides great 
flexibility since integration and algebraic balances are handled simultaneously. 

3.8 TRANSIENT INTEGRATION METHODS 

Rocket engine simulations comprise a set of stiff differential equations that require special 
methods for integration. Integral methods are efficient when the model time increment is small 
relative to the time constant associated with the state being integrated. However, as the model 
time increment is increased, a critical point is reached where convergence failure results. This 
limits the maximum time increment that can be used. 

An alternate method is to use the differential form of the corrector equations instead of the 
integral form, with the error term being formed as the difference between the actual and 
calculated derivative. Scaling by the time constant associated with each state is recommended by 
McLaughlin (Reference 8) . The differential method improves convergence when the model time 
increment is large compared to the state time constant. 

The integration routines used on the ROCETS system automatically uses the appropriate form of the 
corrector equation. The engineering modules approximate the time constant for each state and use 
this information to select the appropriate integration form. Both the integral and differential forms are 
incorporated for trapezoidal, first order Gear, and second order Gear corrector equations. 

3.9 LINEARIZATION 

The ROCETS system was designed to provide accurate linearization about a steady-state or 
transient operating point. Linearization provides state-space matrices of partial derivatives which 
can be used for subset model generation, transfer function creation, or multi-variable control 
analysis. 

Generation of accurate partial derivatives is critical to control design, analysis, and 
development. However, complications arise with large simulations using real properties and many 
dynamics components. These complications are due to changing iteration variables, the necessity 
to close algebraic balances, and from discontinuities associated with thermodynamics properties 
around the saturation dome. 

The ROCETS linearization methodology automatically accommodates any change of iteration 
variables for states and automatically closes all active algebraic balances. With the assumption of 
small perturbations such that the partials represent a linear model, the algebraic balances and 
state iteration parameters can be solved from a linear set of equations after all partials have been 
generated. 

The basic set of equations describing a linear model are: 

Xdot = A * X + B * U 
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Y = C * X + D * U 


where Xdot is a vector of state derivatives, Y is a vector of outputs, X is the state vector, and U is 
a vector of inputs. When using pressure and enthalpy (or internal energy) as iteration variables the 
X’s cannot be directly perturbed, so that the matrices cannot be measured directly. However, equations 
can be measured directly that allows for an analytic substitution of variables and solution of algebraic 
balances. 

Let the nomenclature be that T is a vector of state iteration variables and Z is a vector of 
algebraic balance independent variables. Then the equations that can be directly measured through 
perturbations are: 

Xdot = A1 * T + B$ta1 * U + Alpha2 * Z 
Errors = A3 * T + Beta2 * U + Alpha4 * Z 

Y = Cl * T + Zeta * U + Theta2 * Z 
X = Omega * T 

where the matrices represent the appropriate partial derivatives and Errors is a vector of error terms 
from algebraic balances which are to be zero. The change of variables from T’s to X’s is accomplished 
by solving the last equation for T and substituting Omega-inverse * X for T. Let 

Alphal = A1 * (Omega-inverse) 

Alpha3 = A3 * (Omega-inverse) 

Thetal = Cl * (Omega-inverse) 

so that the equations become; 

Xdot = Alphal * X + Betal * U + Alpha2 * Z 

Errors = Alpha3 * T + Beta2* U + Alpha4 * Z 

Y = Thetal * X + Zeta * U + Theta2 * Z 

Now solve for the algebraic balance parameters by noting that the error terms are to be zero. Then 
Z's are given by: 

Z » -(Alpha4-inverse) * (Alpha 3 * X + Beta 2 * U) 

Substitution yields: 

Xdot = (Alphal - Alpha2 * (Alpha4-inverse) * Alpha3) * X + 

(Betal - Alpha2 * (Alpha4-inverse) * Beta 2 ) * U 

y = (Thetal - Theta2 * (Alpha4-inverse) * Alpha3) * X + 

(Zeta - Theta2 * (Alpha4-inverse) * Beta 2 ) * U 

so that the actual matrices desired are given by; 

A = (Alphal - Alpha2 * (Alpha4-inverse) * Alpha 3) 

B = (Betal - Alpha2 * (Alpha-inverse) * Beta2 ) 

C = (Thetal - Theta2 * (Alpha4-inverse) * Alpha3) . 

D = (Zeta - Theta2 * (Alpha4-inverse) * Beta2 ) 
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In actual practice, an extremely accurate matrix inversion routine is necessary to preserve the 
integrity of the partials. A standard Gauss-Jordan reduction does not have sufficient accuracy. 
Therefore a Gauss-Jordan reduction has been combined with a recursion formula to obtain 
extremely accurate inversion and all matrix operations are performed in double precision. 

3.10 RUN TIME ERROR CHECKING 

Run-time error checking is provided to warn of possible invalid model conditions and run-time 
errors can also be used as a transient termination criteria. When curves are read with 
out-of-range inputs, when internal iterations fail, or any other condition that results in invalid 
conditions, the user is informed by appropriate messages in a "debug" file and a numerical status 
indicator (NSI) is set. 

Each error location is identified by two eight-character names called “module name” and 
“module location". The numerical status indicator is sent to a specific value depending on the 
error severity. Through run-time input, the user can control the NSI at which print will be provided 
and the NSI which is considered fatal. In addition to the NSI for fatal errors, the user supplies the 
number of occurrences of each fatal error before execution terminates. 


A list of error codes and their corresponding errors follows: 


0000 

1000 - 2999 
3000 - 4999 
5000 - 6999 
7000 - 9999 
10000 


3.11 DOCUMENTATION 


No error 

Map Extrapolation 

Input out of Range 

Internal Iteration Failure 

Invalid Solution 

Invalid Option (No Default), 

execution halted immediately by ERCK00 


ROCETS documentation starts in the module source code where the system standards require 
in the comments cards: A list of all module inputs & outputs with their definitions & units, an 
engineering description of the module, a list of sub-modules needed, and a history including 
qualification, author, and revision dates. The next level of documentation is contained in the 
ROCETS System Design Specification (SDS) of Ref. 4 which contains: 

Section 3.4 Documentation 


3.4.1 Standards 

3.4.2 Engineering Manual 

3.4.3 Programmer's Manual 

3.4.4 User’s Manual 

3.4.5 Qualification Test Plans 

3.12 ROCETS SYSTEM STATUS 

The ROCETS system software library contains approximately 100,000 lines of code of the 
executive programs (processors) and engineering modules and sub-modules, numerical utilities, 
and properites. The engineering generic modules represent engine components, and the 
sub-modules in general provide the specific component performance characteristics. Listed below 
are the 24 engineering modules representing the engine components with the corresponding 15 
sub-modules. A brief description of each follows. Engineering write-ups including all equations 
are contained in the SDS (Ref 2). A sample pump module “PUMP01” is presented in Appendix B. 
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ENGINE 

COMPONENTS 

MODULE 

SUB-MODULE 

PUMP 

PUMP01 

PMAP04 

TURBINE 

TURB0 1 

PMAP05 

PMAP06 

PMAP07 

PMAP08 

TBMP03 

TURBOPUMP 

TURB02 

ROTROO 

TBMP04 

TBMP05 

TBMP06 

PREBURNER 

ROTROI 

PBRN01 


MAIN CHAMBER 

MCHB01 


NOZZLE 

QCHM01 

NOZLOO 

CDNZOO 


NCLVOO 

QN0Z01 
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MODULE 


PLUMBING PIPEOO 

PIPE01 

PIPE02 FLPM02 

PIPE03 

PIPE04 

PIPE05 PRFP04 

PRFP06 

PRFP07 

PIPE06 

VOLMOO 

VOLMOI 

VOLM02 


VALVES VALVOO 

POGO SUPPRESSOR POGOOO 

GENERAL 

HEAT TRANSFER METLOO PRPM01 
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Module General Description Summary: 


MCHB01 H2/02 COMBUSTION AND VOLUME DYNAMICS WITH UNBURN CAPABILITY AND 
HELIUM DILUTION. 

METLOO ROUTINE IS A LUMPED MASS ANALYSIS OF A METAL WITH MULTIPLE HEAT 
TRANSFER NODES. 

NCLVOO ROUTINE FOR THE ENERGY ANALYSIS OF A LUMPED COOLING VOLUME USING 
DENSITY AND INTERNAL ENERGY AS STATES. 

NOZLOO CALCULATES FLOW AND THRUST FOR ISENTROPIC EXPANSION NOZZLE. 

PBRN01 PERFECT GAS COMBUSTION (H2/02) WITH VOLUME DYNAMICS AND HELIUM 
PURGE. 

PIPEOO CALCULATES THE FLOW DERIVATIVE AND CRITICAL TIME FOR INCOMPRESSIBLE 
FLUID FLOW IN PIPE WITH INERTIA AND LOSS. 

PIPE01 CALCULATE INCOMPRESSIBLE FLUID FLOW THROUGH A PIPE WITH A LOSS. 

PIPE02 CALCULATES COMPRESSIBLE FLUID FLOW THROUGH AN ORIFICE. 

PIPE03 CALCULATES THE FLOW DERIVATIVE AND CRITICAL TIME FOR INCOMPRESSIBLE 
FLUID FLOW IN PIPE WITH LOSS, INERTIA, AND CHANGE IN ELEVATION. 

PIPE04 CALCULATES UPSTREAM PRESSURE FOR LIQUID FLOW. 

PIPE05 CALCULATES UPSTREAM PRESSURE FOR COMPRESSIBLE FLOW. 

PIPE06 CALCULATES UPSTREAM PRESSURE FOR LIQUID FLOW. 

POGOOO MODELS THE PRIMARY DYNAMICS OF THE POGO SUPPRESSOR. 

PUMP01 ROUTINE REPRESENTS A CONSTANT DENSITY PUMP. 

QCHM01 CALCULATES HEAT TRANSFER RATE BETWEEN MULTIPLE METAL NODES AND THE 
HOT GAS FLOW PATH FOR ROCKET MAIN CHAMBER COOLING USING BARTZ 
CORRELATION. 

QN0Z01 CALCULATES HEAT TRANSFER RATE BETWEEN MULTIPLE METAL NODES AND THE 
HOT GAS FLOW PATH FOR ROCKET NOZZLE COOLING USING BARTZ 
CORPORATION. 

ROTROO CALCULATES THE ROTOR SPEED DERIVATIVE FOR A ROTOR SYSTEM. 

ROTROI ROTOR WITH BREAKAWAY TORQUE FOR STARTING SIMULATION. 

TURB0 1 ROUTINE IS AN ISENTROPTIC ANALYSIS OF A TURBINE THAT IS DRIVEN BY AN 

IDEAL GAS. 

TURB02 ROUTINE IS AN ISENTROPTIC ANALYSIS OF A TURBINE THAT IS DRIVEN BY A 
SINGLE CONSTITUENT FLUID. 

VALVOO CALCULATES INCOMPRESSIBLE FLUID FLOW THROUGH A VALVE USING LIQUID 
FLOW CORRELATIONS. 

VOLMOO ENERGY AND CONTINUITY ANALYSIS OF A VOLUME WITH ONE INLET MASS FLOW, 
ONE EXIT MASS FLOW AND ONE HEAT FLOW. 

VOLMOI GENERAL MULTI-FLOW LUMPED VOLUME FOR SINGLE CONSTITUANT FLUIDS USING 
DENSITY AND INTERNAL ENERGY AS STATES. 

VOLM02 ROUTINE IS USED FOR VOLUMES WITH MULTI-FLOWS, (BOTH IN AND OUT). IT • 
ASSUMES PERFECT GAS PROPERTIES AND CAN HANDLE FLOW REVERSALS. 
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CDNZOO ROUTINE CALCULATES VARIOUS PARAMETERS FOR A CONVERGENT-DIVERGENT 
NOZZLE. 

FLPM02 FLOW PARAMETER BASED ON TOTAL TO TOTAL PRESSURE RATIO AND NUMBER 
OF “VELOCITY HEADS” LOST. 

MACH03 CALCULATES MACH NUMBER FROM FLOW PARAMETER AND GAMMA FOR 
ADIABATIC FLOW OF A PERFECT GAS. 

MACH04 CALCULATES MACH NUMBER FROM AREA RATIO AND GAMMA USING AN 
ISENTROPIC RELATIONSHIP. 

PMAP04 ROUTINE DETERMINES THE PUMP CHARACTERISTICS FROM A MAP FOR THE 
ROCKETDYNE HIGH PRESSURE FUEL .PUMP. 

PMAP05 ROUTINE DETERMINES THE PUMP CHARACTERISTICS FROM A MAP FOR THE 
ROCKETDYNE LOW PRESSURE FUEL PUMP. 

PMAP06 ROUTINE DETERMINES THE PUMP CHARACTERISTICS FROM A MAP FOR THE 
ROCKETDYNE HIGH PRESSURE OXIDIZER PUMP. 

PAMP07 ROUTINE DETERMINES THE PUMP CHARACTERISTICS FORM A MAP FOR THE 
ROCKETDYNE LOW PRESSURE OXIIDIZER PUMP. 

PMAP08 ROUTINE DETERMINES THE PUMP CHARACTERISTICS FROM A MAP FOR THE 
ROCKETDYNE PREBURNER OXIDIZER PUMP. 

PRFP04 ROUTINE GIVES PRESSURE RATIO (TOTAL TO TOTAL) FROM FLOW PARAMETER, 
RKLS AND GAMMA USING TOTAL TEMPERATURE AND TOTAL PRESSURE. 

PRFP06 ROUTINE GIVES PRESSURE RATIO (TOTAL TO STATIC) FROM FLOW PARAMETER 
USING TOTAL TEMPERATURE AND STATIC PRESSURE. 

PRFP07 ROUTINE GIVES PRESSURE RATIO (TOTAL TO TOTAL) FROM FLOW PARAMETER, 
RKLS AND GAMMA USING TOTAL TEMPERATURE AND TOTAL DOWNSTREAM 
PRESSURE. 

TBMP03 ROUTINE DETERMINES THE TURBINE CHARACTERISTICS FROM MAPS FOR THE 
ROCKETDYNE HIGH PRESSURE FUEL TURBINE. 

TBMP04 ROUTINE DETERMINES THE TURBINE CHARACTERISTICS FROM MAPS FOR THE 
ROCKETDYNE LOW PRESSURE FUEL TURBINE. 

TBMP05 ROUTINE DETERMINES THE TURBINE CHARACTERISTICS FROM MAPS FOR THE 
ROCKETDYNE HIGH PRESSURE OXIDIZER TURBINE. 

TBMP06 ROUTINE DETERMINES THE TURBINE CHARACTERISTICS FROM MAPS FOR THE 
ROCKETDYNE LOW PRESSURE OXIDIZER TURBINE. 
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Listed below are the 29 ROCETS utilities subroutines which perform functions like integration, 
table reads, and error checks. A brief description of each follows the list. 


UTILITY FUNCTION 

ROUTINE 

DIRECT MODEL EXECUTION 

XPR000 

ERROR CHECKS 

ERCKOO 


ERCK01 


ERCK02 


ERCK03 

WRITES DUMMY GUESS ROUTINE 

DMGS01 

INPUT SCHEDULE PROCESSOR 

SPROOO 

INTEGREATION 

RINT01 

MATRIX OPERATIONS 

DPVN01 


DPVN03 


MTMU02 

OUTPUT 

LWRT01 


PRPL01 


WRIT01 

PARTIAL GENERATION 

LMRD01 


PART01 


PRTB01 

SOLVER 

ITER05 


SMIT03 


SSBL04 

SORT LIST 

SORTA4 

TABLE READ 

SUNBOO 


SUNB01 


SUNB03 


CPMR02 


CPMR04 


CPMR05 


CPMR06 

UNITS CONVERSION 

UNITOO 
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Numerical/System Utilities General Descirption Summary 


CPMR02 CORRESPONDING POINT BI-VARIENT MAP READER. 

CPMR04 CORRESPONDING POINT BI-VARIENT MAP READER THAT ALSO RETURNS PARTIALS. 

CPMR05 CORRESPONDING POINT BI-VARIENT MAP READER WITH NEW MAP INDEX POINTER. 

CPMR06 CORRESPONDING POINT BI-VARIENT MAP READER THAT ALSO RETURNS PARTIALS 

WITH NEW MAP INDEX DEFINITIONS. 

DMGS01 CREATES A DUMMY GUESS ROUTINE WITH LOCATIONS FOR THE APPROPRIATE 
REQUIRED GUESSES FOR A NEWLY CONFIGURED MODEL, 

DPVN01 DOUBLE PRECISION MATRIX INVERSION USING GAUSS-JORDAN ELIMINATION WITH 
PARTIAL MAXIMUM PIVOTING. 

DPNV03 DOUBLE PRECISION MATRIX INVERSION USING COMBINATION GAUSS-JORDAN AND 
RECURSION. FORMULA. 

ERCKOO ONE OF A PACKAGE OF FOUR ROUTINES TO PROVIDE RUN-TIME ERROR CHECKING 
(SEE ERCK01 , 02, AND 03). THIS ROUTINE IS CALLED AT THE POINT OF A ERROR 
TO PASS IN THE ERROR NUMBER ALONG WITH INFORMATION TO IDENTIFY THE 
ERROR. ERCK01 IS USED TO PROVIDE ERROR PRINT AND ERCK02 IS USED TO 
SPECIFY PRINT AND KILL LEVELS. 

ERCK01 ONE OF A PACKAGE OF FOUR ROUTINES TO PROVIDE RUN-TIME ERROR CHECKING 
(SEE ERCKOO, 02, AND 04). THIS ROUTINE IS CALLED TO PROVIDE OUTPUTING 
ERROR STATUS IN NORMAL TRANSIENT PRINT/PLOT AND FOR PRINTING 
END-OF-RUN ERROR STATUS. ERCKOO IS USED TO ENTER ERRORS AND ERCK02 
IS USED TO SPECIFY PRINT AND KILL LEVELS. 

ERCK02 ONE OF A PACKAGE OF FOUR ROUTINES TO PROVIDE RUN-TIME ERROR CHECKING 
(SEE ERCKOO, 01 , AND 03) . THIS ROUTINE IS CALLED BY THE USER TO SPECIFY 
THE PRINT LEVEL, KILL LEVEL, AND NUMBER OF FATAL ERRORS ALLOWED. 

ERCKOO IS USED AT THE POINT OF AN ERROR TO SET THE STATUS AND ERCK01 
IS USED TO PROVIDE ERROR PRINT. 

ERCK03 ONE OF A PACKAGE OF FOUR ROUTINES TO PROVIDE RUN-TIME ERROR CHECKING 
(SEE ERCKOO - ERCK02) . THIS ROUTINE IS CALLED TO SEARCH THE NAME 
ARRAYS AND RETURN THE LOCATION. IT IS A UTILITY FOR THE OTHER ROUTINES 
AND IS NOT USER CALLABLE. 

ITER05 SECANT METHOD ITERATION WHICH CAN BE USED WITH NESTED LOOPS 

LMRD01 ROUTINE WRITES OUT THE RESULTS OF THE LINEARIZATION OF ROCETS 

SIMULATION. 

MTMU02 DOUBLE PRECISION MATRIX MULTIPLICATION. 

OPCK01 TRANSFORMS THE PROPERTY OPTION CHARACTER FIELD INTO A STANDARD FORM 
AND TO READ THE INDEPENDENT & DEPENDENT PROPERTY NAMES. 

PART01 MEASURES AND CHECKS THE PARTIALS FOR THE ROCKET ENGINE TRANSIENT 
SIMULATION SYSTEM. 
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PRPL01 A GENERALIZED PRINT/PLOT ROUTINE FOR TRANSIENT DECKS PROVIDING 

COLUMNAR PRINT, INTERFACING FOR PLOTS, AND OPTIONAL' USER PRINT/PLOT 
HEADER SPECIFICATION AS WELL AS TAILORED PRINT FORMAT. 

PRTB01 ROUTINE PERTURBATES THE INDEPENDENT VARIABLE (X) AND MEASURE THE 
PARTIAL OF THE VECTOR OF DEPENDENT VARIABLE (Y) WITH RESPECT TO 
INDEPENDENT VARIABLE. 

RINT01 ROUTINE PERFORMS THE CLOSED LOOP INTEGRATION FOR THE ROCKET ENGINE 
TRANSIENT SIMULATION SYSTEM. A CHOICE OF THREE IMPLICIT INTEGRATION 
TECHNIQUES (TRAPEZOIDAL, FIRST ORDER GEAR, AND SECOND ORDER GEAR) IS 
INCLUDED. AN EULER INTEGRATION IS ALSO AVAILABLE. 

SMIT03 ROUTINE SOLVES A SET OF SIMULTANEOUS NONLINEAR EQUATIONS USING 
NEWTON’S METHOD WITH BROYDEN’S INVERSE JACOBIAN UPDATE SCHEME. 

SORTA4 ROUTINE SORTS A LIST OF CHARACTER WORDS. 

SPR000 SCHEDULE PROCESSOR WHICH PROCESSES RUN-TIME SCHEDULES FOR ROCTES 
SIMULATION SYSTEM. 

SSBL04 THE COMPANION ROUTINE TO THE TRANSIENT INTEGRATION ROUTINE (RINT01). 
SSBL04 IS USED TO ACHIEVE A STEADY STATE BALANCE OF STATES AND 
CONVERGENCE OF ADDITIONAL BALANCES. 

SUNBOO UNI-VARIENT OR BI-VARIENT SEPARATE INTERPOLATION MAP READER WITH 
OPTIONAL EXTRAPOLATION FOR OUT-OF-RANGE DATA. 

SUNB01 UNIVARIENT OR BIVARIENT SUNBAR-TYPE MAP READER WITH OPTION TO READ 

MAP IN ANY DIRECTION AND EITHER EXTRAPOLATE OR RETURN CORNER VALUES. 

SUNB03 TRI-VARIENT SEPARATE INTERPOLATION MAP READER WITH OPTIONAL 

EXTRAPOLATION FOR OUT-OF-RANGE DATA AND MULTI-DIRECTION READ OPTION. 

UNITOO ROUTINE SETS CONVERSION FACTORS AND CONSTANTS IN THE UNITS COMMON 
BASED ON THE SYSTEM REQUESTED. 

WRIT01 INTERFACE ROUTINE FOR PRP01 TO WRITE PLOT FILE FOR MSFC IBM 3080 
SYSTEM IN THE ORIGINAL UNIVAC FILE FORMAT. 

XPR001 PURPOSE: INTERFACES WITH SYSTEM COMMONS TO DIRECT MODEL EXECUTION. 
THIS ROUTINE CALLS NECESSARY NUMERICAL ROUTINES AND DIRECTS 
EXECUTION IN THE CALLING PROGRAM BY MEANS OF AN OUTPUT SIGNAL. 

The properties contained in the ROCETS system include combustion for H 2 /O 2 , hydrogen, 
oxygen, helium, nitrogen, methane, and various metals. They are listed below and followed with a 
brief description. 
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PROPERTY 

MODULE 

SU9-IVIQPULE 

COMBUSTION 

HGPROP 

COMB0 1 



COMB02 

ZGASOO 

ZZH201 

HYDROGEN 

H2PROP 

HHCP05 



HHCV05 
HHPK05 
. HHPS05 
HHPU05 
HPSH01 
HPUT05 
HRHP05 



HRHP06 

HRUP05 

OXYGEN 

02PROP 

OHCP05 

OHCV05 

OHPK05 

OHPS05 

HHPU05 

HPSH01 

HPUT05 

HRHP05 

HRHP06 

HRUP05 

HELIUM 

HEPROP 

EHPS05 

EHPT05 

ERHP05 

NITROGEN 

N2PROP 

NHPS05 

NPHT05 

NRHP05 

METHANE 

N2PROP 

MCPT05 

MCVT05 

MHPS05 

MHPT05 

MRHP05 

METALS 


PRPM01 
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Properties General Description Summary: 

COMBOI OBTAINS HOT GAS TRANSPORT PROPERTIES FOR H2/02 COMBUSTION PRODUCTS'. 

COMB02 PERFECT GAS COMBUSTION (H2/02) WITH HELIUM DILUTION. 

EHPS05 OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR HELIUM FROM AN ENTHALPY, 
PRESSURE, ENTROPY MAP. 

EPTH05 OBTAINS THEMOPHYSICAL FLUID PROPERTIES FOR HELIUM FROM A PRESSURE. 
ENTHALPY TEMPERATURE MAP. 

ERHP05 OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR HELIUM FROM A DENSITY, 
ENTHALPY, PRESSURE MAP. 

HEPROP SUBROUTINE ACCESSES HELIUM PROPERTIES VIA HELIUM PROPERTY MAPS. 

HGPROP MAIN DRIVER FOR HOT GAS PROPERTIES (H2/02 COMBUSTION PRODUCTS) FROM 
MAPS (INCLUDES HELIUM PURGE EXCEPT FOR Z, MU, AND K) 

HHCP05 OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN FROM AN 
ENTHALPY, PRESSURE, CONSTANT PRESSURE SPECIFIC HEAT MAP 

HHCV05 OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN FROM AN 
ENTHALPY, PRESSURE, CONSTANT VOLUME SPECIFIC HEAT MAP 

HHPK05 OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN FROM AN 
ENTHALPY, PRESSURE, THERMAL CONDUCTIVITY MAP 

HHPS05 OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN FROM AN 
ENTHALPY, PRESSURE, ENTROPY MAP. 

HHPU05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN 
FORM AN ENTHALPY, PRESSURE, VISCOSITY MAP. 

HPSH01 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN 
FROM A PRESSURE, ENTROPY, ENTHALPY MAP. 

HPUT05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN 
FROM A DENSITY, ENTHALPY, PRESSURE MAP. 

HRHP05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN 
FROM A DENSITY, ENTHALPY, PRESSURE MAP. 

HRHP06 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN 
FROM A DENSITY, ENTHALPY, PRESSURE MAP. 

HRUP05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR PARA-HYDROGEN 
FROM A DENSITY, INTERNAL ENERGY, PRESSURE MAP. 

H2PROP SUBROUTINE ACCESSES HYDROGEN PROPERTIES VIA HYDROGEN PROPERTY 
MAPS. 

MCPT05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR METHANE FROM 
A TEMPERATURE, PRESSURE, CONSTANT PRESSURE SPECIFIC HEAT MAP. 

MCVT05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR METHANE FROM 
A TEMPERATURE, PRESSURE, CONSTANT VOLUME SPECIFIC HEAT MAP. 

MEPROPR SUBROUTINE ACCESSES METHANE PROPERTIES VIA METHANE PROPERTY MAPS. 

MHPS05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR METHANE FROM 
AN ENTHALPY, PRESSURE, ENTROPY MAP. 


25 



MPHT05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR METHANE FROM 
AN ENTHALPY, PRESSURE, TEMPERATURE MAP. 

MRHP05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR METHANE FROM 
A DENSITY, ENTHALPY, PRESSURE MAP. 

NHPS05 SUBROUTINE OBTAINS OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR 
NITROGEN FROM AN ENTHALPY, PRESSURE. ENTHROPY MAP. 

NPHT05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR NITROGEN FROM 
A PRESSURE, ENTHALPY TEMPERATURE MAP 

NRHP05 SUBROUTINE SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR 
NITROGEN FROM A DENSITY, ENTHALPY, PRESSURE MAP 

N2PROP SUBROUTINE ACCESSES NITROGEN PROPERTIES VIA NITROGEN PROPERTY MAPS. 

OHCP05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR OXYGEN FORM 
AN ENTHALPY, PRESSURE, CONSTANT PRESSURE SPECIFIC HEAT MAP 

OHCV05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR OXYGEN 
FROM AN ENTHALPY, PRESSURE, CONSTANT VOLUME SPECIFIC HEAT MAP 

OHPK05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR OXYGEN FROM 
AN ENTHALPY, PRESSURE, ENTROPY MAP. 

OHPS05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR OXYGEN FROM 
AN ENTHALPY, PRESSURE, ENTROPY MAP. 

OHPU05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR OXYGEN FROM A 
PRESSURE, ENTROPY, ENTHALPY MAP. 

OPUT05 SUBROUTINE THERMOPHYSICAL FLUID PROPERTIES FOR OXYGEN FROM A 
PRESSURE, INTERNAL ENERGY, TEMPERATURE MAP 

ORHP05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR OXYGEN FROM A 
DENSITY, ENTHALPY, PRESSURE MAP 

ORHP06 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES AND PARTIALS FOR 
OXYGEN FROM A DENSITY, ENTHALPY. PRESSURE MAP 

ORUP05 SUBROUTINE OBTAINS THERMOPHYSICAL FLUID PROPERTIES FOR OXYGEN FROM A 
DENSITY, INTERNAL ENERGY, PRESSURE MAP 

02PROP SUBROUTINE ACCESSES OXYGEN PROPERTIES VIA OXYGEN PROPERTY MAPS 

PRPM01 SUBROUTINE ACCESSES OXYGEN PROPERTIES VIA OXYGEN PROPERTY MAPS. 

PRPM01 SUBROUTINE GIVES THE SPECIFIC HEAT AND CONDUCTIVITY OF VARIOUS METALS 
AS A FUNCTION OF TEMPERATURE. 

ZGASOO SUBROUTINE CALCULATES THE REAL GAS COMPRESSIBILITY FACTOR FOR H2/H202 

ZZH201 SUBROUTINE CALCULATES COMPRESSIBILITY FACTOR FOR H2. 
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SECTION IV 
TTBE MODEL 

Two models of the Technology Test Bed Engine (TTBE) were generated under the program. The 
initial model was a simple model without boost turbopumps, and with a simulation complexity of 55 
state variables and 2 algebraic loops. After testing and verification of this model, a detailed TTBE model 
with the boost turbopumps and a POGO system was configured with 122 state variables and 14 algebraic 
loops. 

4.1 SIMPLE TTBE SIMULATION 

A simple model of the Technology Test Bed Engine (TTBE) was generated as the initial system 
verification vehicle for the simulation system. Figure 4-1 shows a schematic of the simple TTBE 
along with the 42 specific stations in the simulation. By using generic code, only the following 13 
component modules were required by the simulation: 

1 . INJTOO - Main Injector 

2. MCHBOO - Main Chamber 

3. MIXROO - Flow Mixer 

4. NOZLOO - Nozzle Thrust Calculations 

5. PBRNOO - Preburner 

6. PIPEOO - Incompressible flow with inertia 

7. PIPE01 - Incompressible flow without inertia 

8. PIPE02 - Compressible flow without inertia 

9. PUMOO - Polytropic Pump 

10. ROTROO - Rotor Torque Balance/Speed Derivative 

11. SPLT00 - Flow Splitter 

12. TURB00 - Turbine 

13. VOLMOO - Volume 

The modules described above were configured into the simple TTBE simulation along with 
required property relationships and numerical utilities. There were 55 state variables and 2 
algebraic loops required in the simulation as shown in Table 4-1. State derivatives and outputs are 
calculated form model inputs and states. 

Using initial guesses from data of the Digital Transient Model (DTM) of Reference 2 at 100% 
RPL, SMITE successfully obtained all TTBE model state derivatives and algebraic loop parameters 
to within specific tolerances. This demonstrated the capability of the ROCETS system to converge 
a rocket simulation to a steady-state point without running a transient. Transient capability was 
demonstrated by running the simple TTBE simulation with small perturbations of valve areas about 
the 100% RPL point. The results of these tests are presented in Section 5.0 - System testing and 
verification. 
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ROCETS SIMPLIFIED TTBE MODEL 
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Figure 4-1. ROCETS Simplified TTBE Model 





TABLE 4-1 - Simple TTBE Simulation States and Algebraic Loops 


STATE 

DESCRIPTION 

OFRFPRB 

Fuel Preburner Oxidizer Fraction 

OFRMCHB 

Main Chamber Oxidizer Fraction 

OFRMFI 

Main Fuel Injector Oxidizer Fraction 

OFROPRB 

Oxidizer Preburner Oxidizer Fraction 

RHOFMCO 

Volume FMCO Density 

RHOFPBI 

Fuel Preburner Injector Density 

RHOFPRB 

Oxidizer Preburner Injector Density 

RHOFIOS 

Volume FI OS Density 

RHOMCHB 

Main Chamber Density 

RHOMFI 

Main Fuel Injector Density 

RHOMOI 

Main Oxidizer Injector Density 

RHOOPBI 

Oxidizer Preburner Injector Density 

RHOOPRB 

Oxidizer Preburner Density 

RHOPBSF 

Preburner Fuel Splitter Density 

RHOPBSO 

Preburner Oxidizer Splitter Density 

RH010 

Volume 10 Density 

RH012 

Volume 12 Density 

RH015B 

Volume 15B Density 

RH08 

Volume 8 Density 

SNF2 

Fuel Turbomachinery Speed 

SN02 

Oxidizer Turbomachinery Speed 

TTFPRB 

Fuel Preburner Temperature 

TTMCHB 

Main Chamber Temperature 

TTMFI 

Main Fuel Injector Temperature 

TTOPRB 

Oxidizer Preburner Temperature 

UTFMCO 

Volume FI OS Internal Energy 

UTMOI 

Main Oxidizer Injector Internal Energy 

UTOPBI 

Oxidizer Preburner Injector Internal Energy 

UTPBSF 

Preburner Fuel Splitter Internal Energy 

UTPBSO 

Preburner Oxidizer Splitter Internal Energy 

UT10 

Volume 10 Internal Energy 

UT12 

Volume 12 Internal Energy 

UT15B 

Volume 15B Internal Energy 

UT7S 

Volume 7S Internal Energy 

UT8 

Volume 8 Internal Energy 

WCCV 

Coolant Control Valve Flow 

WFPHI 

Fuel Preburner Fuel Flow 

WFPOV 

Fuel Preburner Oxidizer Flow 

WFTI 

Line FTI Flow 

WFTID 

Line FTID Flow 

WF10 

Line F10 Flow 

WMFV 

Main Fuel Valve Flow 

WOPHI 

Oxidizer Preburner Fuel Flow 

WOPOV 

Oxidizer Preburner Oxidizer Flow 

W11 

Line 1 1 Flow 

W16 

Line 16 Flow 
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W20 

W21 

W4 

W7 

W9 


Line 20 Flow - 
Line 21 Flow 
Line 4 Flow 
Line 7 Flow Flow 
Line 9 Flow 


ITERATION VARIABLE DESCRIPTION 


WFTRB Fuel Turbine Flow - Iterated until equal to 

calculated value 

WOTRB LOX Turbine flow - Iterated until equal to 

calculated value 

4.2 DETAILED TTBE SIMULATION 

After successful verification of the simple TTBE model, a detailed TTBE model simulation was 
developed. The approach was to model the lox side and test, then the fuel side and test, then 
the hot gas system and test, and finally connect the three sub-systems and test. A schematic of 
the entire engine simulation is presented on Figure 4-2. There are 122 states and 14 additional 
balances in the simulation. Each of the station names are labeled on the schematic. 

A description of the modules used to configure the TTBE along with a list of the TTBE 
schematic names that use that particular module follows. 
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Figure 4-2. ROCETS TTBE Detailed Configuration 












Pipe Modules 


Six pipe modules were used to configure the TTBE. They are PIPEOO, PIPE01 , PIPE02, PIPE03, 
PIPE05, PIPE06. Following is a list of the TTBE schematic pipe names that use the corresponding 
pipe routines, a description of each module, and the states of each module. 

Module PIPEOO 

Schematic Names : FI, F2, F5, F6, F7, F8, FFPB, FOPB, 02, 03, 05 

OFPB, OOPB, OINJ 

Description: 

This module calculates the flow derivative for incompressible fluid flow in a pipe with 
inertial effects. 

States and State Derivatives: 

PIPEOO has flow as a state and calculates the flow derivative. 

Module PIPE01 

Schematic Names: F3, F4, F11, FI 2, F13, F14, F15, FSLV, FIG. FTC. 

OTC, 04, 06, 07, MOV, FPOV, OPOV 


Description: 

This module calculates flow through a pipe for an incompressible fluid. 

States and State Derivatives: 

PIPE01 calculates flow but it is not treated as a state. This module does not treat any 
parameters as states and therefore does not calculate any state derivatives. 

Module PIPE02 

Schematic Names: HG3, HG4, OLK 

Description: 

This module calculates flow through a pipe for a compressible fluid. 

States and State Derivatives: 

PIPE02 calculates flow but it is not treated as a state. This module does not treat any a 
parameters as states and therefore does not calculate any state derivatives. 

Module PIPE03 

Schematic Names: 01 

Description: 

This module calculates the flow derivative for incompressible fluid flow in a pipe with 
inertial effects and elevation change. 
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States and State Derivatives; 


PIPE03 has flow as a state and calculates the flow derivative. 

Module PIPE05 

Schematic Names: HG1, HG2, HG3, HG6, FINJ, HG5 

Description: 

This module calculates upstream pressure for compressible flow. 

States and Sta te Derivatives: 

PIPE05 uses flow but it is not treated as a state. This module does not treat any 
parameters as states and therefore does not calculate any state derivatives. 

Module PIPEQ6 

Schematic Names: F9, F10 

Description: 

This module calculates upstream pressure for incompressible flow. 

States and State Derivatives: 

PIPE06 uses flow but it is not treated as a state. This module does not treat any 
parameters as states and therefore does not calculate any state derivatives. 

Yalve Mo dule 

One valve module, VALVOO, was used to configure the TTBE. Following is a list of the TTBE 
schematic valve names that use the VALVOO routine, a description of the VALVOO module, and the 
states of the VALVOO module. 

Module VALVOO 

Schematic Names : MFV, CCV 

Description: 

This module calculates flow through VALVE for an incompressible fluid. 

* 

States and State Derivatives : 

VALVOO calculates flow but it is not treated as a state. This module does not treat any a 
parameters as states and therefore does not calculate any state derivatives. 
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Volume Module 


Module VOLMQO 

Schematic Names: VL1 , VL2, VL5, VL9, VL10. VL14, VL15, VL17, VL18, VL22 

Description: 

This module performs a continuity and energy analysis of a volume for pure fluids with 
one inlet flow, one exit flow, and one heat flow. 

States and State D erivatives: 

VOLMOO has density and internal energy as states. Corresponding derivatives are 
calculated for each state. 

Module VOLMQ1 

Schematic Names: VL3, VL4, VL8, VL13, VL16, VL19, VL20, VL21, PBSF, PBSO 

Description: 

The module performs a continuity and energy analysis of a volume for pure fluids with 
multiple inlet flows, multiple exit flows, and multiple heat flows. 

States and State Derivatives: 

VOLMOI has density and internal energy as states. Corresponding derivatives are 
calculated for each state. 

Mo dul e YQ LMQ2 

Schematic Names: FTBP, OTBP, FSF, OSF, MFI 

Description: 

This module performs a continuity and energy analysis of a volume for perfect gases with 
oxygen, hydrogen, and helium as possible constituents. The analysis is performed with 
multiple inlet flows, multiple exit flows, and multiple heat flows. 

States and State Derivatives: 

VOLM02 has pressure, temperature, oxidizer fraction, and helium fraction as states. 
Corresponding derivatives are calculated for each state. 
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Module NCLVQQ 


Schematic NAMES; VL6, VL7, VL11, VL12 
Description: 

This module models the cooling of the chamber and nozzle. The module performs a 
continuity and heat transfer analysis of a volume of pure fluids with one inlet flow, one 
exit flow, multiple node metal temperatures, and multiple node heat transfer surface 
areas. The heat transfer rate is calculated for each node. 

States and State Derivatives: 

NCLVOO has density and internal energy as states. Corresponding derivatives are 
calculated for each state. 

Rotor Module 

The inertial and transient speed effects for turbopumps are modeled by the ROTROO 
module which mechanically links each turbopump together. Following is a list of the TTBE 
schematic turbopump names that use the rotor routine, a description of the ROTROO 
module, and the states of the ROTROO module. 

Module ROTROO 

Schematic Names: LPFP/LPFT, HPFP/HPFT, LPOP/LPOT, (HPOP and PRBP)/HPOT 
Description: 

Given supply torques, required torques, rotative speed and the overall polar moment of 
inertia, this routine calculates the speed derivative for the given system. 

States and state Derivatives: 

ROTROO has rotative speed as a state and calculates the corresponding speed derivative. 
Pump Module 

One pump module, PUMP01 , was used to configure the TTBE. Following is a list of the TTBE 
schematic pump names that use the pump routine, a description of the PUMP01 module, and the 
states of the PUMP01 module. 

« 

Modul e PU MP Q1 

Schematic Names: LPFP, HPFP, LPOP, HPOP, PRBP 
Description: 

By assessing the appropriate pump performance map, this module calculates exit 
enthalpy, exit pressure, and required torque for a constant density pump. 


35 



States and State Derivatives: 


PUMP01 uses rotative speed but it is not treated as a state (see module ROTROO above) . 
This module does not treat any parameters as states and therefore does not calculate 
any state derivatives. 

Turbine Modules 

Two turbine modules were used to configure the TTBE and they are TURB01 and TURB02. 
Following is a list of the TTBE schematic turbine names that use the corresponding turbine 
routines, a description of each module, and the states of each module. 

Module TURB0 1 

Schematic Name: LPFT, HPFT, HPOT 
Description; 

By accessing the appropriate turbine performance map, using isentropic efficiency this 
module calculates exit enthalpy, supply torque, and required turbine flowrate for a turbine 
driven by a perfect gas. 

States and State Derivatives: 

TURB01 uses rotative speed but it is not treated as a state (see module ROTROO above) . 
This module does not treat any parameters as states and therefore does not calculate 
any state derivatives. 

Module TURB02 

Schematic Name: LPOT 

. Description: 

By accessing the appropriate turbine performance map, using intropic efficiency this 
module calculates exit enthalpy, supply torque, and required turbine flowrate for turbine 
driven by a liquid. 

States and State Derivatives: 

TURB02 uses rotative speed but it is not treated as a state (see module ROTROO above). 
This module does not treat any parameters as states and therefore does not calculate 
any state derivatives. 

EQQQ Module 

A POGO module, POGOOO, was used to configure the POGO suppression system for the TTBE. 
Following is a list of the TTBE schematic component names that use the POGO routine, a 
description of the POGO module, and the states of the POGO module. 
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Module POGOQQ 


Schematic Name: POGO 
Description; 

This module models the POGO suppression system. Given the oxygen-side conditions, 
this module calculates the required exit oxygen flowrate and appropriate derivatives. 

States and State Derivatives; 

POGOOO has pressure, liquid oxygen flowrate, liquid oxygen mass, and helium fraction as 
states. Corresponding derivatives are calculated for each state. 

Main Chamber Combustion Module 

One main chamber combustion module, MCHB01 , was used to configure the TTBE. Following 
is the TTBE schematic name for the main chamber combustion, a description of the MCHB01 
module, and the states of the MCHB01 module. 

Module MCHB01 

Schematic Names: MCHB 

Description: 

This module models perfect gas hydrogen/oxygen combustion with helium dilution, unburn 
capability, and volume dynamics. 

States and State Derivatives: 

MCHB01 has pressure, temperature, oxidizer fraction, and helium fraction as states. 
Corresponding derivatives are calculated for each state. 

Preburner Module 

One preburner module, PBRN01 , was used to configure the TTBE. Following is a list of the 
TTBE schematic preburner names that use the PBRN01 routine, a description of the PBRN01 
module, and the states of the PBRN01 module. 

Module PBRNQ1 

Schematic Names: FPRB, OPRB 
Description: 

This module models perfect gas hydrogen/oxygen combustion with helium dilution, and 
volume dynamics. 

States and State Derivatives; 

PBRN01 has pressure, temperature, oxidizer fraction, and helium fraction as states. 
Corresponding derivatives are calculated for each state. 
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METAL MODULE 


The transient metal temperature effectors are modeled by the METLOO module. Following is a 
list of the TTBE schematic component names that use the METLOO routine, a description of the 
METLOO module, and the states of the METLOO module. 


Module METLOO 

Schematic Names: VL6/QDOTNZ1 , VL7/QDONTNZ2. VL1 1/QDOTCHM1 , VL1 2/QDOTCHM2. 

BfiSgtifltiflni 

Given the mass of the metal and the temperature of the metal, this routine calculates the 
metal temperature derivative for each of the given multiple nodes. 

States and State Derivatives: 

METLOO has metal temperature as a state and calculates the corresponding derivative. 
Chamber Hot Side Heat Transfer Module 

One chamber heat transfer module, QCHM01 , was used to configure the TTBE. Following is a 
list of the TTBE schematic names that used the QCHM01 routine, a description of the QCHM01 
module, and the states of QCHM01 module. 

Module Q CH MQ1 

Schematic Names: QDOTCHM1 , QDOTCHM2 
Description: 

This module calculates the multiple node heat flowrate through the chamber wall from the 
combustion gases using a Bartz empirical correlation. 

States and State Derivatives: 

WCHM01 uses metal temperature but it is not treated as a state (See METLOO above) . 
This module does not treat any parameters as states and therefore does not calculate 
any state derivatives. 

Nozzle Hot Side Heat Transfer Module 

One nozzle heat transfer module, QNOZOI , was used to configure the TTBE. Following is a 
list of the TTBE schematic names that use the QNOZOI routine, a description of the QNOZOI 
module, and the states of QNOZOI module. 
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Module QN0ZQ1 


Schematic Names: QD0TNZ1 , QD0TNZ2 
Description: 

This module calculates the multiple node heat flowrate through the nozzle wall from the 
combustion gases using a Bartz empirical correlation. 

States and State Derivatives: 

QNOZOI uses metal temperature but it is not treated as a state (See METLOO above) . 

This module does not treat any parameters as states and therefore does not calculate 
any state derivatives. 

Module NOZLOQ 

Schematic Names: NOZL 

Description: 

This module calculates the gross thrust, flow through the nozzle, and the exit mach 
number using isentropic relations. 

States and St ate Derivatives: 

NOZLOO calculates flow but it is not treated as a state. This module does not treat any 
parameters as states and therefore does not calculate any state derivatives. 

TTBE States 

The TTBE simulation has 122 states. The states, a description of the states, the module 
names where the states are differentiated, and the corresponding schematic names are listed for 
the fuel side in Table 4-1, for the oxidizer side in Table 4-2, and for the hot gas side in Table 4-3. 

The thermodynamic states, density and internal energy are difficult parameters to iterate. To 
overcome this difficulty pressure and enthalpy are iterated to solve the density and internal energy 
corrector equations. The thermodynamic states, a description of the thermodynamic states, the 
corresponding iteration parameters, and a description of the corresponding iteration parameters 
are listed for the fuel side in Table 4-4 and for the oxidizer side in Table 4-5. 

TTBE Additional Required Balances 

Fourteen additional balances are required to close the loop on pressures and temperatures to 
achieve a power balance. The iteration parameters and the two balance parameters along with 
their descriptions are listed in Table 4-6. 
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Table 4-1 . Fuel Side States 


State 

Description 

Module 

Name 

Schematic 

Name 

SNFL 

Low-pressure fuel turbopump speed 

ROTROO 

LPFP/LPFT 

RHOVL1 

Density Of Vol. 1 

VOLMOO 

VL1 

UTVL1 

Internal Energy Of Vol. 1 

VOLMOO 

VL1 

WF1 

Flow Rate Through Fuel Line 1 

PIPE00 

FI 

RHOVL2 

Density Of Vol. 2 

VOLMOO 

VL2 

UTVL2 

Internal Energy Of Vol. 2 

VOLMOO 

VL2 

SNFH 

High-Pressure Fuel Turbopump 
Speed 

ROTROO 

HPFP/HPFT 

RHOVL3 

Density of Vol. 3 

VOLM01 

VL3 

UTVL3 

Internal Energy of Vol. 3 

VOLM01 

VL3 

RHOVL4 

Density of Vol. 4 

VO.M01 

VL4 

UTVL4 

Internal Energy of Vol. 4 

VOLMOI 

VL4 

WF2 

Flow Rate Through Fuel Line 2 

PIPE00 

F2 

RHOVL5 

Density of Vol. 5 

VOLMOO 

VL5 

UTVL5 

Internal Energy of Vol 5 

VOLMOO 

VL5 

RHOVL6 

Density of Vol. 6 

NCLV00 

VL6 

UTVL6 

Internal Energy of Vol. 6 

NCLV00 

VL6 

TMMTL1 

Metal Temp, for Qdotnozl 

METL00 

Qdotnozl 

RHOVL7 

Density of Vol. 7 

NCLV00 

VL7 

UTVL7 

Internal Energy of Vol. 7 

NCLV00 

VL7 

TMMTL2 

Metal Temp, for Qdotnoz2 

METL00 

Qdotnoz2 

WF5 

Flow rate through fuel line 5 

PIPE00 

F5 

RHOVL8 

Density of vol. 8 

VOLMOI 

VL8 

UTVL8 

Internal energy of vol. 8 

VOLMOI 

VL8 

RHOVL9 

Density of vol. 9 

VOLMOO 

VL9 

UTVL9 

Internal energy of vol. 9 

VOLM009 

VL9 

WF6 

Flow rate through fuel line 6 

PIPE00 

F6 

WF8 

Flow rate thorough fuel line 8 

PIPE00 

F8 

RHOVLIO 

Density of vol. 10 

VOLMOO 

VL10 

UTVL10 

Internal energy of vol. 10 

VOLMOO 

VL10 

RHOVL1 1 

Density of vol. 11 

NCLV00 

VL1 1 

UTVL11 

Internal energy of vol 1 1 

NCLV00 

VL1 1 

TMMTL3 

Metal temp, for Qdotchml 

METL00 

Qdotchml 

RHOVL12 

Density of vol. 12 

NCLV00 

VL1 2 

UTVL12 

Internal energy of vol. 12 

NCLV00 

VL1 2 

TMMTL4 

metal temp, for Qdotchm2 

METL00 

Qdotchm2 

RHOVL13 

Density of vol. 13 

VOLMOI 

VL13 
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Table 4-1. Fuel Side States (Continued) 




Module 

Schematic 

State 

Description 

Name 

Name 

UTVL13 

Internal energy of vol. 13 

VOLMOI 

VL13 

RHOVL14 

Density of vol. 14 

VOLMOO 

VL14 

UTVL14 

Internal energy of vol. 14 

VOLMOO 

VL14 

RHOVL15 

Density of vol. 15 

VOLMOO 

VL15 

UTVL15 

Internal energy of vol. 15 

VOLMOO 

VL15 

RHOVL16 

Density of vol. 6 

VOLMOI 

VL16 

UTVL16 

Internal energy of vol. 16 

VOLMOI 

VL16 

WF7 

Flow rate through fuel line 7 

PIPEOO 

F7 

PHOPBSF 

Density of PB fuel splitter vol. 

VOLMOI 

PBSF 

UTBSF 

Internal energy of PB fuel splitter vol. 

VOLMOI 

PBSF 
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Table 4-2. Oxidizer Side States 


State 

Description 

Module 

Name 

Schematic 

Name 

WOI 

Flow rate through oxid. line 1 

PIPE03 

01 

RHOVL17 

Density of vol. 17 

VOLMOO 

VL1 7 

UTVL17 

Internal energy of vol. 17 

VOLMOO 

VL1 7 

W02 

Flow rate through oxid. line 2 

PIPE00 

02 

RHOVL18 

Density of vol. 18 

VOLMOO 

VL18 

UTVL18 

Internal energy of vol. 18 

VOLMOO 

VL18 

SNOL 

Low-pressure oxid. turbopump speed 

ROTROO 

LPOP/LPOT 

RHOVL19 

Density of vol. 19 

VOLM01 

VL19 

UTVL19 

Internal energy of vol. 19 

VOLM01 

VL19 

W03 

Flow rate through oxid. line 3 

PIPE00 

03 

RHOVL20 

Density of vol. 20 

VOLM01 

VL20 

UTVL20 

Internal energy of vol. 20 

V0LM01 

VL20 

SNOH 

High-pressure oxid. turbopump 
speed 

ROTROO 

HPOP/HPOT 

RHOVL21 

Density of vol. 21 

VOLM09 

VL21 

UTVL21 

Internal energy of vol. 21 

VOLM01 

VL21 

W05 

Flow rate through oxid. line 5 

PIPE00 

05 

RHOVL22 

Density of vol. 22 

VOLMOO 

VL22 

UTVL22 

Internal energy of vol. 22 

VOLMOO 

VL22 

WPOGO 

Liquid oxid. flow into POGO 

POGOOO 

POGO 

PTPOGO 

Pressure in POGO vol. 

POGOOO 

POGO 

HFRPOGO 

Helium frac. in POGO vol 

POGOOO 

POGO 

RMLPOGO 

Mass of liquid in POGO vol. 

POGOOO 

POGO 

W09 

Flow rate through oxid. line 9 

PIPE00 

09 

RHOVL23 

Density of vol. 23 

VOLMOO 

VL23 

UTVL23 

Internal energy of vol. 23 

VOLMOO 

VLWE 

RHOMOI 

Density of main oxid. inj. vol. 

VOLMOO 

MOI 

UTMOI 

Internal energy of main oxid. inj. vol. 

VOLMOO 

MOI 

WOINJ 

Flow rate through oxid. line OINJ 

PIPE00 

OINJ 

RHOPBSO 

Density of PB oxid. splitter vol. 

VOLM01 

PBSO 

UTPBSO 

Internal energy of PBoxid. splitter vol. 

VOLM01 

PESO 

WOIO 

Flow rate through oxid. line 10 

PIPE00 

010 

won 

Flow rate through oxid. line 11 

PIPE00 

011 

RHOVL24 

Density of vol. 24 

VOLMOO 

VL24 

UTVL24 

Internal energy of vol. 24 

VOLMOO 

VL24 

RHOVL25 

Density of vol. 25 

VOLMOO 

VL25 

UTL25 

Internal energy of vol. 25 

VOLMOO 

VL25 
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Table 4-3. Hot Gas Side States 




Module 

Schematic 

State 

Description 

Name 

Name 


RHOFPBI 

Density of fuel PB injector vol. 

VOLMOO 

FPBI 

UTFPBI 

Internal energy of fuel PB injector vol. 

V0LM00 

FPBI 

RHOOPBI 

Density of oxid. PB injector vol. 

VOLMOO 

OPBI 

UTOPBI 

Internal energy of oxid. PB injector vol. 

VOLMOO 

OPBI 

WFFPB 

Fuel flow rate to fuel PB 

PIPEOO 

FFPB 

WFOPB 

Fuel flow rate to oxidizer PB 

PIPEOO 

FOPB 

WOFPB 

Oxid. flow rate to fuel PB 

PIPEOO 

OFPB 

WOOPB 

Oxid. flow rate to oxid. PB 

PIPEOO 

PPPB 

PTFPRB 

Pressure in fuel PB vol. 

PBRN01 

FPRB 

TTFPRB 

Temperature in fuel PB vol. 

PBRN01 

FPRB 

OFRFPRB 

Oxid. fraction in fuel PB vol. 

PBRN01 

FPRB 

HFRFPRB 

Helium fraction in fuel PB vol. 

PBRN01 

FPRB 

PTOPRB 

Pressure in oxid. PB vol. . 

PBRN01 

OPRB 

TTOPRB 

Temperature in oxid. PB vol. 

PBRN01 

OPRB 

OFROPRB 

Oxid. fraction in oxid. PB vol. 

PBRN01 

OPRB 

HFROPRB 

Helium fraction in oxid. PB vol. 

PBRN01 

OPRB 

PTFTBP 

Pressure in fuel turb. bypass vol. 

VOLM02 

FTBP 

TTFTBP 

Temperature in fuel turb. bypass vol. 

VOLM02 

FTBP 

OFRFTBP 

Oxid. fraction in fuel turb. bypass vol. 

VOLM02 

FTBP 

HFRFTBP 

helium fraction in fuel turb. bypass 
vol. 

V0LM02 

FTBP 

PTOTBP 

Pressure in oxid. turb. bypass vol. 

VOLM02 

OTBP 

TTOTBP 

Temperature in oxid. turb. bypass vol. 

VOLM02 

OTBP 

OFROTBP 

Oxid. fraction in oxid. turb. bypass vol. 

V0LM02 

OTBP 

HFROTBP 

Helium fraction in oxid. turb. bypass 
vol. 

V0LM02 

OTBP 

PTFSF 

Pressure in fuel secondary flow vol. 

VOLM02 

FSF 

TTFSF 

Temperature in fuel secondary flow 
vol. 

V0LM02 

FSF 

OFRFSF 

Oxid. fraction in fuel secondary flow 
vol. 

V0LM02 

FSF 

HFRFSF 

Helium fraction in fuel secondary flow 
vol. 

V0LM02 

FSF 

PTOSF 

Prewssure in oxid. secondary flow vol. 

V0LM02 

OSF 

TTOSF 

Temperature in oxid. secondary flow 
vol. 

VOLM02 

OSF 

OFROSF 

Oxid. fraction in oxid. secondary flow 

VOLM02 

OSF 


vol. 
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Table 4-3. Hot Gas Side States (Continued) 


State 

Description 

Module 

Name 

Schematic 

Name 

HFROSF 

Helium fraction in oxid. secondary flow 
vol. 

V0LM02 

OSF 

PTMFI 

Pressure in main fuel inj. vol. 

VOLM02 

MFI 

TTMFI 

Temperature in main fuel inj. vol. 

VOLM02 

MFI 

OFRMFI 

Oxid. fraction in main fuel inj. vol. 

VOLM02 

MFI 

HFRMFI 

Helium fraction in main fuel inj. vol. 

VOLM02 

MFI 

PTMCHB 

Pressure in main chamber vol. 

MCHB01 

MCHB 

TTMCHB 

Temperature in main chamber vol. 

MCHB01 

MCHB 

OFRMCHB 

Oxid. fraction in main chamber vol. 

MCHB01 

MCHB 

HFRMCHB 

Helium fraction in main chamber vol. 

MCHB01 

MCHB 
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Table 4-4. Fuel Side Thermodynamic States and Their Corresponding Iteration Parameters 


Thermo- 

dynamic 

State 

Thermodynamic 
State Description 


Corresp. 

Iteration 

Parameter 

Corresponding Iteration 
Parameter description 

RHOVL1 

Density 

of vol. 1 


PTVL1 

Pressure in vol. 

1 

UTVL1 

Internal 

energy of vol. 

1 

HTVL1 

Ethalpy of vol. 1 


RHOVL2 

Density 

of vol. 2 


PTVL2 

Pressure in vol. 

2 

UTVL2 

Internal 

enegy of vol. 

2 

HTVL2 

Enthalpy of vol. 

2 

RHOVL3 

Density 

of vol 3 


PTVL3 

Pressure in vol. 

3 

UTVL3 

Internal 

energy of vol. 

3 

HTVL3 

Enthalpy of vol. 

3 

RHOVL4 

Density 

of vol. 4 


PTVL4 

Pressure in vol. 

4 

UTVL4 

Internal 

energy of vol. 

4 

HTVL4 

Enthalpy of vol. 

4 

RHOVL5 

Density 

of vol. 5 


PTVL5 

Pressure in vol. 

5 

UTVL5 

Internal 

energy of vol. 

5 

HTVL5 

Enthalpy of vol. 

5 

RHOVL6 

Density 

of vol. 6 


PTVL6 

Pressure in vol. 

6 

UTVL6 

Internal 

energy of vol. 

6 

HTVL6 

Enthalpy of vol. 

6 

RHOVL7 

Density 

of vol. 7 


PTVL7 

Pressure of vol. 

7 

UTVL7 

Internal 

energy of vol. 

7 

HTVL7 

Enthalpy of vol. 

7 

RHOVL8 

Density 

of vol. 8 


PTVL8 

Pressure of vol. 

8 

UTVL8 

Internal 

energy of vol. 

8 

HTVL8 

Enthalpy of vol. 

8 

RHOVL9 

Density 

of vol. 9 


PTVL9 

Pressure of vol. 

9 

UTVL9 

Internal 

energy of vol. 

9 

HTVL9 

Enthalpy of vol. 

9 

RHOVLIO 

Density 

of vol. 10 


PTVL10 

Pressure of vol. 

10 

UTVL10 

Internal 

energy of vol. 

10 

HTVL10 

Enthalpy of vol. 

10 

RHOVL1 1 

Density 

of vol. 11 


PTVL11 

Pressure of vol. 

11 

UTVL11 

Internal 

energy of vol. 

11 

HTVL11 

Enthalpy of vol. 

11 

RHOVL12 

Density 

of vol. 12 


PTVL12 

Pressure of vol. 

12 

UTVL12 

Internal 

energy of vol. 

12 

HTVL12 

Enthalpy of vol. 

12 

RHOVL13 

Density 

of vol. 13 


PTVL13 

Pressure of vol. 

13 

UTVL13 

Internal 

energy of vol. 

13 

HTVL13 

Enthalpy of vol. 

13 

RHOVL14 

Density 

of vol. 14. 


PTVL14 

Pressure of vol. 

14 

UTVL14 

Internal 

energy of vol. 

14 

HTVL14 

Enthalpy of vol. 

14 

RHOVL15 

Density 

of vol. 15 


PTVL15 

• 

Pressure of vol. 

15 

UTVL15 

Internal 

energy of vol. 

15 

HTVL15 

Enthalpy of vol. 

15 

RHOVL16 

Density 

of vol. 16 


PTVL16 

Pressure of vol. 

16 

UTVL16 

Internal 

energy of vol. 

16 

HTVL1 6 

Enthalpy of vol. 

16 

RHOPBSF 

Density 

vol. 

of PB fuel splitter 

PTPBSF 

Pressure in PB fuel 
splitter vol. 

UTPBSF 

Internal energy of PB fuel 
splitter vol. 

HTPBSF 

Enthalpy of PB fuel 
splitter vol. 

RHOFPBI 

Density of fuel PB injector 
vol. 

PTFPBI 

Pressure in fuel PB 
Injector vol. 

UTFPBI 

Internal energy of fuel PB 
injector vol. 

HTFPBI 

Enthalpy of fuel PB 
Injector vol. 
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Table 4-5. Oxidizer Side Thermodynamic States and Their Corresponding 
Iteration Parameters 


Thermo- 

dynamic 

State 

Thermodynamic 
State Description 


Corresp. 

Iteration 

Parameter 

Corresponding Iteration 
Parameter description 

RHOVL17 

Density of vol. 17 


PTVL17 

Pressure in vol. 17 

UTVL17 

Internal energy of 

vol. 17 

HTVL17 

Enthalpy of vol. 17 

RHOVL18 

Density of vol. 18 


PTVL18 

Pressure in vol. 18 

UTVL18 

Internal energy of 

vol. 18 

HTVL18 

Enthalpy of vol. 18 

RHOVL19 

Density of vol. 19 


PTVL19 

Pressure in vol. 19 

UTVL19 

Internal energy of 

vol. 19 

HTVL19 

Enthalpy of vol. 19 

RHOVL20 

Density of vol. 20 


PTVL20 

Pressure in vol. 20 

UTVL20 

Internal energy of 

vol. 20 

HTVL20 

Enthalpy of vol. 20 

UTVL21 

Density of vol. 21 


PTVL21 

Pressure in vol. 21 

RHOVL21 

Internal energy of 

vol. 21 

HTVL21 

Enthalpy of vol. 21 

UTVL22 

Density of vol. 22 


PTVL22 

Pressure in vol. 22 

RHOVL22 

Internal energy of 

vol. 22 

HTVL22 

Enthalpy of vol. 22 

UTVL23 

Density of vol. 23 


PTVL23 

Pressure in vol. 23 

RHOVL23 

Internal energy of 

vol. 23 

HTVL23 

Enthalpy of vol. 23 

RHOMOI 

Density of main oxid. injector 
vol. 

PTMOI 

Pressure in main oxid 
injector vol. 

UTMOI 

Internal energy of main oxid. 
injector vol. 

HTMOI 

Enthalpy of main oxid 
injector vol. 

RHOPBSO 

Density of PB oxid 
vol. 

. splitter 

PTPBSO 

Pressure in PB oxid. 
splitter vol. 

UTPBSO 

Internal energy of 
splitter vol. 

PB oxid 

HTPBSO 

Enthalpy of PB oxid. 
splitter vol. 

RHOVL24 

Density of vol. 24 


PTVL24 

Pressure in vol. 24 

UTVL24 

Internal energy of 

vol. 24 

HTVL23 

Enthalpy of vol. 24 

RHOVL25 

Density of vol. 25 


PTVL25 

Pressure in vol. 25 

UTVL25 

Internal energy of 

vol. 25 

HTVL23 

Enthalpy of vol. 25 

RHOOPBI 

Density 


PTOPBI 

Pressure in oxid. PB 
injector vol. 

UTOPBI 

Internal energy of 
injector vol. 

oxid. PB 

HTOPBI 

Enthalpy of oxid. PB 
injector vol. 
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Table 4-6. TTBE Model Additional Required Balances 


iterated 

Parameter 

Iterated 

Parameter 

Description 

Balance 
Parameter t 

Balance 
Parameter 1 
Description 

Balance 
Parameter 2 

Balance 
Parameter 2 
Description 

WLPFP 

LPFP flow 

PTVL1 

Vol. 1 pressure 

PTLPFD 

LPFP Dlsch. Pressure 

WHPFP 

HPFP Flow 

PTVL3 

Vol. 3 Pressure 

PTHPFD 

HPFP Dlsch. Pressure 

WLPOP 

LPOP Flow 

PTVL19 

Vol. 19 Pressure 

PTLPOD 

LPOP Dlsch. Pressure 

WHPOP 

HPOP Flow 

PTVL21 

Vol. 21 Pressure 

PTHPOD 

HPOP Dlsch. Pressure 

WPRBP 

PBRP Flow 

PTPBSO 

Vol. PBSO 
Pressure 

PTPBPD 

PRBP Dlsch. Pressure 

TTHTFD 

HPFT Dlsch. 
Temp. 

TTHTFD 

HPFT Dlsch. 
Temp. 

TTHTFDC 

Calc. HPFT Dlsch. Temp. 

TTHTOD 

HPOT Dlsch. 
Temp. 

TTHTOD 

HPOT Dlsch. 
Temp. 

TTHTODC 

Calc. HPOT Dlsch. 
Temp. 

WHG1 

FPRB Dlsch. 
Flow 

PTFPRB 

FPRB Pressure 

PTFPRBC 

Calc. FPRB Pressure 

WHG2 

OPRB Dlsch. 
Flow 

PTOPRB 

OPRB Pressure 

PTOPRBC 

Calc. OPRB Pressure 

WHG5 

FSF to MFI 
Flow 

PTFSF 

FSF Pressure 

PTFSFC 

Calc. FSF Pressure 

WHG6 

OSF to MFI Flow 

PTOSF 

OSF Pressure 

PTOSFC 

Calc. OSF Pressure 

WF1NJ 

FINJ Flow 

PTMFI 

MFI Pressure 

PTMFIC 

Calc. MFI Pressure 

WF9 

FL10 to VL1 1 
Flow 

PTVL10 

Vol. 10 Pressure 

PTVL10C 

Calc. Vol. 10 Pressure 

WF10 

FL1 1 to VL12 

PTVL11 

Vol. 11 Pressure 

PTVL11C 

Calc. Vol. 11 Pressure 


Flow 
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SECTION V 

SYSTEM TESTING AND VERIFICATION 


System simulations were generated and operated to verify the proper operation of the 
ROCETS system. Math model simulations utilized for this testing included: 

• Simple TTBE Model 

• Detailed LOX side Model 

• Detailed TTBE Model 

• Sub-set Model 

Because the models did not use all the same equations and calculations of the SSME DTM 
(Reference 2), the resulting predictions were not expected to reproduce exactly the DTM. 

However, comparisons to the DTM results were used as a guide that the models were functioning 
properly in the ROCETS system. 

5.1 SIMPLE TTBE MODEL TEST 

The simple TTBE Model (as defined in Section 4.1) had 55 state variables and 2 algebraic 
loops. The initial test was obtaining a steady state balance at 100% RPL by driving the state 
derivatives* to zero and closing the algebraic loops (to within a specific tolerance) . This 
demonstrated the capability of the ROCETS system to converge a rocket simulation to a 
steady-state point without running a transient. Transient capability was demonstrated by running 
the simple TTBE simulation with small perturbations of preburner valve areas about the 100% RPL 
point. 

Additional algebraic loops (balances) were placed in the model to set preburner valve 
coefficients at points other than 100% RPL. The fuel preburner valve coefficient was iterated until 
chamber pressure (PTMCHB) was equal to the request: 

PTMCHB ReqU est = PTMCHB 100 RPL X %RPL 

The LOX preburner valve coefficient was iterated until chamber oxidizer fraction (OFRMCHB) 
was equal to a constant value of 0.865. (this is equivalent to a mixture ratio of 6.407). Figure 
5-1 shows main chamber pressure and oxidizer fraction as a function of RPL. A series of 
steady-state points between 60% and 119% RPL were then run with the solver iterating on valve 
coefficients until chamber pressure and LOX fraction were satisfied, this demonstrated the ability 
of the model to use the solver as a means to set a model parameter based on an input constraint. 
Output of the run gives a reference steady-state operating characteristic for the model . and 
provides data for SMITE guess curves. Figure 5-2 shows turbine speeds as a function of RPL. 
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SIMPLE TTBE SIMULATION 
STEADY-STATE CHARACTERISTICS 


MAIN CHAMBER OXIDIZER FRACTION 


1.4 1.2 1.0 0.8 0 . 6 . 0,4 
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Rated Power Level (0-100%) 

Figure 5-1 - Main Chamber Conditions as a Funciton of RPL 
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Figure 5-2 









The steady-state data as a function of RPL was tabularized and used to construct an 
open-loop control with RPL request input and valve data area requests calculated from the table. 
The valve request were put through a first order lag to simulate actuator dynamics. Gross throttle 
transients were run by inputting an RPL request as a function of time and using the open-loop 
control to provide valve areas. Figures 5-3 and 5-4 show results of a transient run from 100% to 
65% RPL decel and figures 5-5 and 5-6 show a transient run from 65% to 109% RPL. These test 
with the simple TTBE model demonstrated the ability of ROCETS to obtain a steady-state balance, 
obtain steady-state valve schedules based on imposed constraints, and to operate transiently. 
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SIMPLE. TTBE SIMULATION 
STEADY-STATE CHARACTERISTICS 
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TIME (SEC) 

Figure 5-3 - Turbine Speeds as a Function of Time 
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Figure 5-4 - Main Chamber Conditions as a Function of Time 
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Figure 5-6 - Main Chamber Conditions as a Fucntion of Time 






5.2 DETAILED LOX SIDE MODEL TEST 


The Oxidizer side of the detailed Technology Test Bed Engine (TTBE) model was configured to 
verify component fidelity and basic model definition. Figure 5-7 shows a schematic of the model 
with the station names labeled. This version of the model did not include a pogo system, which 
was added later. 

In order to verify the configured model, a shutdown transient was executed by giving the 
detailed lox side model transient inputs from the Rocketdyne Digital Transient Model (DTM) of 
reference 2. These inputs were the high-pressure rotor speed, the low-pressure pump inlet 
pressure, the pogo flowrate, and the oxidizer flowrates to the preburners and main chamber. 
Figures 5-8 and 5-9 show these inputs as a function of time. Some results from operating the 
lox-side TTBE simulation are shown in Figures 5-10 to 5-13 with overlays to the DTM predictions. 

A 10ms time step, which is approximately 50 times larger than the DTM, was used with the implicit 
integration scheme of ROCETS. 

With 5 passes, or less, at each time step to obtain the implicit integration, an order of 
magnitude savings exist in computer calculations for a simulated transient. Figure 5-10 & 5-11 
presents the low pressure pump and line 2 flowrates for both simulations showing excellent 
agreement. The comparisons between the two simulations of low pressure rotor speeds (Figure 
5-12) and mixer 2 total pressure (figure 5-13) are also very good transiently with only some initial 
steady-state level differences. 
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Figure 5-7. ROCETS TTBE Model Loxside Schematic 
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Transient Tima (sec) Transient Time (sec) 

Figure 5-8. ROCETS TTBE Loxside Shutdown Transient - Transient Inputs 
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Transient Time (sec) 

Figure 5-11. ROCETS TTBE Loxside Shutdown Transient - Line 2 Flowrate 
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Transient Time (sec) 

Figure 5-12. ROCETS TTBE Loxside Shutdown Transient - Low Pressure Rotor Speed 
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Figure 5-13. ROCETS TTBE Loxside Shutdown Transient - Mixer 2 Total Pressure 



















5.3 DETAILED TTBE MODEL TEST ' 

The complete detailed TTBE model was balanced at 100% RPL and compared to the DTM. As 
shown on Table 5-1 , the balance condition was close to the values of the DTM parameters. 

Various open loop transients in main stage operation were exercised to verify proper operation of 
the model operating in the ROCETS system. Figures 5-14 through 5-21 present the results of one 
of these experiments. The engine was initiated at 100% RPL, and then the fuel preburner valve 
was closed 10% (Figure 5-14). The transient response to selected parameters of flows, pressure, 
speeds, and pogo flowrate and pressure are presented on Figure 5-15 through 5-21. 
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Table 5 -1 . Comparison Between DTM & TTBE 
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Figure 5-1 
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Transinet Time (Sec) 

Figure 5-17 - TTBE Model Response - Fuel and Oxidizer Preburner Pressures 
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Figure 5-19 - TTBE Model Response - Speeds of Oxidizer Pumps 
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Transinet Time (Sec) 

Figure 5-20 - TTBE Model Response - Speeds of Fuel Pumps 
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Figure 5-21 - TTBE Model Response - POGO Flowrate and Pressure 












5.3.1 Shutdown Transient 


The detailed TTBE model exercised a shutdown transient from 100% RPL. The open loop 
valve schedules were taken from a DTM shutdown and imposed on the TTBE model as presented 
on Figures 5-22 and 5-23. Comparison plots of TTBE and DTM predictions of selected model 
parameters are presented on Figures 5-24 thru 5-31. Parameters presented are main chamber 
and preburner pressures and temperatures along with the rotor speeds of the four turbopumps. In 
general the TTBE decelerated faster than the DTM, but no attempt was made to tune the TTBE 
model. The significant verification from the test was the model could operate successfully, 
including implicit integration, through this transient of such drastic operating changes. 
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Figure 5-22 - Shutdown Open-Loop Valve Schedules 
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Figure 5-23 - Shutdown Open-Loop Valve Schedules 


PRATT & WHITNEY - ROCKET PERFORMANCE 



79 


Figure 5-24 - Shutdown Chamber Pressure 
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Figure 5-25 - Shutdown Pressures of Both Preburners 
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Figure 5-26 - Shutdown Fuel Turbopump Speeds 
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TIME (SEC) 

Figure 5-27 - Shutdown Low Fuel Turbopump Speed 
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Figure 5-28 - Shutdown Oxidizer Turbopump Speed 
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Figure 5-29 - Shutdown Low Oxidizer Turbopump Speed 
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Figure 5-30 - Shutdown Main Chamber Temperature 
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Figure 5-31 - Shutdown Preburners Temperatures 


5.3.2 Start Transient 

The start transient was simulated by providing the following changes in the TTBE MODEL. 1 ) 
Test 'ROTROI ’ was created with a minimum break-away torque requirement before allowing pump 
rotation. 2) Heat transfer Q’s representing the latent heat of the nozzle were input as schedules 
of times. 3) The temperatures of the hardware metals were not integrated. 4) Calculations were 
added to simulate the filling, or priming, of the LOX injectors for the two preburners and the main 
chamber. The filling representations were the simple models taken from the DTM, but not the 
detailed, multi-volumes models of the preburners which were used in the DTM predictions. 

The start predictions of the TTBE model are presented on Figures 5-32 to 5-40 with 
predictions of the DTM for reference. The rotor speeds and pressure/temperature of the main 
chamber and preburners are presented. Some of the differences between the predictions of the 
two models is due to the LOX injectors filling times (see Table 5-2). With the earlier fuel 
preburner priming the fuel speed of the TTBE leads the DTM at the 1 .5 sec time (Figure 5-32) . 

The higher chamber pressure results from the higher temperature of the TTBE after ignition 
(Figures 5-36 and 5-38). While other differences exist between the predictions of the two models, 
the verification test was to show the TTBE model in ROCETS could operate through all the transient 
phases. This was successfully accomplished including operation with the implicit integration 
scheme. 


Table 5-2. LOX Injector Priming Times (Sec) 


LOX Injector Priming Times (SEC) 


TTBE 

DTM 

Fuel P/B 

1.22 

1.40 

Main Chamber 

1.55 

1.50 

Oxidizer P/B 

2.00 

1.60 

Combustor Ignition Times (SEC) 


Fuel P/B 

0.45 

0.45 

Oxidizer 

0.90 

0.90 

Main Chamber 

1.45 

1.45 
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Figure 5-32. Start Fuel Turbopump Speed 





ooooe 


00002 


0000 1 


(Wdd) peeds dujndoqjni jsz|p|xo 


89 



PRATT & WHITNEY - ROCKET PERFORMANCE 



(Wdfcl) peeds dwndoqjni |onj moh 


90 


Figure 5-34. Start Low Fuel Turbopump Speed 
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Figure 5-35. Start Low Oxidizer Turbopump Speed 
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Figure 5-36. Start Main Chamber Pressure 


DTM0889: RO START TO 100% RPL 8/16/89 105818R033 
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Figure 5-37. Start Preburner Pressure 
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Figure 5-38. Start Main Chamber Temperature 



DTM0889: RD START TO 100% RPL 8/16/89 105818R033 

TTBE ST TM ON .005 PRPL .01 TOL .001 GEAR1ST TTBE100 


mi 

ii i 

i|i 

•hi 

m 

mi 

BSSSSSSSShhhhhSSii 

L mmmmmaaa Sg amamaam ^^mmurnmmmmm 


III 


miiimaagBBaga 



(d o ) 0 jrnBJQdLU 0 i joujnqejd |enj 


95 


Figure 5-39. Start Fuel Preburner Temperature 
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Figure 5-40. Start Oxidizer Preburner Temperature 



5.3.3 Closed-Loop With Control 


A NASA-MSFC control Model was interfaced in the ROCETS system (Appendix C) and used for 
closed-loop operation with the TTBE Model. Figures 5-41 through 5-44 present results for a 
throttle transient from 100% power to 65% power and back to 100% power. Parameters shown are 
chamber pressure, mixture ratio, and the four rotor speeds. 

5.4 SUB-SET MODEL GENERATION TEST 

To verify the generation of the linear model partials, a linear, sub-set model of the detailed 
TTBE model was created. Then the linear model time domain response was compared to the 
non-linear model predictions. 

The detailed TTBE model had 122 states and 14 algebraic balances. Of the 122 states, 60 
are using iteration parameters other than the states (i.e., 30 volumes are using pressure and 
enthalpy as the iteration variables to close the density and internal energy corrector equations) . 

To reduce the linear model order to a manageable size for the verification test of the new partial 
generation technique, all states were set to be driven to their steady-state values except for the 
four rotor speeds. Thus, of the 136 TTBE simulation equations, 132 were analytically eliminated 
leaving a 4 state model. 

The linear model was generated at 100% RPL with a 0.1% perturbation size. The oxidizer 
preburner oxidizer valve area was used as the model input, and pressure at the low pressure fuel 
pump discharge was the model output. Time domain results were obtained using approximately a 
1 .25% step on valve area by first generating transfer functions from the linear model matrices and 
performing an inverse Laplace transform. 

The non-linear model was executed using the same constraints (i.e., all states forced to 
steady-state except for the four rotor speeds) for comparison to the linear model results. It 
should be noted that a steady-state balance was not performed prior to initiating the time 
transient, so some initial drift is observed. 

Figures 5-45 through 5-48 present comparisons of the linear model to the non-linear model. 
Excellent agreement is observed, especially considering that the time response has an order of 
magnitude larger step than the perturbation size used to generate the partials. 

The excellent agreement is verification of the linear model generation method, it involves a 
change-of-variables for 60 states and analytically eliminating the 14 algebraic balances and the 
118 states which were set to steady-state. The partial generation technique provides a powerful 
tool for performing linear analysis and generation of reduced-order models. 
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Figure 5-41. Chamber Pressure During Throttle Transient 
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Figure 5-42. Mixture Ratio During Throttle Transient 
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Figure 5-43. High Spool Speeds During Throttle Transient 
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Figure 5-44. Low Spool Speeds During Throttle Transient 
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Figure 5-47. Linear and Non-Linear Model Comparison for 
High Pressure Fuel Pump Speed 
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SECTION VI 

CONTRACT END ITEMS 


This Section presents the 23 contract end items which in general provided the overall 
requirements for the simulation system. Each contract end item is presented, followed by the 
program accomplishments. 

1 . The simulation system must have the inherent capability to be applied to current and 
projected liquid rocket engine cycles including: Staged combustion, Expander, Gas 
Generator, and Tri-Propellant. Projected engine applications include: TTBE, STBE, STME, 
and OTVE. 

ROCETS is a very flexible simulation system without any built-in rocket engine cycle or 
configuration. The user can select component modules and designate the interfaces to 
configure virtually any conceivable rocket engine. 

2. The simulation system must simulate start, main stage and shutdown phases of engine 
operation. 

The TTBE model was operated in all three modes as presented in Section 5 - System 
Testing and Verification. 

3. The system must include the capability to operate the plant in both open and closed loop 
control. This requires that provision must be made for attaching a control submodel which 
may be either sampled data or continuous. 

The TTBE model start and shut-down transients were operated open-loop. A NASA-MSFC 
FORTRAN control model (Appendix C) was operated closed-loop with the TTBE model in 
main stage operation as presented in Section 5 - System Testing and Verification. 

4. Methodology must be created to allow representation of various failure modes and off 
nominal operating conditions including random parameter variations within each submodel. 

ROCETS was designed with generic component modules which (based on user input) call a 
designated sub-module which provides the specific component performance characteristic. 
Component failure implies the component performance (map characteristic) changes 
drastically. ROCETS can utilize modules which accept a failure flag to switch from a normal 
operating characteristic to a failed performance characteristic. 

5. The simulation system shall be organized so that multiple levels of detail may be user 
selected for each sub-model where appropriate. This requires that both highly detailed 
simulation modes and “quick and dirty" simulation modes may be selected at user 
discretion. 

ROCETS is designed to configure an engine simulation based on user defined modules which 
leaves the amount of model detail up to the user. During the program, a simple TTBE model 
and a detailed TTBE model were generated and operated. 

6. The simulation system shall be designed with a minimum bandwidth goal of 300 Hz for the 
detailed simulation mode. In some instances it might be appropriate to model higher 
frequency dynamics. 

ROCETS has several features to enhance transient operation. Integration methods include 
trapezoidal and Gear (first and second order) , while other methods can be adapted if 
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required in the future. Implicit (closed-loop) integration is recommended for cost-effective 
computer operation, but open-loop Euler integration is also available. A relative time 
constant is calculated in each module and compared to the simulated time step, to 
automatically select integration or differentiation to be used. Therefore, ROCETS can not only 
operate models up to 300 Hz, it also operates them in an efficient manner. 

7. The simulation system shall be designed so major system components, i.e., the generic 
submodel library, the engine specific data, the generic data, and the simulation experiment 
data, will be separate and distinct. 

The ROCETS system has component modules with generic calculations which call 
sub-modules with specific component performance characteristics and data. Properties are 
called when required through a module for each fluid and sub-modules which contain the 
various property maps. After a simulation has been configured, the user defines the 
simulation experiment with inputs to the run processor. 

8. Input data to define a particular engine shall be defined in terms of design data as opposed 
to model parameter data. Likewise, the empirical data necessary to characterize should be 
defined in terms of industry standard practice. For example, a dynamic model of a sensor 
will usually be given in transfer function form. Turbine performance maps will be 
nondimensionalized. This requires that a design data to model data to model data (sic) 
translator component be developed. 

The ROCETS system was designed with a component-by-component module and 
sub-module performance characteristic concept. This allows the user to build-in conversions 
of design data to model parameter data as required. 

9. A consistent set of nomenclature, model generation coding style, and documentation 
requirements shall be defined and adhered to. This requires that the code must be self 
documenting to the extent possible. 

The ROCETS system software standards are presented in the SDS, P&W FR-20284 
(Reference 4). An example of self-documenting code based on these standards is 
presented in Appendix B. 

10. The simulation system shall be designed so that subset simulations may be readily derived 
from the transient simulation of an engine. These subset simulations include linear operating 
point simulations for controls design, fast operating nonlinear simulations for controls analysis 
and parametrics, and real time simulations for hardware-in-the-loop- testing. 

ROCETS provides the capability to generate linear partial derivatives around transient, or 
steady-state operating points, The matrices of these partials are output by the system for 
use in subset simulations or linear control analysis. Because ROCETS can quickly eliminate' 
states in the non-linear simulations by forcing the derivatives to zero, it can be used to 
develop real-time models which require limited number of states. 

1 1 . The simulation system will provide some method of warning the user when a simulation run 
uses out of range data, such as requesting thermodynamic property routines to extrapolate 
to 6000 psia when the data is good to 5000 psia. It also shall be the user’s option to limit 
the warning and/or utilize it as a stopping condition. 

This was accomplished with good traceability and warnings arranged in different levels of 
severity as discussed in the User’s Manual (Appendix A). 
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12. All generic data and mathematical models utilized in the simulation system shall be 
documented in the code such that the user will know the source of the data and will know 
the limitations and assumptions under which the data was generated and employed in the 
system. Specifically, internal documentation shall include: precise explanation of program 
and subprogram purpose, identification of version data and number, identification and 
description of all inputs and outputs, and identification of all blocks of mathematical 
calculations. 

This was accomplished and can be viewed in the example pump module (Appendix B) and in 
the other system modules and sub-modules of the SDS (Reference 4) . 

13. The simulation system shall be generated in the “Advanced Continuous Simulation 
Language” and in FORTRAN 77 unless an overriding justification can be made for an 
alternate approach. Such justification would be if an alternative were shown to be obviously 
and substantially superior to ACSL, or if a necessary capability were identified which would 
be prohibitive to develop in ACSL. 

The ACSL requirement was eliminated at the Critical Design Review at MSFC on 21 July 
1988, because of the following justification: The ACSL system uses a FORTRAN labeled 
common structure to communicate between the ACSL FORTRAN modules. These common 
statements are built through an internal algorithm and are not structured in a predictable 
format, making user interfacing with ACSL modules very difficult. On the other hand, ACSL 
as a system is not structured to generate large, detailed rocket simulations from user 
supplied FORTRAN modules and operate the simulation in an efficient manner. Therefore, the 
ROCETS system should not be generated in ACSL. 

14. The approach to be taken in mathematical modelling shall always give preference to first 
principals models first, empirical correlations second, and a transfer function approach third. 
For example, it is important to use first principals models of volume filling and gross heat 
transfer when modelling an injector prime. On the other hand, turbomachinery performance 
can be obtained by nondimensional performance maps so that simulation run time may be 
kept reasonable. Likewise, a sensor model need only be in transfer function form since any 
increase in detail would greatly encumber the simulation. 

In general, these guidelines were utilized in generating the modules to represent the TTBE 
model. The module building-block architecture of ROCETS allows component models with 
different levels of detail to be substituted if required for particular application. 

15. In general the detailed mode of simulation should be sufficient to reflect the influence as 
would be measured by performance instrumentation and reflected in aggregate internal 
parameters of the following: design changes, property changes, start phenomena, shutdown 
phenomena, control logic performance, key parameters that limit operation like turbine 
temperature limits, instrumentation performance and location effects, engine performance 
variation, interface condition changes, and purge effects. This list is not all inclusive. The 
detail generally required is that reflected in the SSME DTM. 

The SSME DTM (Reference 2) was used as a guide to provide the amount of detail in the 
TTBE simulation. 

16. The acceptance test of the simulation system shall be a complete simulation of the 
Technology Test Bed Engine. All thermodynamic and thermophysical property data generated 
for the simulation system must reflect the requirements that the TTBE has for such data. 
Likewise, heat transfer correlations must be valid in TTBE operating ranges. To provide 
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capability for the modeling system to be utilized in the study of hydrocarbon engines, 
thermodynamic and thermophysical data must also be supplied for at minimum the most 
likely hydrocarbon propellant candidate. NASA will specify the choice during Phase II efforts. 
These statements require that the data is for characterization of liquid hydrogen, liquid 
oxygen, hydrocarbon fuel, purge gasses, and their materials utilized in the TTBE. 

The TTBE simulation was generated and used to verify the simulation system. Because of 
decreased interest in tri-propellant engines, NASA did release P&W from the hydrocarbon 
requirement at the 14 November 1989 meeting. The P&W system to be delivered will include 
methane thermodynamic properties as part of the property package, but will not include 
combustion properties of methane. The simulation system will accept data tables of 
combustion properties, and NASA can generate the properties in data table format if 
required for tri-propellant simulations in the future. 

17. All typical liquid rocket engine components such as turbines, pumps, valves, ducts, 
accumulators, etc., shall be defined in generic fashion such that they can be connected in 
any user desired manner to simulate any of the engines or engine types listed earlier in this 
document. 

The ROCETS configuration processor allows the flexibility to generate simulations of any 
engine. 

18. To verify proper operation, all normal operating modes of the TTBE will be simulated in both 
the detailed and the quick and dirty modes. In addition, the subset simulation generation 
capability must be exercised. 

As discussed in this report, a simple TTBE model and a detailed TTBE model were generated 
and operated. Linear partials were generated and verified by comparing a linear model 
prediction with the non-linear model prediction in the time domain. 

19. To verify submodel operation, test requirements defined in task must include testing the 
operation of the submodel against known analytical solutions and experimentally verified data, 
when available in open literature. 

The system qualification test plans are written to verify module code by specifying tests to 
be performed and the required evaluation, including comparison source and acceptance 
criteria. As an example, the values in the property tables were compared to National Bureau 
of Standards data. 

20. The simulation system shall be installed and proper operation verified on the MSEC EADS IBM 
3083 computer system. 

This was accomplished. 

21 . At the completion of each sub-model or component, the code must be delivered to NASA 
MSFC for testing and utilization, all submodels and components must be delivered at least 3 
months prior to contract completion in order to assure timely testing. 

The initial software delivery to NASA-MSFC was 27 December 1989, with updates on 5 
March 1990 and 10 August 1990. 

22. A review visit to MSFC will occur on or about six months intervals. A Critical Design Review 
will be performed as a part of the first review, with MSFC concurrence required for work to 
proceed. The results of tasks 1, 2 and 3, in Phase I of the activities shall be delivered as a 
document to be utilized in Critical Design Review. 
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The Critical Design Review was conducted 21 July 1988. Other reviews occurred on 9 
December 1988, 27 July 1989, and 16 May 1990. 

23. The final report will include a section listing the equations utilized with all ROCETS code. 

All of the equations of ROCETS are presented in the SDS, P&W FR-20284 (Reference 4) . 
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SECTION VII 
CONCLUSIONS 


1 . The ROCETS system is a valuable new tool which will save time and money in developing 
and using liquid rocket engine transient simulations. 

2. The implicit integration scheme saves computing calculation time, and has been used 
successfully with the detailed TTBE model in simulating start, main stage, and shut-down 
transients. 

3. The same simulation can be used for steady-state cycle balance as well as transient 
operation. 

4. FORTRAN models developed outside the ROCETS system can easily be interfaced with and 
operate in the ROCETS system. 
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SECTION VIII 
RECOMMENDATIONS 


1. The detailed TTBE model should be enhanced by verification with engine data. 

2. The ROCETS system should be maintained with future changes and enhancements. 

3. Potential ROCETS enhancements include: 

• All-electronic documentation and on-line user assistance 

• Improved linear partial generation technique 
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Apppendix A 
User’s Manual 


The User’s Manual for ROCETS is contained in the SDS, P&W FR-20284 (Reference 4). It is 
reproduced in this report for reference. 


a a_ I ft 


.INKNIKWAIU BUM 


119 PRECEDING PAGE BLANK NOT FILMED 



Pratt t Whitney 
FR - 20284 
31 May 1990 


ROCETS USER'S MANUAL 


31 May 1990 


United Technologies 
Pratt & Whitney 
Government Engine Business 
West Palm Beach. Florida 


120 




Pratt l Whitney 
FR • 20284 
31 May 1990 


Table of Contents 


3.4.4. 1 INTRODUCTION TO ROCETS . 1 

3.4.4. 2 GENERAL DESCRIPTION . 2 

3.4 4.2.1 The Module 2 

3.4.4 2.2 The Sub-Module 2 

3.4.4.2 3 Variable Tags 2 

3.4 4.2.4 Executive Programs 3 

3.4 4.3 CONFIGURATION INPUT 4 

3 4 4.3 1 OPTIONS Keyword 4 

3.4 4.3.2 INSTREAM Keyword 4 

3 4 4 3.3 EXTERNALS Keyword 5 

3 4.4.3.4 INTEGRATION Keyword S 

3 4.4.3 5 BALANCES Keyword I 

3.4 4.3.6 SYSTEM A80VE, INSIDE BELOW Keywords 6 

3 4 436.1 MODULE Sub-Block r 7 

3 4.4 3 6 2 PROPERTY Sub-Block I 

3.4.4.3 8 3 EQUATION Sub-Block 10 

34.4.4 RUN INPUT 11 

3.4.4.4.1 SCHEDULES Keyword 11 

3 4 4.4 2 INPUT Keyword 12 

3.4 4.4.3 INTEGRATION DEFAULTS Keyword 13 

3.4.4.4 4 INTEGRATION EXCEPTIONS Keyword 14 

3.4 4.4.5 BALANCES Keyword 14 

3.4.4.4 6 BALANCE DEFAULTS Keyword 15 

3.4 4.4.7 BALANCE EXCEPTIONS Keyword 15 

3 4. 4.4.6 LINEARIZATION Keyword 16 

3 4.4.4 9 LINEARIZATION DEFAULTS Keyword 16 

3.4.4.4.10 LINEARIZATION EXCEPTIONS Keyword . 17 

3.4.4.4.11 RESTART Keyword 17 

3.4.4.4.12 OUTPUT Keyword 18 

3.4.4 4.13 RUN Keyword 19 

3.4.4.$ BUILDING A MOOULE 21 

3.4.4.5.1 Module Communication . 21 

3 4.4.5 2 Interface Data Section 21 

3 4.4.5 3 Module Print/Plot Output . . . . 25 

34.4.6 SIMULATION DEBUG 26 

34 4.0.1 Error flags 26 

3.4.4,7 RUNNING A MODEL 27 

3.4.4.7.1 Running the Configuration Processor (TSO) 27 

3.4 4.7 2 Running a configured model (TSOI 28 


34.4 

iii 


121 




Krau a wnitnay 
FR • 20284 
31 May 1990 


3.4.4.1 Introduction to ROCETS 

ROCETS is an acronym for ROCket Engine Transient Simulation. The objective of the ROCETS system 
Is to apply the state-of-the-art in modeling and simulation technology to simulating liquid rocket 
engines. The versatility of this system makes it ideal for the performance engineer, the system 
engineer, and the control engineer. Also, the structure of ROCETS makes it highly adaptable to 
simulate any type of rocket engine cycle with varying levels of modeling detail as desired by the user. 

The ROCETS system is designed for use by engineers with average experience. While extensive 
modeling experience is not required, it is assumed that (he user is familiar with modeling practices 
and techniques. The goal of ROCETS is to aid the user in creating and using a simulation by 
automatically generating an executable model from input, scanning the model for undefined variables 
or variables which require algebraic loops, and supplying state-of-the-art numerical techniques A 
flexible run-time processor aids in defining inputs for a particular model experiment. In addition, the 
ROCETS system makes available fully verified engineering representations of most rocket engine 
components. The modules in ROCETS which implement non-linear engineering representations are 
written In structured FORTRAN77. The system also has provisions lo generate linear partial 
derivatives at user selected points for subset models 
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3.4A.2 General Description 

The ROCETS system implements engineering representation in the form of FORTRAN subroutines 
called 'modules* The modules are stored on the ROCETS library and are accessible for generating 
simulations. A configuration processor is used to generate an executable simulation from user input. 
Once a simulation is generated, input is supplied to a run processor to execute a particular simulation 
experiment. 


3.4.4. 2.1 The Module 

A module in the ROCETS system is a stand alone FORTRAN 77 subroutine which implements the 
engineering equations to represent a particular engine component A module is distinct from other 
types of subroutines in that only modules communicate directly with the main program All 
communication between modules is via the main program using named variables 

As part of required user input when defining a particular simulation, each selected module must 
include a character name to distinguish the variables associated with that module from other 
variables in the simulation. The name can be up to lour characters The actual variable names are 
formed by concatenating the module name with pre- defined system names for each type of variable. 
As an example, consider the variable name for density Inside a volume Let the volume module name 
supplied by the user be 'VOL1\ The system name for density is RHO' so that the actual variable 
name is 'RHOVOU*. 

In addition to the system names to be used for the variables of each module a variable 'tag* is 
contained in the comment cards at the beginning of each module The variable lag is used to group 
all the variables comprising the model into several categories depending on their function in the 
simulation. The categories are important because only certain variable types can be used for various 
functions. 


3.4.4. 2.2 The Sub-Module 

Sub-modules are called by modules or other sub-modules by a FORTRAN subroutine call list. They 
are stand alone subroutines but. unlike modules, they do not communicate with the main program. 
Sub-modules are divided into map sub-modules and utility sub-modules. Map sub-modules are 
performance characteristics representing a particular component The user selects which map to use 
for a given component along with the ability to 'scale* the map Utility sub-modules implement 
generalized functions and are typically analytic engineering representations or mathematical 
operations. 


3.4. 4.2.3 Variable Tags 

States are variables for which derivatives are calculated and whose values will normally be obtained 
through numerical integration using a predictor-corrector scheme In addition to stales, there can be 
State Iteration Variables. These are variables used as the independent variables for the iteration to 
close the corrector equations. In particular for rocket applications, it is useful to use pressure and 
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enthalpy as the iteration variables for the density and internal energy states First guesses for the 
states or the state iteration variables must be supplied 

External inputs are variables that are used but never calculated. They are tagged as external inputs 
because they must be supplied externally in some form by the user 

Design variables represent variables that are normally fixed for a given engine cycle. Examples are 
volume sizes, line lengths, etc. They have the same function as external inputs but a distinction is 
made for future system enhancements. 

Outputs are simply variables that. are output of a module and may be used as input to other modules 
downstream. No action concerning outputs is necessary by the user 

Independent Balance Parameters are variables that are used as the independent iteration parameter 
for an algebraic loop. Dependent Balance Parameters are the variables that form the error term. 

J.4.4.2.4 Executive Programs 

Four processors are used in the ROCETS system A configuration processor reads the user 
configuration input, retrieves the specified modules and sub- modules from the ROCETS library 
system and builds the simulation. In setting up the simulation, the processor builds the 
communication structure along with global commons It also builds the main program (subroutine 
ROCETS) with the calls to the engineering modules, and any property calls or equations as specified 
In the configuration input 

An input processor is used to interpret user input specifying parameters to define a particular mode* 
experiment, tt consists of a sel of callable FORTRAN subroutines that read user input, interpret the 
input, can load input variables Into the commons, establish balances, and set necessary flags for 
model execution. 

Execution control is provided by an execution processor W controls looping, print, balancing, and 
linearization. Within the execution processor are calls to the numerical utilities that provide 
steady-state balancing, transient integration, and linearization It provides a centralized location, for 
alt numerical operations so that adding new features to the system is simplified 

Output processing is controlled by an output processor that accepts input as to what parameters are 
to be printed and plotted It has an interface routine for plot information so that it can be used with 
a variety of plotting software simply by changing the interface routine However, linearization output 
is not controlled by the outpul processor, but ralher all necessary information is passed to a separate 
interface routine for linearization output. This feature allows tailoring of the linearization output by 
changing only the interface routine 
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3.4.4.3 Configuration Input 


Configuration Input consists of user commands to build a particular simulation. It consists^ required 
information regarding the algebraic engineering modules to be used and their placement. The sysle 
to be modeled is defined using the following keywords: 

OPTIONS 

INSTREAM 

EXTERNALS 

INTEGRATION 

BALANCES 

SYSTEM ABOVE. INSIDE, or BELOW 

3.4.4.3.1 OPTIONS Keyword 

The options block contains optional input to the processing program and should be located at the top 
of the configuration input file Included are engineering units options, provision for a title provision 
lor specifying a POS for those operating on MVS/TSO. and a cross reference option. The form of the 
OPTIONS block is: 

DEFINE OPTIONS 

UNITS * < ENGLISH or 51 ) 1 . . . 

TITLE i ( The title of the node! to be confined ) i 
PDS * ( Data dictionary file nano ) i 
CROSS * { ON / OFF > ; 

END OPTIONS 

The UNITS keyword is used to specify the engineering units to be used Units must be either ENGLISH 
or SI 

The TITLE keyword is used to specify a 50 character title that wilt be placed at the top of the mam 
FORTRAN program. 

The POS keyword is used to specify the name of the file which contains the data dictionary This We 
should contain the INTERFACE. UNITS, and KEYWORDS blocks for all of the modules that are in the 
ROCETS System 

The CROSS keyword is used to turn the cross reference output from the configuration processor on 
or off The cross reference output contains an alphabetized list of every ^ f 

parameter in the model (with the exception of some global variables like IPRPL and lUPOA it) . ™e 
cross reference output can be very useful in debugging a simulation and in verifying that the 
configuration processor produced the desired simulation. 


3.4.4. 3.2 INSTREAM Keyword 

The define instream block contains a list of modules to be used which are not in the system The list 
includes both the module name and a designation of the file in which the module source code r *f , “** % 
For CMS users the file designation consists of a file name, a file type and a file mode. For MVS/TSO 
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users the file designation is the complete file name The processor will read each file and interpret 
the Interface information. The format is: 

DEFINE INSTREAN 

< Nodule n*m* ) ■ {File designation ) t 

( Module near > > (File designation I ; 

END INSTREAN 

For example. If you wanted to test a new heat transfer module that is named HEATOS and is in a 
FORTRAN file named NEWHEAT on your D-disk Ihe DEFINE INSTREAM block would have the form: 

DEFINE INSTREAN 

HEATOS • NEWHEAT FORTRAN Dr 
END INSTREAN 

For MVS/TSO. the OEFINE INSTREAM block might have Ihe form- 
DEFINE INSTREAN 

HEAT05 i ABCD12J. NEWHEAT. FORTRAN; 

END INSTREAN 

3.4.4.3.3 EXTERNALS Keyword 

The define externals block contains a list of external inputs to the simulation. The variables must be 
separated by a comma and Ihe list must end with a semicolon. These are variables which are used 
but never calculated. The processor requires this information to scan for undefined variables and 
required balances. 

DEFINE EXTERNAIS 

(Variable name list) 

END EXTERNAIS 

For example, to have a tank pressure and enthalpy as inputs to the model the OEFINE EXTERNALS 
block would have Ihe form: 

DEFINE EXTERNAIS 
FTHTNK. 

HTHTNKr 
END EXTERNALS 

3.4A3.4 INTEGRATION Keyword 

The integration block allows a change of iteration variables tor implicit integration With 
predictor-corrector methods, it is not necessary that Ihe state be the iteration parameter. The most 
common example of this for rocket applications is the choice of iteration variables tor thermodynamic 
states. The engineering states are generally densily and internal energy. However, density and 
internal energy are difficult to provide first guesses tor, and more importantly, they are difficult 
parameters to iterate upon Better convergence is achieved by iterating on pressure and enthalpy to 
satisfy the density and internal energy corrector equations. The torm of the INTEGRATION block is: 

DEFINE INTEGRATION 

ITERATE* ( ver ) for { sisie ) l 

END INTEGRATION 

If a change of iteration variables for stales is being used, it is specified by Ihe ITERATE keyword, 
followed by the desired iteration parameter and the state for which it is to be used If there is no 
change of iteration variables then nothing needs to be be specified in the OEFINE INTEGRATION block 
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An example of an integration block for a volume that models the dynamics tor a single constituent fluid 
and uses a change of iteration variables follows: 

DEFINE INTEGRATION 

ITERATE i HTVOll for UTVOll l 
ITERATE i PTV0L1 for RHOVOli ; 

END INTEGRATION 


3.4.4.3.S BALANCES Keyword 

This block is used to define algebraic balances at configuration time. It is normally used tor required 
algebraic balances, however it is often useful to set-up other commonly used balances at 
configuration time. The format is. 

DEFINE BALANCES 

balance ( bolnow# ) » < ver } until C v*r > * { ver } i 


END BALANCES 

The functional form has up to an 8 character name 1o uniquely identify the balance followed by an 
Independent parameter name to be iterated until two dependent variables will be equal A problem 
arises when a module output is used as another module's iripul before the output has been calculated 
and the parameters have the same name tl would be desirable to simply concatenate a tag. a C' for 
example, on to the end of one of the parameter names However, for most systems the maximum 
number of characters a variable can have Is eight, all of which could be used when following the 
ROCETS nomenclature. So a 'C' cannot be concatenated on to the end of the variable name 
Therefore the configuration processor assigns a system defined parameter name (of the form 
SYBL000N) to the second occurrence of the parameter. Once the configuration processor has 
assigned the system name, a balance can be set up to drive the parameter in question to be equal 
to the system defined parameter 

For example, suppose you wanted to calculate pressure upstream or a pipe but that pressure has 
been defined as a state for the volume upstream of the pipe The configuration processor would tag 
the pressure upstream of the pipe as requiring an additional balance and would rename one of the 
pressures with a system defined name To achieve a balance, the independent parameter WPtPE is 
varied until the upstream pressure is equal to the volume pressure A balance can be set up in the 
following form and the model reconfigured 

Balance PUPIAl * HPIPE Until PUP * SYBLOOOl j 

If it is desired to establish a balance to drive a dependent variable to a constant value, a name should 
be assigned to the requested value and then used as the second dependent value. The requested 
value must also be added to the external input list. An optional method would be to set the balance 
up at run time instead of configuration time 


3.4.4.3.6 SYSTEM ABOVE, INSIDE, BELOW Keywords 

The engineering representation for the simulation is specified in the SYSTEM blocks Sub-blocks 
within the SYSTEM blocks are used in specifying engineering modules, thermodynamic properties and 
equations. The SYSTEM ABOVE keyword directs the various sub-blocks within the SYSTEM ABOVE 
block to be placed above the iteration loop Likewise, the SYSTEM INSIDE keyword directs the various 
sub-blocks within the SYSTEM INSIDE block to be placed inside the iteration loop Finally, the SYSTEM 
BELOW keyword directs the various sub-blocks within the SYSTEM BELOW block to be placed below 
the iteration loop. The form of the SYSTEM blocks is: 
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DEFINE SYSTEM (ABOVE* or INSIDE* or BCIOH) 
( Sub-Block ) 

( enLryfsJ ) 

END ( Sub-Block } 

END SYSTEM (ABOVE* or INSIDE* or BEL OH) 

3.4.4.3.6.1 MODULE Sub-Block 


The MODULE sub-block is used to specify the engineering module to represent a given component 
along with necessary information concerning the module The form of the MOOULE sub-block is: 


MODULE > ( Modulo Subroutine Ka»e ) 


NAME 
I/O LIST 
DESIGN VALUES 
MAP 
CRT 

END MODULE 


( Component designation ) ; 

(..Node keyword ) 1 { n+me(s) ) 
( none ) * { value ) * . . . i 
( none ) j 

( 65 character Message ) ; 


The module subroutine name appears after the MODULE keyword and specifies the engineering 
module desired to represent a component of the physical system. The NAME keyword is used to input 
a 4 character alphanumeric component designation that is specific to the particular engine 
component 

Nodal connections for the modules are specified by the I/O LIST keyword Node keywords are pari 
of the Interface cards for each module and are used to specify input and output locations. 

Design values for the component can be input by the DESIGN VALUES keyword This will set the 
default value for component design parameters but they may also be input at run time or can be an 
iteration parameter. 

)f the module uses any external maps, this must be entered with the MAP keyword For readability 
of the final main program, a SS character commenl can be input with the CMT keyword The comment 
wilt be pfaced at a comment card prior to the module call. 

As an example, consider a multi-node-volume (VL3) for single constituent fluids with two intet flows 
from pipe R1 and pipe R2. with corresponding inlet properties from volume VL1 and volume VL2. two 
exit flows to pipe R3 and pipe R4. with corresponding exit properties from volume VL4 and volume 
VL5. and one heat now from COOL A schematic and the configuration input for this volume follows 


COOL 


1 Vll 


si 


VIA I 


HI 


R3 

v 


1 

" — — A” — 

VI 3 

x 

1 




| V12 

.v 



VIS J 

A *■ 

“ 
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HOOULEi VOLMOl 
NAME i VI i 

I/O l 1ST i INtET FlOH * Rl . 

EXIT FlOH « R5 , 
1NIET PROP * VI 1 . 
EXIT PROP * VIA . 
9 DOT * COOL l 

DESIGN VALUES ■ VOl * 508. > 

CNT'i VOLUME THREE; 

END MODULE 


R2 . 
R4 . 
VIE . 
VL5 . 


3.4.4. 3. 6. 2 PROPERTY Sub-Block 

The properly sub-block is used to obtain thermodynamic properties as a function of two other 
thermodynamic properties The form of the property block is: 

PROPERTY PACKAOEi < Fsckeoe I ...... i \ 

LOCATION ( Node Heme I i (dependent) * fl Itndepl, indepZ 1 h 


{dependent) * f< (indepl. indep2 )>J 


END PROPERTY 


The package name appears after the PROPERTY PACKAGE keyword 

Within the block, the keyword LOCATION followed by the node name is used to specify at which 
location the properties are requested The particular properties are obtained by specifying the 
dependent properly type as well as (wo independent properly types 

Currently the packages and corresponding options allowed are: 

H2PR0P - Pere-Hydrogen Properties From Maps 

P*F(RH0#H> - Pressure as a function of density and enthalpy 

RH0*F(P#H) - Density as a function of pressure and enthalpy 

T«F(P#H) - Temperature es a function of pressure and enthalpy 

MP(RHOiU) - Pressure as a function of density and internal energy 

RH0*F(P«U) - Density as a function of pressure and internal enargy 

T*F(P* U) - Temperature es e function of pressure end internal energy 
3*F(H,P) - Entropy es e function of onthalpy and prassura 

H*F($,P) - Enthalpy as a function of entropy and pressure 

CP*F(H,P) - Constant Pressure Specific Heat as a function of 
enthalpy end pressure . . . d 

CV*F(H#P> - Constant Volume Specific Heat as a function of 
enthalpy and pressure 

GAMA*F(H#P> - Gamma as a function of enthalpy and 
HU*F(H,P) - Viscosity as a function of enthalpy and prassure 
IC*F<H,P> - Thermal Conductivity as a function of onthalpy 
end pressure 


02PR0P • Oxygen Properties From Haps 


P*F<RH0,H> 

RH0*F<P,H> 

T*FCP#H) 

P*F(RH0#U) 

RH0*F(P,U) 

T*F<P,U> 

S 2 F(H*P) 

H«F(S,P> 

CP*F<H» P> 

CV*F(H#P) 

GAHA*F(H.P) 

HU-F(H,P) 

*«FCH.P> 


Pressure as a function of dansity and enthalpy 
Densi ty as a function of prassura and enthalpy 
Temperature as a function of pressure and an thalpy 
Pressure es e function of density and intarnal enargy 
Density es • function of pressure and intarnal energy 
Temperature es a function of pressuro and internal energy 
Entropy as a function of enthalpy and pressure 
Enthalpy as a function of entropy and pressure 
Constant Prassure Specific Heat as a function of 
enthalpy and prassura . 

Constant Volume Specific Heat as a function of 
enthalpy and pressure 

Gamma as a function of enthalpy and pressure 
Viscosity as a function of enthalpy and pressure 
Thermal Conductivity as a function of enthalpy 
end pressure 
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HEPROP - Helium Properties From Haps 


P*FC RHO.H) 
RHO*FCP, H> 
T»FCP,H) 
5*FCH,P> 
H*FC3,P> 


- Pressure as a function of density and enthalpy 

- Density as a function of pressure and enthalpy 

- Temperature as a function of pressure and enthalpy 

- Entropy as a function of enthalpy and pressure 

- Enthalpy as a function of entropy and pressure 


N2PR0P - Nitrogen Properties Fro* Haps 


P*FCRH0* M) 
RH0»F(P,H> 
$*F<H, P) 
H*F< 5# P ) 
T-FCP.H) 


- Pressure as a function of density and enthalpy 

* Density as a function of pressure and enthalpy 

- Entropy as a ^function of enthalpy and pressure , 

* Enthalpy as a function of entropy and pressure 

- Temperature as a function of pressure and enthalpy 


HEP HOP * He thane Properties From Haps 


P«FCRH0,H> - Pressure as a function of density and enthalpy 
RH0*FCP,H) - Density as a function of pressure and enthalpy 
T*FCH,P) - Temperature as a function of pressure and enthalpy 

H*FCT,P) - Enthalpy as a function of temperature and pressure 

5*F(HiF) - Entropy as a function of enthalpy and. pressure 
H*F($ f P) - Enthalpy as a function of entropy and pressure 
CP=F(T,P> • Constant Pressure Specific Heat as a function of 
temperature and pressure 

CV*F(T,P) - Constant Volume Specific Heat as a function of 
temperature and prassura 

HGPROP - Ideal Gas H2/02 Combustion Properties From Haps 

CP*F(P,T) - Constant Pressure Specific Heat as e function of 
pressure and tamparatura 

GAMA*F(P,T) - Gamma as a function of pressure and tamparatura 
R*F(P,T) - Gas Constant as a function of prassura and temperatura 
RHO-F(PrT) - Density as a function of prassura and tamparatura 
K=FCP,T> - Thermal conductivity as a function of prassura 
end temperature 

MU S FCP,T1 * Viscosity as a function of pressure and temperature 

Z*F(P,T) - Compressibi li ty Factor as a function of pressure 

and temperature 


Additional inputs for combustion properties are the oxygen fraction (OFR) and helium fraction (HFR). 
However* these are always required inputs to the HGPROP property package and do not need to be 
specified by the user within the PROPERTY PACKAGE block 

Examples: 


Using H2PROP for hydrogen, obtain density and temperature as a function of pressure and 
enthalpy at several locations. 

PROPERTY PACKAGE » H2PR0P 

LOCATION 10 * RHO * FCPT, NT), TT * FCPT, NT) ; 

LOCATION PISF i RHO * FCPT. NT), TT » FCPT. NT) ; 

LOCATION FHCO i RHO * FCPT, NT), TT « FCPT, HT) i 

END PROPERTY 

Using HGPROP. obtain gas constant, specific heal, and specific heat ratio Note that the oxygen 
and helium fractions do not have to be specified since ihey are always required and can 
therefore be included in the call list automatically by the processor. 

PROPERTY PACKAGE * HGPROP 

LOCATION OPRO t RGAS * FCPT, TT), CP * FCPT , TT). 

CAHA * FCPT . TT); 

LOCATION FPRS . RGAS « FCPT, TT). CP * FCPT , TT), 

CAHA » FCPT , TT); 

END PROPERTY 
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A property variable for (he executable code is created from the the property type specified for the 
given location, concatenated with the given location. For the H2PROP example the following variables 
would be created for LOCATION PBSF: PTPBSF. HTPBSF. RHOPBSF. TTPBSF. 

3.4.4.3.B.3 EQUATION Sub-Block 

The equation sub-block is used to enter FORTRAN equations into the simulation. The format is: 
EQUATION » ( Fortran Equation ) ; 

Most standard FORTRAN mathematical symbols and intrinsic functions are allowed. Note also that 
the equation input may be continued on up to one subsequent tine and is closed by a semicolon. 

An example of the use of the EQUATION sub-block follows: 

EQUATION i RHO * PT / RGAS / TT ; 
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3.4.4.4 Run Input 

Run Input consists of user commands (o execule a configured simulation. It consists of required 
information to set inputs, define algebraic loops, specify output, and control execution. The following 
set of keywords accomplish this: 

I. SCHEDULES 

2 INPUTS 

3 INTEGRATION DEFAULTS 

4 INTEGRATION EXCEPTIONS 

5. BALANCES 

6. BALANCE DEFAULTS 

7 BALANCE EXCEPTIONS 

9. LINEARIZATION 

9. LINEARIZATION DEFAULTS 

10. LINEARIZATION EXCEPTIONS 

II. RESTART 

12. OUTPUT 
13 RUN 

If a fine within the run input is to be ignored, this can be accomplished by placing an asterisk {*) in 
column one. 

Debug output for the run input will be generated when then following tine is located on the firsl tine 
starling In the first column of the run input file 

I<D€1UG> 

The blocks within the run input are processed as they are encountered, thus the order in which the 
blocks are arranged is important. 

3.4.4.4.1 SCHEDULES Keyword 

The define schedules block is used to input univariant or bivariant curves representing a functional 
relationship for a model input. For steady-stale schedules a system counter named POINT counts 
points for reading schedules and TIME is available for reading transient schedules 

Schedule dependent parameters can be single precision real or integers but the table itself is single 
precision real Schedule independent parameters must be POINT. TIME, an exlernal input, a state. 
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or an Independent balance parameter. A model output cannot be used since Ibis would require an 
Implied algebraic loop. 

The SCHEDULE block consists of two parts The first is the schedule definition and the second is the 
loading of data into the schedule. The format is: 

DEFINE SCHEDULES , , ^ % 

Schedule i ( naif ) is ( d#p ) * F< ( fnd > ) J 

Set schedule * ( naif ) * C data ) l 

END SCHEDULES 

Schedules have a unique eight character name identifier The functional relation can be either 
univarlant or bivariant. If the schedule is bivariant. two independent parameters are required 
separated by a comma The schedule data are in standard map-reader format The first two numbers 
of the schedule data should be set to zero They are used as pointer storage locations by the table 
reader. The third and fourth numbers of the schedule dala indicate the number of data points 
contained in the schedule for the two independent parameters. If Ihe schedule is univariant, the fourth 
number must be set to zero. If extrapolation of the schedule is desired, the third and/or fourth (if 
bivariant) schedule data points should be made negative The resl of the schedule data consist of a 
list of data separated by commas for the first independent parameter in ascending order, followed by 
a fist of data separated by commas for the second independent parameter in ascending order for a 
bivariant schedule, followed by a list of data separated by commas for the dependent parameter. The 
dependent parameter data is arranged with the dependent data corresponding to Ihe first dala point 
for the first independent parameter and all of the corresponding second independent parameter data 
points, followed by the second data point for the first independent parameter and all of the 
corresponding second independent parameter data points and soon. 

Example: Set up a schedule to vary tank pressure (PTANK) from 200 to 100 over 10 seconds, 
represented as a linear variation with a two point curve. 

DEFINE SCHEDULES 

Schedule i htnkpres is PTANIC « F(TIHE) ; 


Set schedule i htnkpres 
END SCHEDULES 

* 0. t 

0., 
200.* 

0 .. 
10.. 
100. i 

2.* 

0 .* 

If extrapolation is desired the schedule should be set up as follows: 
DEFINE SCHEDULES 

Schedule * htnkpres is PTANK * F(TIHE) ; 


Set schedulei htnkpres 
END SCHEDULES 

* 0. * 

0.* 

200.* 

0.* 

10.* 

100.; 

-2.* 

0 .* 


Example: Set up a schedule of drag coefficient (DRAG) as a function of atlitude (ALT) and vehicle 
mach number (VM) with a bivariant schedule 

DEFINE SCHEDULES 

Schedule « vdreg is DRAG = F(ALT.VH) ; 


Set schedulei vdreg 

9 

0 .* 

0 .* 

3 . * 

3., 

m 



0. * 

10000 .* 

100000 . * 




0 .. 

1 .* 

5., 




0 .* 

.7, 

1.4* 




• .* 

.5* 

1.2* 




• .* 

.3* 

1 . 0 ; 


END SCHEDULES 


3.4.4.4.2 INPUT Keyword 

The define Input block is used to define model input values for a particular run. Generally the inputs 
will have been previously defined as external inputs during configuration. However, there are no 
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restrictions on what may be actually input Inputs from a schedule are specified by entering the 
keyword SCHEDULE and the schedule name The format is: 

DEFINE INPUT 

( vsrneme ) * ( data ) # ' 

( verneme } * schedule ( schedule heme } t 
END INPUT 

The following example shows how to define inputs from a schedule and various model input values 
for a particular run: 

DEFINE INPUT . 

PUN > schedule PI PEA 1C , 

CF * 95.0 . 

DPI * 45.0 ; 

END INPUT 


3.4.4.4.3 INTEGRATION DEFAULTS Keyword 

The integration defaults block set$-up default Integration information at run time. It is generally easier 
to set-up default information which is adequate for most stales and then to override the defaults for 
specific states when necessary. The form is: 

DEFINE INTEGRATION DEFAULTS 
Method * 

Convergence Criteria « 

Tolerance * 

Perturbation * 

Allowed Change * 

State Bias * 

State Noriealizer i 

END INTEGRATION DEFAULTS 

The method keywords allow user selection of the Integration technique within the limits allowed by the 
Integration method selected at configuration. Currently there are no limits, but as new integration 
routines are added limits will be necessary Current methods are: EULER. TRAPEZOIDAL, GEAR 1ST 
and GEAR 2ND WARNING: If the balances are on or the state iteration variables are active when 
using the Euler method an error will occur. 

For predictor-corrector schemes, the corrector equation is iterated to convergence using a 
multi-variable Newton- Raph son method The mulli-variable iteration routine includes internal 
Jacobian scaling to improve convergence with stiff systems. In effect both the rows and columns art 
normalized by the maximum element in the .row and column The error term used for convergence 
lesting can be the actual error or it can be normalized as part as pad of the Jacobian conditioning 

Convergence criteria, tolerance, perturbation, and allowed change apply only to iterative methods. 
The keywords to specify the convergence criteria are A8SOLUTE ERROR or NORMALIZED ERROR. 

Tolerance is defined as how close to zero the error term must be before solution is considered 
converged Experimentation to determine a good tolerance value is usually necessary. 

The perturbation is the amount each independent variable is moved for generating a Jacobian ft is 
specified as a percentage of the biased state, input as a decimal fraction. The allowed change is the 
amount an independent variable Is allowed to change each pass. This is necessary in non-linear 
systems to prevent excessive movement leading to exceeding map bounds, etc The allowed change 
Is also a percentage of the biased stale, input as a decimal fraction 

The state bias is the value that is to be added to the slate to bias the state This is necessary if the 
state Is going to change sign or approach zero during the run The state normalizer is the value that 
the state is to be divided by to normalize the stale for the first point If the state normalizer ts set to 


( keyword ) ; 
( keyword ) t 
( velue } i 
( velue } i 
( velue } i 
{ velue > i 
{ velue ) ; 
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zero the state is normalized by the initial state guess plus the state bias foe the first point. For all 
successive points, the stale is normalized by the previous converged value of the biased state. 

The defaults for the various keywords and values are: Method defaults to TRAPEZOIDAL. 

Convergence Criteria defaults to ABSOLUTE ERROR. Tolerance defaults lo 0001, Perturbation 
defaults to 0.01. Allowed Change defaults to 0 1. State Bias defaults to 00. and State Normalizer 
defaults to 0.0. 

3.4.4.4.4 INTEGRATION EXCEPTIONS Keyword 

This block defines exceptions lo Ihe default integration sel-up With a large number of stales, it is 
convenient to set up defaults which will lake care of most slates and override Ihe defaults for specific 
states. The form and adored items are: 


PCFIWC INTEGRATION EXCEPTIONS 

Activation for J ***** ! 1 

Convergence Criteria for < state ) * 

Tolerance for ( state > * 

State lias for ( state > » 

Independent Bias for ( state ) * 

Perturbation for ( ***** ) i 

Allowed Change for < state ) * 

State Nornalizer for C state ) « 

Independent Normalizer for ( state ) 

END INTEGRATION EXCEPTIONS 


{ on# off# steady-state > I 
("keyword ) * 

( value ) t 
( value ) ; 

( value ) i 
( value ) ; 

{ value ) ; 

< value ) ; 
i ( value ) ; 


The default lor all states is to be active always. However, it is often convenient lo iurn states off at 
various times. The activation keyword has options for when the stale is active 

steady-state * the state is always Iterated to steady-state, thereby removing the 
dynamic effect of the state 
on - the state is active 

off * the state is inactive and held constant 

The independent bias and the independent normatizer refer to the state iteration variables defined In 
the DEFINE INTEGRATION block All of the other items have already been discussed. 

3.4.4.4.S BALANCES Keyword 

This block is used to define algebraic balances at run time. The independent variable can be a model 
Inpat or a module design parameter. The dependent variable can be either model output, a state, or 
a constant The maximum number of balances that can be defined at run time is ten. The form is: 

DEFINE BALANCES . f 

balance C b«ln««# > * ( v#r ) until ( v*r > • ( var J i 

END BALANCES 

For example, suppose you wanted to balance the flow exiting a pump. WPUMP, to a flow that is 
calculated downstream of the pump. WCALC. by varying the speed of Ihe pump SNPUMP. 


DEFINE BALANCES 
BALANCE WBAL 
END BALANCES 


SNPUMP until WPUMP = WCALC i 


The DEFINE BALANCES block can also be used to allow a parameter to be one DT out of phase This 
is accomplished by setting up a balance in the form: vary X until X * XC and then turning the balance 
off in the DEFINE BALANCE EXCEPTIONS block This allows X to be used during the convergence 
attempt and then be updated to XC after convergence has been achieved If you want lo turn the 
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balance off and not run one 0T out of phase, the balance should be turned off and the balance 
description should be rearranged in the form: vary X until XC = X. 


3.4.4.4.S BALANCE DEFAULTS Keyword 

The define balance defaults block is similar to the define Integration defaults section. H is used to 
define parameters for configuration or run time defined balances. 

DEFINE 1ALANCE DEFAULTS 

Convergence Criteria i ( keyword ) f 
Dependent Normal izer « { value ) i 

Independent Normal izer i < value ) J 

Tolaranca * < value ) l 

lies * < value ) ; 

Perturbation * ( value ) i 

Allowed Change * 1 value ) ; 

END BALANCE DEFAULTS 


The value of the normalizers are set according lo the following tables 


steady-state point 

Normal izer 


For the first point in a transient run or a 
Condi tion 

Dependent Normalizar Set to 0. 

Dependent Normalizar Set to 0. 
•nd 2nd Dependent Variable * 0. 

Independent Normalizar Set to 0. 

For the successive points of a transient run: 
Condi tion 

Dependent Normalizar Set to 0. 

or • veluo 

Dependent Normalizar Set to 0. 
or a value and 2nd Dependent 
Variable * 0. 

Independent Normalizar Set to 0. 

or a value 


Initial guess for. the 2nd Dependent 
Variable 

Initial guess for the Independent 
Variable plus the bias 

Initial guess for the Independent 
Variable plus the bias 


Mormalizer 

Previous converged value of the 
2nd Dependent Variable 

Previous converged value of tha 
Independent Variable plus the bias 


Previous converged value of tht 
Independent Variable plus the bias 


The defaults for the various keywords and values are Convergence Criteria defaults to ABSOLUTE 
ERROR, Dependent NormaUzer defaults to 0 0. Independent Normalizar defaults to 0 0. Tolerance 
defaults to 0 0001. Bias defaults to 0 0, Perturbation defaults to 0 001. and Allowed Change defaults to 
0 . 1 . 


3.4.4.4.7 BALANCE EXCEPTIONS Keyword 

The define balance exceptions block is similar to the define integration exceptions section. It is used 
to define exceptions to the defined balance defaults for configuration or run time defined balances 
The default for alt balances is to be active always. 
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( on, oft ) J 
( keyword ) J 
( value ) J 
{ value ) i 
( value ) i 
( value ) i 
{ value ) i 
{ value ) j 


3.4.4.4.S LINEARIZATION Keyword 

This block it used to define linearization options and set-up at run time. Note w |\ ei " J ,near jf , l 1£ 

th« slates are controlled through the INTEGRATION DEFAULTS and EX^PTIONS blocks and the 
balances are controlled through the BALANCE DEFAULTS and EXCEPTIONS blocks 


DEFINE LINEARIZATION 

REPEATABILITY CHECK 
LINEARITY CHECK 
INPUTS 
OUTPUTS 

END LINEARIZATION 


{ value ) 2 
( value ) 2 
{ variable list ) l 
{ variable list ) 2 


Typically, the base point {sleady-state or transient) about which linearization is desired is run first ^and 
then the linearization is performed Multiple linearization runs about different base points can be 
accomplished by setting up multiple pairs of base point and linearization runs m senes 

When linearizing stales and balances which are not turned off. and variables which are defined as 

perturbed the other variables which are to be perturbed are held constant The Pe*"*" 1 ?" 
procedure is to first make a positive perturbation, then make a negative perturbation, then repeal 
positive perturbation To perform the repeatability check, the two positive perturbations are used to 
calculate a percent difference which Is compared to the repeatability check value. If the p 
difference is greater than the repeatability check value then a message is printed that d **cnbes _the 
nonrepeatability To perform the linearity check, the positive and negative P art ' a '* ?*“**?. 
calculate an average partial. Ihen a percent difference for the part.als is calculated which is compared 
to the linearity check value If the percent difference is greater than the linearity check value Ihen a 
message if printed that describes the nonlinearity 

INPUTS are variables that are to be inputs lo the linear model (maximum of 15) 

OUTPUTS are non-stale variables lhat are lo be oulputs from the linear model (maximum of 20} 

The defaults tor the values are: REPEATABILITY CHECK defaults to 001. and LINEARITY CHECK 

defaults lo 0.1. 


3.4. 4. 4.9 LINEARIZATION DEFAULTS Keyword 

The define linearization defaults btock is similar to the define integration d * fau ”* **? too N '* U5 * d 
to define values that control the perturbation size and bias of the linear model inputs. 

DEFINE LINEARIZATION DEFAULTS 
Perturbation i ( value ) J 
lias t ( value ) » 

END LINEARIZATION DEFAULTS 

The defaults for the values are: Perturbation defaults to 0 01. and Bias defaults to 0.0. 
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3.4.4.4.10 LINEARIZATION EXCEPTIONS Keyword 

The define linearisation exceptions block is simitar to the define integration exceptions section. It is 
used to define any exceptions to the LINEARIZATION DEFAULTS values for perturbation and bias of 
the linear model INPUTS. 

DEFINE LINEARIZATION EXCEPTIONS 

Perturbation for ( var } i C value ) j 
lias for { var ) » ( value ) i 


END LINEARIZATION EXCEPTIONS 


The following example shows a typical set-up for the linearization block: 

DEFINE LINEARIZATION 

REPEATABILITY CHECK « 8.1 ; 

LINEARITY CHECK * 8.1 ? 

INPUTS* AREAFPOV. AREA0P0V ; 

OUTPUTS* F0 > 

ENO LINEARIZATION 

DEFINE LINEARIZATION DEFAULTS 
PERTURBATION * 0.001 i 
BIAS * 0.0 ; 

END LINEARIZATION DEFAULTS 

DEFINE LINEARIZATION EXCEPTIONS 

PERTURBATION for AREA0P0V * 0.0005 ; 

BIAS for ARCA0P0V * 0.1 ; 

END LINEARIZATION EXCEPTIONS 


3.4.4.4.11 RESTART Keyword 

The RESTART block is used to specify information necessary for eilher restarting a simulation from 
a previously saved balanced point, or for specifying the time and successive lime increment at which 
a restart file is to be written during a run. or for both restarting and writing restart Tiles It is also used 
1o specify a value that is to be passed into the GUESS routine and to determine if the GUESS routine 
Is to be called. (A blank GUESS routine is generated by the configuration processor and must be 
completed by the user) The form is. 

DEFINE RESTART 
INPUT FILE 
OUTPUT FILE 
BEGIN TIME 
DT 

GUESS 

END RESTART 


The INPUT FILE file designation allows the user to specify the file from which the restart information 
Is to be read. Likewise, the OUTPUT FILE file designation allows the user to specify the file to which 
the restart Information is to be writlen. For CMS users the file designation consists of a file name, a 
file type and a file mode, while for MVS/TSO users the file designation is the complete file name. 
NOTE: The first restart file that is to be written will not write over an existing file of Ihe same name. 
If this Is attempted, execution will be halted by the run time reader and a warning message will be 
issued. 


The location of the RESTART block within the run input file can affect the set-up of the run. If the 
RESTART block is located at Ihe top of the run input file, then the remaining blocks can change the 
set-up of the run However, if the RESTART block is located somewhere else within the run Input file, 
then the run input specified in the preceding blocks could be overridden 


i ( file designation ) i 

i ( file designation ) ; 

i ( tie# ) t 

t ( delta tine ) l 

> ( guess value ) J 
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The BEGIN TIME time specifies the lime at which the first restart file is lo be written The DT delta time 
specifies the time increment at the end of which a restart file is to be written Restart files will 
continue to be written over the previously outputted restart file until the lenQth of the run has been 
completed. If only one restart file is desired, the sum of the BEGIN TIME and the DT must be greater 
than the length of the run. 

The guess value is a R # 4 value that will be passed into the GUESS routine. The user can then utilize 
this value to access different sets of guess data within the guess routine. If a guess value is not 
entered tor GUESS, then the GUESS routine will not be called. 

The defaults for the various keywords and values are - BEGIN TIME defaults to 99 999. 0T defaults to 
0.0. and the default is for the GUESS routine not lo be called. 

3*4.4.4.12 OUTPUT Keyword 

The output block Is used to specify output desired from a simulation run. The formal is: 


DEFINE OUTPUT 

TRANSIENT PRINT « ( on/off > / 

LINEARIZATION PRINT « < on/off ) i 

STEADY-STATE PRINT i ( on/ott ) l 

PRINT t { option ) , < print parameter list ) ; 

PLOT * { option ) , ( plot parameter Jut I ) 

PLOT FILE » C file designation ) ; 

PLOT TITLE * < 4d cherecter title ) i 
ERROR HANDLING for < modnsme modloc ) t 
PRINT LEVEL ■ ( val > , 

DIELEVEL * ( val ) , 

DIECOUNT « ( val > , 

END OUTPUT 

If the TRANSIENT PRINT is on and convergence is not achieved, the last pass of the convergence 
attempt is output. If the TRANSIENT PRINT is off. no convergence information is output. U the 
STEADY-STATE PRINT is on a full print of both Ihe Jacobian evaluation and each convergence attempt 
is provided. If the STEADY-STATE PRINT is off, a shod message is printed that specifies if 
convergence was achieved and how many passes were made. 

If the linearization print Is on. exceptions to Ihe linearity check and'or repeatability check are printed 

A variety of options control print and plot outpul for the PRINT and PLOT keywords. The options art' 

NOPRINT - No print output 

JCOPLOT - No plot output 

DUMP ALL - Output ell occurrence* of oil peremoter* 

DUHPONCC - Output the fir*t occurrence* of ell peremeter* 

AIL - Output ell occurrence* of specified peremetor* 

ONCE - Output first occurrencoa of specified peremeter* 

OMITAIL - Omit ell occurrence* of specifiod peremetor* end output ell 
other* . A . 

0HIT0NCE - Omit ell occurrences of specified peremetor* end output 
the first occurrence of ell other* 

If specified parameter* are necessary, the parameter list follows the option keyword 

The defaults for the various on/off Hags and options are: TRANSIENT PRINT defaults to on. 

LINEARIZATION PRINT defaults to on. STEADY-STATE PRINT defaults to ofl PRINT defaults lo 
OUMPALL. and PLOT defaults lo NOPLOT 

The plot file keyword is used to specify the file designation for the file that is lo contain Ihe plot data. 
For CMS users the file designation consists of a file name, a file type and a file mode, while for 
MVS/TSO users the file designation is Ihe complete file name. 
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The ERROR HANDLING keyword allows selection of a print error level and an error level and couni 
at which to atop execution 

The following example will print out the first occurrence of the specified parameters and produce plot 
output of the first occurrence of all the parameters for a transient run using the CMS file designation 

DEFINE OUTPUT 

TRANSIENT PRINT i ON; 

PRINT * ONCE , TIME, PTMCHB, $N0H, SNFH, PTOPRB, PTFPR1 i 

PLOT * DUKPONCC ; 

PLOT FILE * CAPAFILE BINARY 0; 

PLOT TITLE * TRANSIENT RUN NUMBER ONE; 

DEFINE OUTPUT 

The following example will print out all occurrences of all parameters and and produce plot output of 
the first occurrence of FG for a steady-state run with error handling designations using the MVS/TSO 
file designation. 

DEFINE OUTPUT 

STEADY STATE PRINT t ON; 

PRINT * DUMPALl; 

PLOT * ONCE# POINT, FG; 

PLOT FILE * ABCD123 .CAPAFIl E . UNARY; 

PLOT TITLE * STEADY STATE BALANCE; 

ERROR HANDLING for C0MB02 MCHB 1ST; 

PRINT LEVEL « 3000, 

DIELEVEL * 10000. 

DIEC0UNT » 1, 

END OUTPUT 


3.4.4.4,13 RUN Keyword 

The run block is used to define necessary simulation control inputs for a particular simulation run 
The syntax is : 

DEFINE RUN 

( STEADY STATE, TRANSIENT, or LINEARIZE ) * < options ) ; 

END RUN 

The options for STEADY-STATE are 

POINTS * ( value ) 

MAXPASS * ( value ) 

For POINTS, value is the number of consecutive points to be run A system variable POINT will be 
set to one and incremented by one on each steady-state balance for use in schedules. 

For MAXPASS. value is the maximum number of iteration passes that wifi be made before a 
convergence attempt is halted 

The default value for POINTS is Y and the default value for MAXPASS is SO 
The following Is an example of a steady-state run block: 

DEFINE RUN 

STEADY STATE « POINTS ■ 3, MAXPASS * 35; 

END RUN 
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The following is a list of the TRANSIENT options: 

STOP TIME - ending time for transient operation 

DT - transient time increment 

PRINT TIME - time increment for print 

PLOT TIME - time increment for plot 

MAXPASS - maximum number of convergence passes 

Note that If the PLOT TIME is an even multiple of the model DT t « will be changed to an odd multiple 
to avoid masking numerical instabilities. 

The default value for OT is 0 0001 and the default value for MAXPASS is 20. 

N it is required to start the POINT count or TIME at some value other than one or zero respectively, 
this can be accomplished by setting POINT or TIME to the desired value in the INPUT block. 

The following is an example of a transient run block 

DEFINE RUN 

TRANSIENT > DT > .001* STOP TIME • .01, 

PRINT TIME .001, PLOT TIME * .001, MAXPASS * 50 ; 

END RUN 

Currently there are no options for the LINEARIZE keyword. 

The following is an example of a linearize run block 

DEFINE RUN 

LINEARIZE i ; 

END RUN - 
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3.4A.5 Building a Module 

A FORTRAN subroutine can easily be converted 1o the ROCETS system. The following sections must 
be Included in the system module 

1. Subroutine call list 

2. Interface section 

3. History of the module including author and dated list of revisions 

4. Listing of all subroutines and commons required by the module 

Information required to interface a module into the ROCEtS system will be contained in comment 
cards within the prologue of the module The interface section will be in three pads 


1. Interface information 

XBEGlM INTERFACE ( Module Heme ) 
XEND INTERFACE ( Nodule Heme) 

2 Keyword information 

X BEGIN KEYWORDS ( Nodule Heme) 
XEND KEYWORDS ( Nodule Heme ) 

3. Units information 

XBEGIN UNITS { Module Heme ) 
XEND UNITS < Module Heme) 


3.4.4.5.1 Module Communication 

Modules may only communicate to the ROCETS system through the subroutine call list of the module. 
Commons cannot be used to communicate with the main or other modules However, common blocks 
may be used in certain cases for communication between a module and a sub-module 


3. 4. 4.5. 2 Interface Data Section 

The interface data seclion of the module allows the configuration processor to create the 
communication link with the ROCETS system Specific standards for the interface section of modules 
follow: 

M INTERFACE BLOCK € ” 

Tht interface block relates call list names to system names, defines the slalus of each variable for 
system operation, defines the I/O status of each variable, and the FORTRAN variable type The set-up 
consists of 6 pieces of information for each variable: 

1. Call list name 

2. System name 
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3 System tag 

4. Array status 

5. I/O status 

8 . FORTRAN variable type 
A sample interface block is shown below. 

XBEGIN INTERFACE V01MXX 


1 CALL LIST I 

SYSTEM 

SYSTEM TAG I 

ARRAY 

1 I/O 

VAR 1 

I NAME 1 

NAME | 

1 

STATUS 

1 STATUS 1 

1 TYPE 1 


IPRPL 

IPRPL 

GLOBAL 


IN 0 

IM 

I UP DAT 

IUPDAT 

GLOBAL 


IN 8 

IM 

MODN 

MODN 

NAME 8 


IN 8 

CM 

NODI 

NODI 

NAME 1 


IN • 

CM 

N0D2 

N0D2 

NAME 2 


IN 0 

CM 

NODS 

NODS 

NAME 3 


IN 0 

CM 

N0D4 

N0D4 

NAME A 


IN 0 

CM 

N0D5 

N0D5 

NAME 5 


IN 0 

CM 

VOl 

VOL 

DESIGN 


IN « 

RM 

HTIN 

HT 

VARIABLE 


IN 1 

RM 

HIM 

M 

VARIABLE 


IN 2 

RM 

NOUT 

M 

VARIABLE 

* 

IN 3 

RM 

NTOUT 

MT 

VARIABLE 


IN 4 

RM 

ROOT 

QDOT 

VARIABLE 


IN S 

RM 

RH0V0L 

RHO 

STATE 


IN 0 

RM 

UTV01 

UT 

STATE 


IN 0 

RM 

NTVOL 

HT 

VARIABLE 


IN 0 

RM 

DRDT 

DRDT 

DERIVATIVE 


OUT 0 

RM 

DUDT 

DUDT l 

DERIVATIVE 


OUT 0 

RM 

TAUCR 

TAU1 

TAUC 


OUT 0 

RM 

TAUCU 

TAU2 

TAUC 

1 

OUT 0 

RM 


XCND INTERFACE V01MXX 

The first entry is the name in the module cell list, tt should follow ROCETS naming convention, but this 
is not required for proper operation within the system The call list must include as input 4 character 
user defined names for the module name and any associated nodes 

The second entry is the system name Unlike the call list name a consistent naming convention must 
be used tor proper operation in the system The actual variable name wilt be constructed at 
configuration time by concatenating the system name with the proper module / node name Unless 
consistent nomenclature is used, the constructed variable name wilt not match the names constructed 
lor other modules For parameters that are arrays (except for nodes), the first four characters of the 
call list name of the array must be different than the system name This is required to avoid duplicate 
names for arrays that have the same system name The third entry is the system lag. This informs 
the configuration processor what the parameter is used for so proper action can be taken The current 
system tags are* 

GLOBAL 

A constructed name wilt not be generated ~ the call list name will be the actual name. This is 
used for standard flags which are the same for all modules At present, the flag IPRPL used for 
enabling print and the nag IUP0AT used for initialization are global flags. 

NAME n 

The keyword name followed by an integer value is used for module/node names. The integer 
value specifies the particular node for concatenating names. In all cases the integer value 0 
should be used for the module name If the module has configuration dependent arrays, the 
node name wilt also be an array. 
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DESIGN 

Specifies that the variable is a design value for the module. At present no special action is taken 
with the DESIGN keyword (ie.. the keyword VARIA8LE will also work), but it is included as a 
special tag for future enhancements. 

VARIABLE 

Specifies that the parameter is an input/output and has no special significance to the system. 
STATE 

Specifies that the parameter is a slate variable. 

DERIVATIVE 

Specifies that the parameter is a state derivative. States and derivatives hold special 
significance in that pointers must be constructed to locate the states and derivatives within the 
global commons. NOTE it Is required that when a module has multiple states and derivatives, 
they must be ordered in the call list That is. Ihe first derivative must match the first state, and 
so on. Also, any module that calculates a derivative must have the state in the call list even if 
the state is not required. 

DSC 

Specifies that the parameter is a discrete flag Oiscrete flags are used to "freeze' operation 
aboul a discontinuity 

DSCR 

Specifies that the parameter is a discrete flag request. The request is used to inform the system 
on which side of a discontinuity the module should be operating Some action will be required 
when the discrete and discrele request are different after a converged point 

TAUC 

Specifies that the parameter is a critical time associated with a stale As with derivatives, if a 
module has multiple states, the number of critical times must equal the number of states and 
be ordered in the calf list. However, it is not required that a module output critical times. The 
configuration processor will assign default values for any stales for which the critical time has 
not been defined 

MAP 

Specifies that the parameter is the external name of a map subroutine 

The fourth entry is array status. Tor non-array parameters this field is left blank For parameters that 
arrays, the word ARRAY is entered followed by either an asterisk or an integer number An 
asterisk specifies that the array si zt is configuration dependent and the configuration processor will 
count the number of elements In the array Additionally, tor configuration dependent arrays, Ihe 
processor will put the number of elements in the first location and dimension the array to the number 
of elements plus one. An integer number specifies that the array is not configuration dependent and 
the processor takes no special action other than to dimension the array to the specified value 

The fifth entry Is the I/O status. Each parameter must be tagged as either an input (IN), output (OUT), 
input/output (I/O), or oufput/input (O/l) In addition, each parameter must include an integer number 
corresponding to the named node with which it is associated. 
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The final entry Is the FORTRAN variable type. The types are specified by 

R*4 = Single precision real variable 

1*4 = Integer variable 

R*$ =* Double precision real variable 

X*S *= Single precision complex variable 

C*4 m 4 byte character Longer character strings must be treated as C*4 arrays. 


••• KEYWORD BLOCK — 

Keywords may be defined tor inputting" node information to the configuration processor and are 
required for all configuration dependent arrays This allows the user to use engineering terms instead 
of system terms when specifying the inputs/outputs for a module Keywords can be specified only for 
node names 

Keywords are contained on comment cards within a %BEGIN KEYWOROS / %EN0 KEYWORDS block. 
The format of the keywords is 

{ call list name for node > * ( kayword Ifst ) i 

The keyword list may be one or more keyword phrases separated by commas. A semicolon is used 
to terminate the list. The list for any node may be on more than one line. A sample keyword block 
if shown below: 

X BEGIN KEYWOROS V01HXX 

UPSTREAM PROP, INLET PROP . ‘ _ _ * 

UPSTREAM FLOW, WIN, ENTERING FLOW, INLET FLOW i 

DOWNSTREAM FLOW, WOUT, EXIT FLOW ; 

DOWNSTREAM PROP, EXIT PROP ’ 

ROOT ' 

XEND KEYWORDS VOLMXX 

A standard set of keywords is necessary to avoid confusion and promote consistency and readability 
tor configuration input. A preliminary set has been defined and as the ROCETS system begins to be 
used, user comments on appropriate keywords will be used to refine the keyword list. A preliminary 
set of standard keywords follows: 

Properties: 

• Intel Prop. Upstream Prop 

- Exit Prop. Downstream Prop 

- Fuel Prop, Oxidizer Prop. Helium Prop 

Flows 

• Intel Flow. Upstream Flow 

• Exit Flow. Downstream Flow 

Heat Transfer: 

• Qdot 

- Metal Temp, Tmetat 
Shafts 


N001 
NODE 
NODS 
NO DA 
NODS 


- Torq 

• Shaft. Rotor 
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UNITS BLOCK 


Call list parameter units are contained in a % BEGIN UNITS / %ENO UNITS block Both English and 
St units are required Twelve characters are allotted for each set of units. A sample units block is 
shown below: 

XBEGIN UNITS V01HXX 


CAU LIST 
NAME 


IPRPl 

IUPDAT 

HO ON 

NODI 

NO 02 

NODS 

N0D4 

N005 

HI 

HO 


U 
0 
D 

TAUCR 

TAUCU 


ENGLISH 


D f LESS 

O'LESS 

O'LESS 

O'LESS 

O'LESS 

O'LESS 

O'LESS 

on^Ess 

LBH/S 

LIK/S 

ITU/LBH 

ITU/LBH 

IN**3 

ITU/S 

ITU/LBH 

LBH/JN*«3 

BTU/LBH 

LBH/IN**3/S 

BTU/LBH/S 

S 

s 


O'LESS 

O'LESS 

O'LESS 

O'LESS 

O'LESS 

O'LESS 

O'LESS 

O'LESS 

KG/S 

KG/S 

J/KG 

J/KG 

H*»3 

J/S 

J/KG 

KG/H**3 

J/KG 

KG/Mi*3/S 

J/KG/S 

S 

s 



XENO UNITS VOLHXX 


3.4. 4.5.3 Module Print/Plot Output 

Print and plot output is handled by the utility module PRPL01. Parameter output is controlled through 
user requests in the run (execution) input file but is actually output from within modules and 
sub-modules The struclure of a PRPL01 call within a routine is 


CALL PRPL01 < 


( outnaiit ) 


Where n is the number of parameters to be output for this call If you are outputting an array, this is 
the number of elements of the array that you want to output Outname is the eight character name 
that will represent the parameter in the print and plot output Parmunit is the 12 character string of 
the units of the output parameter, however, parmunit is left blank in the PRPL01 calls that are within 
a module or sub-module The proper units are passed to PRPLOt by the system and then output 
Parmname is the parameter real type variable name as it appears in the module or sub-module This 
call wilt be made every time there is a request for plot or print output. A detailed description of 
PRPL01 usage can be found in the SDS Section 3 2 3. 
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3.4.4.6 Simulation Debug 


FORTRAN subroutines ERCK00. ERCK01/ERCK02 arc called each time there is the possibility for an 
error in modules and sub-modules. This run- time error checking aids the user in pinpointing PO«' b '« 
fatal errors and debugging them ERCKOO is called from the routines, and then in turn calls ERCK01 
and ERCK02. The call list for ERCKOO is as follows 

CALL ERCKOO < IUPDAT. MODHAH, M0D10C, IERC0DE. MSC ) 


Hhtrci 


IUPDAT * Update flag <I»«> 

HODNAM * Calling prograa name <C*0) 

HOD IOC » location in calling prograa <C**J 

IERCODE * Error lovol coda <1*0 

HSO * Error message <C*50) * 


0000 


1000 

- 2999 

3000 

- 4999 

5000 

- 4999 

7000 

- 9999 

10000 



3.4.4.6.1 Error flags 

Tht error coda and MODLOC allow for bolh multiple error checks per module and multiple error 
Igvala. The error code currently Is set between 0 and 10000. with 0 being no error and 10000 being 
the fatal error kill level. The error codes are saved in arrays for print/kill checking A list of error 
codes and their corresponding errors follows: 

No error 

Map Cxtrapola tion 
Input out of Rang# 

Internal Iteration Failures 
Invalid Solution 
Invalid Option <No Default)' 
execution halted immediately by ERCKOO 

Note that if an error level of 9000 to 9999 is encountered, execution is halted after the current pass is 
completed 

For multiple error calls within a module. MOOIOC must be unique lor each call. This can be 
accomplished by using the four character module name concatenated with a string that denotes the 
order of error occurrence within the module, 1st. 2nd. 3rd .... for example 

The error print level may be set by the user at run-time. This allows the suppression of lower level 
errors that may have little or no effect upon the overall solution By using the error checking routine 
wisely the user can ensure the model is fully debugged. 
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3. 4. 4. 7 Running a Mode I 

Both the Configuration Processor and the resultant configured model execute in the MVS batch 
environment. A clist. ROCETS. contained in the ROCETS CLIST library is provided to generate the JCt 
submittal dataset necessary to run either the Configuration Processor or a previously configured 
model. The name of this dataset is 'prefix. ROCETS.TEMP JCL". 


3.4.4. 7.1 Running (he Configuration Processor (TSO) 

The ROCETS clist requires no arguments, ail input is obtained through prompts. No validation is done 
on input items A misspelled dataset name will cause subsequent job failure. ROCETS first tries to 
obtain the TSO logon account number to be inserted in the job card If a valid account number cannot 
be obtained then ROCETS wilt prompt for it The following is a sample dialog ROCETS execution to 
run the configuration processor 

1 DO YOU WANT TO ROUTE OUTPUT TO A DIFFERENT NODE {YES/NO)’* * = > 

Answer NO and prompting will proceed to the next topic. Output wilt be held in the output queue 
where it can be viewed via the IOF option of SPF Answer YES and ROCETS prompts for an 
alternate node and userid 

ENTER DESTINATION NODE. (IE PWAGPOH E092928) * - - > 

2 DO YOU WANT TO RECONFIGURE (YES/NO)? - » * > 

Answer YES to run the Configuration Processor. 

3 ENTER DATASET CONTAINING CONFIG INPUT. NO QUOTES * = * > 

Enter here the complete TSO dataset name containing the configuration input without quotes For 
example. ROCETS OATA(CTTBEOOl). 

4. ENTER OATASET FOR CONFIG FORTRAN OUTPUT, NO QUOTES * » * > 

Enter here Ihe complete TSO data name for the configured model output For example 
userid MYMOOEL FORTRAN. 

At this lime the following message is displayed and the job Is submitted to the MVS batch machine 
for execution 

FILE userid ROCETS.TEMP JCL CONTAINS SUBMITTED JCL 

When the job finishes the dataset specified above for configuration output will contain the configured 
fortran model In addition a dataset, userid GUESS FORTRAN, will contain skeleton fortran for initial 
guesses This dataset must be completed and merged with the model fortran prior running the model. 
A sample guess dataset. ROCETS FORTRAN(GTTBEOOI) is provided. If a guess routine is not 
appended to the model fortran dataset the library copy mentioned above will be used 
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3.4.47.2 Running a configured model (TSO) 

Th« ROCETS clisl require* no argumenl*. all inpul is obtained through prompts No validation is done 
on Input items. A misspelled dataset name will cause subsequent job failure ROCETS first tries to 
obtain the TSO logon account number to be inserted in the job card If a valid account number cannot 
be obtained then ROCETS will prompt for it. The following is a sample dialog ROCETS execution to 
run a previously configured model. 

1 DO YOU WANT TO ROUTE OUTPUT TO A DIFFERENT NODE (YES/NO)? * - * > 

Answer NO and prompting will proceed to Ihe next topic. Output will be held in the output queue 
where it can be viewed via the IQf option of SPF. Answer YES and ROCETS prompts for an 
alternate node and userid. 

ENTER DESTINATION NODE, (IE PWAGPOH.E09292S) * - - > 

2 DO YOU WANT TO RECONFIGURE (YES/NO)? = = = > 

Answer NO to run a previously configured model 

3. ENTER DATASET CONTAINING ROCETS FORTRAN, NO QUOTES = - = > 

Enter here Ihe complete TSO dalaset name containing a configured model without quotes. For 
example userid. MYMODEL FORTRAN 

4 ENTER OATASET CONTAINING ROCETS INPUT. NO QUOTES = * * > 

Enter here the complete TSO data name containing run time input For example 
ROCETS DATA(RTTBEOOI) 

5. ENTER DATASET TO CONTAIN LOAD MODULE. NO QUOTES, PRESS ENTER FOR TEMPORARY 
LOAO 

Press enter, no input, for a temporary load dataset. Enter a complete dataset name with the 
correct DC8 attributes for a toad module library to keep Ihe load module. If a dataset with 
incorrect DC8 attributes is specified the job will fail 

At this lime the following message is displayed and the job is submitted to the MVS batch machine 
tor execution. 

FILE userid ROCETS TEMP JCL CONTAINS SUBMITTED JCL 
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Appendix B 

Example Pump Module 


A FORTRAN code listing of a Pump Module (PUMP01) and the engineering documentation is 
presented as an example, all of the engineering modules, sub-modules, and utilities are contained 
in the SDS, P&W FR-20284 (Reference 4). 
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PAGE 1 


SUBROUTINE PUMPOI 


SUBROUTINE PUMPOI C 


CXXMX 


IPRPL 


I UP DAT 

9 

MAP 

9 

MODN 

' 

1 

NODI 

* 

N0D2 

9 

NODS 

9 

N0D4 

9 

2 

HDD 

» 

HTIN 

9 

PTIN 

9 

RHOIN 

9 

3 

RHOOUT 

9 

SNREF 

9 

GRATIO 

9 

SND 

9 

4 

TRQD 

9 

WD 

9 

WIN 

9 

HTOUT 

9 

5 

PTOUT 

9 

TORQ 

) 





6 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


XBEGIN CLASS PUMPOI 


SUBPROGRAM PUMPOI 


XEND CLASS PUMPOI 


UNCLASSIFIED 

UNITED TECHNOLOGIES CORPORATION 
PRATT A WHITNEY 
WEST PALM BEACH, FLORIDA 


SIDi 


x 

x 

E950 x 
x 

X 

X 

X 

X 


CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXIIXIIXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


XBEGIN PURPOSE PUMPOI 

THIS ROUTINE REPRESENTS A CONSTANT DENSITY PUMP, 
XEND PURPOSE PUMPOI 


x 

X 

X 

X 

X 


cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


C XBEGIN HISTORY PUMPOI 

c 

WRITTEN 


- ADDED CHECK FOR FLOW APPROX. EQUAL 
TO ZERO TO ELIMINATE DIVIDE CHECK 
IN DH CALCULATION. 

- ADDED CHECK FOR TORQ APPROX. EQUAL 
TO ZERO TO ELIMINATE DIVIDE CHECK 
IN ETA CALCULATION. 

XEND HISTORY PUMPOI 


08/23/89 M.H.SABATELLA 
AND J.P.SPINN 
S/27/90 T.F. DENMAN 


A/1 0/90 T.F. DENMAN 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C XBEGIN SCHEMATIC PUMPOI 

c 

c /|\ NODS 

C I 

c 

c » 

c « 

c 
c 

c *** 

C NODI > I 

C > | I N0D4 

C N0D2 > I 

C *** 

c * 


x 

x 

x 

x 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


7 

8 

9 

10 

u 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 

0 

1 

Z 

3 

4 

5 

6 

7 

8 
9 
0 
1 
2 

3 

4 

5 

6 

7 

8 
9 
0 
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PRATT A WHITNEY 


PACE 2 


SUBROUTINE PUMP 01 


c 


« 1 

V 

51 

c 


• 1 

X 

52 

c 


• 1 

X 

53 

c 


— I 

X 

5A 

c 



X 

55 

c 

XEND SCHEMATIC PUMP01 


X 

54 

CXXXXXXXXXXXIXXXXXKXKXXXXXXKIXXXXXXXXXXXMXXXXXXKXXXXXXXXXXXXXXXXXXXXXXX 

57 

c 

XBEGIN DESCRIPTION PUMP 01 


X 

5B 

c 



X 

59 

c 

INPUTS b 


X 

40 

c 

IPRPL - 

PRPL01 OUTPUT FLAG 

X 

41 

c 


0 * NO PRINT 

X 

42 

*c 


1 ■ PRINT 

X 

43 

c 

I UP OAT - 

UPDATE FLAG 

X 

4A 

c 


-I * INITIALIZATION'S BALANCE 

X 

45 

c 


0 * TRANSIENT ITERATION PASS 

X 

44 

c 


1 * TRANSIENT UPDATE PASS 

X 

47 

c 

MAP 

EXTERNAL PUMP CHAR. MAP 

X 

41 

c 

MOON 

MODULE NAME (A CHARACTERS) 

X 

49 

c 

NODI 

FLOW NODE (A CHARACTERS) 

X 

70 

c 

NO 02 

INLET THERMAL NODE (A CHAR.) 

X 

71 

c 

NODS 

EXIT THERMAL NODe (A CHAR.) 

X 

72 

c 

NO DA 

SHAFT NODE (A CHARACTERS) 

X 

73 

c 

NOD 

PUMP DESIGN HEAD 

X 

7A 

c 

NT IN 

INLET ENTHALPY 

X 

75 

c 

MIN 

INLET FLOM 

X 

74 

c 

PTIN - 

INLET PRESSURE 

X 

77 

c 

RHOIN - 

INLET DENSITY 

X 

74 

c 

RHOOUT - 

EXIT DENSITY 

X 

79 

c 

SMREF - 

SPEED OF THE SHAFT 

X 

40 

c 

SN 

SPEED OF THE PUMP 

X 

41 

c 

SND 

PUMP DESIGN SPEED 

X 

42 

c 

TROD - 

PUMP DESIGN TORQUE 

X 

43 

c 

MO 

PUMP DESIGN FLOH 

X 

4A 

c 



X 

45 

c 

OUTPUTS! 


X 

44 

c 

NT OUT - 

EXIT ENTHALPY 

X 

47 

c 

PTOUT - 

EXIT PRESSURE 

X 

44 

c 

TORN - 

TORQUE REQUIRED 

X 

49 

c 



X 

90 

C 

INPUTS FROM GUNITS COMMON i 

K 

91 

c 

GC 

UNITS CONVERSION FACTOR 

X 

92 

c 

Gt 

GRAVITATIONAL CONSTANT 

X 

93 

c 

RJ 

PROPORTIONALITY FACTOR J 

X 

9A 

c 



X 

95 

c 

XEND DESCRIPTION PUMP01 


X 

94 

CXXXXXXXXXIXXXXXXKXXXKXKXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

97 

c 

XBEGIN DERIVATION PUMP01 


X 

94 

c 



X 

99 

c 

DERIVATION OF RUHR DISCHARGE PRESSURE CALCULATION 

X 

100 
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PAGE 3 


SUBROUTINE PUHPOI 


FOR A CONSTANT DENSITY PUMP, GIVEN THE FLUID DENSITY. 
ROTATIONAL SPEED. AND FLOW THROUGH THE PUHP. THE 
PUHP HEAD RISE AND REQUIRED TORQUE CAN BE DETERHINED 
FROM THE PUHP HAP. 

FROH CONSERVATION OF ENERGY, THE DISCHARGE PRESSURE 
CAN BE CALCUALTED AS> 


I 0 I 

I POUT » HEAD * RHO * ( — ) ♦ PIN | 
| GC I 


DERIVATION OF EFFICIENCY 


DEFINE EFFICIENCY. ETA 


WORK DONE ON THE FLUID 


ETA ■ 


ENERGY AVAILABLE 


CALCULATE ANGULAR VELOCITY, OHEGA 
2 I PI 

OHEGA * N » < ) 

(I 

WHERE N IS THE ROTATIONAL SPEED IN RPN 


| HEAD * W * (G/GC) I 

| ETA * I 

| TORQUE i OHEGA I 


DERIVATION OF DISCHARGE ENTHALPY CALCULATIQN 

POWER CAN BE DEFINED IN TERHS OF TORQUE AND OHEGA AS* 


TORQUE l OMEGA 


POWER 


POWER CAN ALSO BE DEFINED IN TERMS OF THE CHANCE IN ENTHALPY 
DH AND THE FLOWRATE N AS> 
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PRATT A HHITNEY 


PAGE A 


SUBROUTINE PUMPOl 


POWER » W * DH 


EQUATING (1) AND (2) AND SOLVING TOR DH YIELDS* 
TORQUE ■ OMEGA 


H * RJ 


GIVEN THE INLET ENTHALPY, THE EXIT ENTHALPY IS* 


| HOUT « HIN ♦ DH I 


C XEND DERIVATION PUMPOl 

CX********************™IMW!IMMinni**XXX*K*XX niwncwjnwil XXXXX! x k XXXVXXXX 

C XBEGIN COMMENTS PUMPOl 

C 

c 

C 1. THE EXTERNAL MAP RETURNS EXIT PRESSURE AND TORQUE 

C AS A FUNCTION OF FLOH, EXIT DENSITY, AND SPEED. 

C 

C 2. FOR DESCRIPTION OF COMMON "QUNITS", SEE SUBROUTINE 

C •UNITOO*. 

C 

c 5. THE GEAR RATIO IS DEFINED AS THE SHAFT SPEED DIVIDED 

C BY THE PUMP SPEED. 

C 

C A. THE SIGN OF THE TORQUE FOR A PUMP IS NEGATIVE BY 

C CONVENTION FOR PROPER INTERFACING WITH ROTROO. 

C 

C XEND COMMENTS PUMPOl 

CXXXXXHHXXXXIXXXIXXXXXXXIXIXXXXVXXXKyxXXXXXXKXXMXXUXXXXXXKXXXXMXXKK 

C XBEGIN INTERFACE PUMPOl 
C 


Cl CALL LISTI 
C| NAME I 


SYSTEM 

NAME 


SYSTEM TAG 


1 

ARRAY 

1 I/O 

I VAR 


1 

STATUS 

| STATUS 

I TYPE 

l« 

1 


1 IN 0 

1 I*A 

l« 

1 


1 IN 0 

1 1M 

l« 

1 


1 IN 0 

1 

1* 

1 


1 IN 0 

1 CM 

»« 

1 


I IN 0 

I CM 

1" 

1 


1 IN 0 

I CM 

1* 

1 


1 IN 0 

1 CM 

l« 

1 


I IN 0 

I CM 

1* 

1 


I IN 0 

I RM 

I* 

1 


1 IN 2 

I RM 

1* 

1 


1 IN 2 

1 RM 

1* 


C| IPRPL 
C| IUPDAT 
Cl MAP 
Cl HODN 
Cl NODI 
Cl NODE 
Cl NODS 
Cl NODA 
Cl HDD 
Cl HTIN 


IPRPL 

IUPDAT 

MAP 

MODN 

NODI 

NODE 

NODS 

NODA 

HDD 

HT 

PT 


I GLOBAL 
I GLOBAL 
I MAP 
| NAME 0 
I NAME 1 
I NAME 2 
I NAME S 
| NAME A 
I DESIGN 
I VARIABLE 
| VARIABLE 



155 




PRATT « WHITNEY 


PAGE 5 


SUBROUTINE PUHPOl 


Cl RHOIN 
Cl RHOOUT 
C| SNREF 
Cl GRATIO 
C| SND 
Cl TRQD 
C| HD 
C| HIN 
Cl MTOUT 
C| PTOUT 
Ct TORQ 


| VARIABLE 
I VARIABLE 
| VARIABLE 
I DESIGN 
I DESIGN 
| DESIGN 
| DESIGN 
| VARIABLE 
I VARIABLE 
| VARIABLE 
| VARIABLE 


IN Z 
IN S 
IN A 
IN I 
IN • 
IN 0 
IN 0 
IN 1 
OUT 3 
OUT 3 
OUT D 


C XEND INTERFACE PUMPOI 


C XBEGIN UNITS PUHPOl 
C 


I CALI LIST! OWLISH 

| NAME 1 

1 si 

1 

i 

i 

| IPRPL 

I DHESS 

I DHESS 

i 

| I UP DAT 

| DHESS 

I DHESS 

i 

1 MAP 

| D V LESS 

1 DHESS 

i 

| MODN 

I DHESS 

I DHESS 

i 

1 NODI 

| DHESS 

1 DHESS 

i 

| N0D2 

| DHESS 

1 DHESS - 

i 

| NODS 

1 DHESS 

1 DHESS 

i 

| NO DA 

1 DHESS 

1 DHESS 

i 

1 HDD 

1 IN 

1 M 

i 

1 HTIN 

I BTU/IBM 

I J/KG 

i 

I PTIN 

| LBF/IN«2 

| N/MXK2 

i 

1 RHOIN 

| LBM/INOZ 

| KG/HXXJ 

i 

1 RHOOUT 

| LBft'INXiZ 

| KG/MXX3 

i 

1 SNREF 

1 RPM 

I RPM 

i 

| GRATIO 

| DHESS 

1 DHESS 

i 

1 SND 

1 RPM 

I RPM 

i 

1 TRQD 

I IM-LBF 

| N-M 

i 

1 HD 

1 LBM/S 

1 KG/S 

i 

I HIN 

1 LBH/S 

I KG/S 

i 

1 NT OUT 

| ITU/ LIN 

| J/KG 

i 

I PTOUT 

| LBF/INXKZ 

I N/MXX2 

i 

I TORQ 

I IN-LBF 

I N-N 

i 


C XEND UNITS PUHPOl 

C XBEGIN KEYWORDS PUHPOl 
C 

C NODI i UPSTREAM FLOW# INLET FLOH 
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PACE « 


SUBROUTINE PUMP 01 


C NODE i UPSTREAM PROP, INLET PROP * 1 

C NODS « DOWNSTREAM PROP, EXIT PROP * 1 

C NO DA i SHAFT, ROTOR * 3 

C ! 

C XEND KEYWORDS PUMP01 

CMIHXlIXHHimiMIXnHKMlinnHHHWHIMHHininHHIHmHIl] 

C XBECIN SUBROUTINES REQUIRED PUMP01 1 

C ' 

C SUBROUTINES REQUIRED i PMAPXX MAP (EXTERNAL) 1 

C PRPL01 1 

c ! 

C XEND SUBROUTINES REQUIRED PUHP01 

CXHXXXXXUKKnMMIKHHHIIWHWIHHJtimimilWIIHinHmim*] 

C XBEGIN COMMONS REQUIRED PUMPCI ] 

C 1 

C COMMONS REQUIRED t GUNITS 
C 

C XEND COMMONS REQUIRED PUMPOl 

C***xxx***»*****K********K MM ************** 1KI ***** IK ** ,IKMM * 1,l * R1<I3l,nni11 


CHARACTER** MOON, NODI# N0D2# NODS# NOD* 

EXTERNAL MAP 

COMMON / GUNITS / IUNIT , GC # GR # RJ 

X CIEN , CMASS # CFORCE # CTENP 

X FLOCON 

DATA PI / 5.1*159265* / 


RU 

CENERGY# 


CXXXXVXXXXHXXXXXnXXXlXHHXXXHXI 

C MISCELLANEOUS INITIALIZATIONS « 


CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXII 

SN * SNREF / GRATIO 
SNRAD * <SN*2.*PI/60.) 


CXXXXXXXXXXXXXXXXXXXXXXIXXXUXXXXXXXnXXXXUXXXXIXXXXIIHI 

C READ MAP WITH FLOM AND SPEED FOR HEAD RISC AND TORQUE B 

CHXXxxxxxixxxxxxxxxxxxxxxxxxxxxixixnxxmxixxxxxxxmxw 

CALL MAP C IUPDAT# HDD # RHOOUT# SN # SND . TRQD . 
$ HD # HIN # HD # TORQ ) 

CXXXXXXXXXXXXXXXKXXXXXXXXXXXXIXXXXXXXXXXXHmXXXXX 

C CALCULATE DISCHARGE PRESSURE AND EXIT ENTHALPY * 

CXXXXXXXXXXXXXVXXXXXXXUXXHXHXXXXXXXIXXXIXXXXXIXX 

PTOUT * HD*RHOOUT*GR/GC ♦ PTIN 
IFC TORQ .GT. .01 ) THEN 

ETA * HIN*HD*GR/(SNRAD*TORQ*GC) 

ELSE 

ETA « 0.0 

ENDIF 

POHR * TORQ»SNRAD/RJ 
IF( MIN .QT. .01 ) THEN 


DH 

* POMR/MIN 

ELSE 


DH 

> 0.0 

ENDIF 



251 

252 
255 
25* 

255 

256 

257 
25S 

259 

260 
261 
262 
265 
26* 

265 

266 
267 
261 

269 

270 

271 

272 
275 
27* 

275 

276 

277 
27S 
27* 
2S0 
261 
282 

285 
26* 
265 

286 
267 
288 

289 

290 

291 

292 
295 
29* 

295 

296 

297 

298 

299 
500 
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PAGE 7 


SUBROUTINE PUHPOl 


HTOUT 

• HTIN ♦ 

DH 



SOI 

TORO 

■ - TORQ 




502 

QXXXXXXXXXXXXXXXXXXX 




503 

C PURPLE SECTION x 




304 

Cxxxxxxxxxxxxxxxxxxx 

• 



305 

IF (IPRPL .OT. •) THEN 



304 

CAU 

PRPLOK-9 

,' MODULE SHODM//' OUTPUT 

•#DUHNY) 


307 

CALL 

PRPLOKl# 

• HD '/V HO DM//' 

• #WD 

) 

300 

CALL 

PRPLOKl# 

•TRQDV/MODM ,' 

1 ,TRQD 

) 

309 

CALL 

PRPLOKl# 

'HDO'^/HODM//' S' 

1 # HDD 

) 

310 

« CALL 

PRPLOKl# 

'ETA'/ZMODK//' S' 

•#ETA 

) 

311 

CALL 

PRPLOKl# 

•SM'//NOD«//' S' 

* « SNREF 

) 

312 

CAU 

PRPLOKl# 

•SHV/MODM//' S' 

SSN 

) 

313 

CALL 

PRPLOK 1 # 

•HD'//MODM//' S' 

1 # HD 

) 

314 

CALL 

PRPLOKl# 

•HT'//M0D2//' S' 

9 f HTIN 

) 

315 

CALL 

PRPLOKl# 

•HT'/zNODW S' 

•# HTOUT 

) 

314 

CALL 

PRPLOKl# 

•FT '//MODS//' S' 

•#PTOUT 

) 

317 

CALL 

PRPLOK 1» 

•TORQ'//MODN .' 

9 # TORQ 

) 

310 

CALL 

PRPLOK I, 

'RHO'/vMODZ//' S' 

9 » RHOIN 

) 

319 

CALL 

PRPLOK If 

•RMOV/NODW S' 

SRHOOUT) 

320 

CALL 

PRPLOK If 

•HV/MOD1//' S' 

1 f WIN 

) 

321 

CALL 

PRPLOKl# 

•DH'//MODM//» S' 

*f DH 

) 

322 

CALL 

PRPLOK If 

•POMRV/MODM .' 

•fPOWR 

) 

325 

CALL 

PRPLOKl# 

•SHD'/vMODM/v' S' 

•fSMD 

) 

324 

CALL 

PRPLOKlf 

•PT'//N0D2//' S' 

• #PTIN 

) 

325 

ENDIF 

. 




524 

99 CONTINUE 




327 

RETURN 





320 

END 
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PUMP01 

Thie routine represents a constant density rump. 


/ In NODS 


MODI > 

> 

NO 02 > 


INPUTS i 


OUTPUTS » 


IUPDAT 


HAP 

HOON 

NODI 

N0D2 

NODS 

NO DA 

HDD 

HTIN 

WIN 

PTIN 

UNO IN 

RHOOUT 

SNREF 

SN 

SND 

TROD 

HD 


NT OUT 
PTOUT 
TORO 


PRPL01 OUTPUT FLAG 

0 * NO PRINT 

1 * PRINT 

-i P ? A initiaiization/ss BALANCE 
0* * TRANSIENT ITERATION PASS 
1 * TRANSIENT UPOATE PASS 
external PUMP CHAR. MAP 
MODULE NAME ( A CHARACTERS) 
FLOM NODE (A CHARACTERS) 

INLET THERMAL NODE < A CHAR.) 
EXIT THERMAL NODE < A CHAR.) 
SHAFT NODE (A CHARACTERS) 

PUMP DESIGN HEAD 
INLET ENTHALPY 
INLET FlOH 
I HIE" PRESSURE 
INLET DENSITY 
EXIT DENSITY 
SPEED OF THE SHAFT 
SPEED OF THE PUMP 
PUMP DESIGN SPEED • 

PUMP DESIGN TORQUE 
PUMP DESIGN FLOW 


EXIT ENTHALPY 
EX1 1 PRESSURE 
TORQUE REQUIRED 


INPUTS FROM GUNITS COMMON. ||# rAC70R 

OR - GRAVITATIONAL CONSTANT 

r.I PROPORTIONALITY FACTOR J 


COMMENT Sj. 


The external »np returns exit pressure end tereue 
•s a function of flow. <*xit density, and speed. 


3.4.24 

PUMP01 

1 
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2. For description of common "GUNITS", see subroutine 
"UHIT00". 

3. The gear ratio is defined as the shaft spaod dividod 
by tha pump spaed. 

A. Tho sign of the torque for a pupp is nogativa by 
eonvantion for proper intarfacing with R0TR00. 


KEYHflRPS-* 

No da keywords ara part of tha intarfaca cards for aach module. 

In tha configuration input for a modulo, an 1/0 list containing 
tha noda kaywords is used to spacify tha nodal connections . Tha 
node kaywords for this module are* 

NODI i UPSTREAM FLOW, INLET FLOW 

MOD? i UPStREAM PROP, INLET PROP 

MODS * DOWNSTREAM PROP, EXIT PROP 

NODA * SHAFT, ROTOR 


aEBi y .A.TIQ.MSj. 

A derivation of th# calculations usad in this module follows* 


DERIVATION OF PUMP DISCHARGE PRESSURE CALCULATION 

FOR A CONSTANT DENSITY PUMP . GIVEN THE FLUID DENSITY, 
ROTATIONAL SPEED, AND FLOW THROUGH THE PUMP, THE 
PUMP HEAD RISE AND REQUIRED TORQUE CAN IE DETERMINED 
FROM THE PUMP MAP. 

FROM CONSERVATION OF ENERGY. THE DISCHARGE PRESSURE 
CAN IE CALCUALTED AS » 


! 0 I 

I POUT * HEAD » RHO « (•-) * PIN | 
I GC I 


DERIVATION OF EFFICIENCY 

DEFINE EFFICIENCY. ETA 

WORK OONE ON THE FLUID 

ETA * - 

ENERGY AVAIL All E 

CALCULATE ANGULAR VELOCITY, OMEGA 
2 * PI 

OMEGA * N r ( ) 

A Cl 

WHERE N IS THE ROTATIONAL SPEED IN RPM 


HEAD « H * (G/GC) 

ETA * - 

TORQUE v OMEGA 


I 

I 

I 


DERIVATION OF DISCHARGE ENTHALPV CALCULATION 


34.2.4 

PUMPOt 

2 
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POWER CAN IE DEFINED IN TERMS OF TORQUE AND OMEGA AS> 


POWER « 


TORQUE * OMEGA 


POWER CAM ALSO IE DEFINED IN TERMS OF THE CHANGE IN ENTHALPY 
DH AND THE FLOWRATE M AS* 

POWER * H i DH (2 

EQUATING (I) AND (2) AND SOLVING FOR DH YIELDS* 

TORQUE « OMEGA 

DH 

W * *J 

GIVEN THE INLET ENTHALPY, THE EXIT ENTHALPY IS* 


I HOUT * HIN ♦ DH | 


The configuration processor use* the following intorfaco cards 
to generate the main program ^all list. 


I IPRPL 
I IUPDAT 
I MAP 



IPRPL 

IUPDAT 

MAP 

MODN 

NODI 

NODE 

N0D3 

NOD* 

HDD 

HT 

PT 

RHO 

RHO 

SN 

GEAR 

SND 

TRQD 

HD 

W 

HT 

PT 

TORQ 


GLOIAL 

GLOIAL 

MAP 

NAME 0 

NAME 1 

NAME 2 

NAME 3 

NAME * 

DESIGN 

VARIABLE 

VARIABLE 

VARIABLE 

VARIABLE 

VARIABLE 

DESIGN 

DESIGN 

DESIGN 

DESIGN 

VARIABLE 

VARIABLE 

VARIABLE 

VARI All E 


IN 0 
IN 0 
IN 0 
IN 0 
IN 0 
IN 0 
IN 0 
IN 0 
IN 0 
IN 2 
IN 2 
IN 2 
IN 5 
IN 4 
IN 0 
IN 0 
IN 0 
IN 0 
IN 1 
OUT 3 
OUT 3 
OUT 0 



Tha configurat ion proct*r,or use* the following intorfaco cards 
to provide units for parameter's whose output in requested by e 
PRPLOl call from the main program 
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IPRPL 

D'LESS 

D'LESS 

1UPDAT 

D'LESS 

D'LESS 

hap 

D'LESS 

D'LESS 

HODN 

D'LESS 

D'LESS 

MODI 

D'LESS 

D'LESS 

•M0D2 

D'LESS 

D'LESS 

NODS 

D'LESS 

D'LESS 

MO DA 

D'LESS 

D'LESS 

HDD 

IN 

H 

HTIM 

ITU/LBH 

J/KG 

PTIM 

LIF/INMI2 

H/M**2 

RHOIN 

LlM/IM*i2 

K0/M**3 

RHOOUT 

LSK/IN1W2 

K0/H*»3 

SNREF 

RPM 

RPM 

ORATIO 

D'LESS 

D'LESS 

SND 

RPM 

RPM 

TROD 

IN-LDF 

M-M 

HO 

L1H/S 

KG/S 

HIM 

LIM/S 

KG/S 

HTDUT 

ITU/ L DM 

J/KG 

PTOUT 

LIF/IN«*2 

M/M**2 

TORO 

IN-LBF 

N-M 


Following is o list of the subroutines thot oro required by this Modulo. 

SUBROUTINES REQUIRED : PMAPXX HAP (EXTERNAL) 

PRriOl 

Following is s list of the commons thot or# ronui rod by this modulo. 
COMMONS REQUIRED < C-UNIT3 


3.4.24 

PUMP01 

4 
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Appendix C 

Interfaced NASA Control Model 

Presented is the listing of the NASA MSFC FORTRAN Control Model with the ROCETS interface 
incorporated in the comment cards. 
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SUBROUTINE CNTLOO 


SUBROUTINE CNTLOOdPRPL 

, IUPDAT 

, TINE 

9 

NOON 

, 9N 

i 

t 

PCN 

, PFDN 

, REFN 

9 

QFFN 

, PC , 

z 


PFD1 

, TFP1 

, PCREF 

9 

NRREF 

, DXOPI , 

3 

$ 

DXFPI 

, XMFVC 

, XHOVC 

w 

XCCVC 

, XFPVC , 

4 

$ 

XOPVC 

, EOPI 

, EFPI 

9 

AAA 

) 

5 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

6 

C XBEGIN CLASS CHTLOO 






X 

7 

C 






X 

8 

C SUBPROGRAM CNTLOO 


UNCLASSIFIED 


SID: 

E950 x 

9 


c 


X 

10 

C JCEND CLASS CNTLOO 

X 

11 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX)(XXX 

1Z 

C XBEGIN 

PURPOSE CNTLOO 

X 

13 

C 


X 

14 

c 

SSNE CONTROL FOR MAINSTAGE OPERATION 

X 

15 

C 


X 

16 

C SEND PURPOSE CNTLOO 

X 

17 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

18 

C XBEGIN 

HISTORY CNTLOO 

X 

19 

c 


X. 

ZO 

c 

OBTAINED FROM NASA/MSFC FOR TESTING ROCETS SYSTEM MAY 

1990 x 

Z1 

c 


X 

ZZ 

C XEND HISTORY CNTLOO 

X 

23 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

Z4 

C XBEGIN 

DESCRIPTION CNTLOO 


25 

c 



26 

C IPRPL 

PRPLOl OUTPUT FLAG 

D'LESS 

27 

c 

0 = NO PRINT 


28 

c 

1 * PRINT 


29 

C IUPDAT 

UPDATE FLAG 

D'LESS 

30 

c 

-1 * INITIAL IZATION/SS BALANCE 


31 

c 

0 » CONTROL BYPASSED DURING ITERATION 


32 

c 

1 - CALCULATIONS PERFORMED AFTER CONVERGED POINT 


33 

C TIME 

SIMULATED TIME 

SEC 

34 

C MODN 

MODULE NAME 

D'LESS 

35 

C QN . 

FUEL FLOW FEEDBACK NODE 

D'LESS 

36 

C PCN 

CHAMBER PRESSURE NODE 

D'LESS 

37 

C PFDN 

FUEL PUMD DISCHARGE NODE 

D'LESS 

38 

C REFN 

NODE NAME FOR 'REFERENCE' INPUTS CPCREF, MRREF) 

D'LESS 

39 

C QFFM 

MEASURED FUEL FLOW 

OPM 

40 

C PC 

MEASURED CHAMBER PRESSURE 

PSIA 

41 

C PFD1 

MEASURED LOW PRESSURE FUEL PUMP DISCHARGE PRESSURE 

PSIA 

42 

C TFP1 

MEASURED LOW PRESSURE FUEL PUMP DISCHARGE TEMPERATURE 

DEG. R 

43 

C PCREF 

REFERENCE (REQUESTED) CHAMBER PRESSURE 

PSIA 

44 

C NRREF 

REFERENCE (REQUESTED) MIXTURE RATIO 

D-LESS 

45 

C 



46 

C INPUT ON IUPDAT=-1, OTHERWISE OUTPUT 


47 

C 



48 

C DXOPI 

OPV INTEGRATOR VALUE 

PCT 

49 

C DXFPI 

FPV INTEGRATOR VALUE 

PCT 

50 
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SUBROUTINE CNTLOO 


OUTPUTS 

XNFVC 
XMOVC 

xccvc 

XFPVC 
XOPVC 
EOPI 
EFPI 
AAA 


HFV ACTUATOR POSITION COMMAMD 
HOV ACTUATOR POSITION COMMAND 
CCV ACTUATOR POSITION COMMAND 
FPV ACTUATOR POSITION COMMAND 
OPV ACTUATOR POSITION COMMAND 
OPV INTEGRATOR ERROR 
FPV INTEGRATOR ERROR 
STORAGE ARRAY FOR RESTART 


PCT 

pct 

PCT 

PCT 

PCT 


_ XEND DESCRIPTION CNTLOO 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

XBEGIN COMMENTS CNTLOO * 

X 


C XEND COMMENTS CNTLOO 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C XBEGIN INTERFACE CNTLOO K 

C “ 


Cl 

Cl 

CALL LISTI SYSTEM 
NAME 1 NAME 

| SYSTEM TAG 
1 

| ARRAY 

| STATUS 

| I/O 
| STATUS 

VAR 

TYPE 

1* 

|x 

— X 

c — 

C| 

IPRPL 

I IPRPL 

| GLOBAL 

1 

| IN 0 

1X4 

1* 

C| 

IUPDAT 

| IUPDAT 

| GLOBAL 

1 

1 IN 8 

1*4 

|X 

cl 

TIME 

| TIME 

| GLOBAL 

1 

I IN • 

R*4 

|X 

Cl 

MODN 

| MODN 

I NAME 0 

1 

I IN 0 

C*4 

1* 

Cl 

ON 

1 ON 

| NAME 1 


I IN • 

C*4 

1* 

C| 

PCN 

I PCN 

1 NAME 2 

1 

I IN 0 

CX4 

1* 

C| 

PFDN 

| PFDN 

| NAME 3 

1 

I IN B 

CX4 . 

1* 

C| 

REFN 

| REFN 

| NAME 4 

1 

I IN • 

C*4 

1* 

Cl 

QFFM 

1 0 

| VARIABLE 

1 

| IN 1 

R*4 

1* 

Cl 

PC 

1 PT 

| VARIABLE 

1 

I IN 2 

RX4 

1* 

Cl 

PFD1 

1 PT 

| VARIABLE 

1 

I IN 3 

RX4 

1* 

Cl 

TFP1 

1 TT 

| VARIABLE 

1 

I IN 3 

RX4 

1* 

Cl 

PCREF 

I PC 

| VARIABLE 

1 

I IN 4 

Rx4 

1* 

Cl 

MRREF 

1 MR 

| VARIABLE 

1 

| IN 4 

R*4 

1* 

Cl 

DXOPI 

| DXOP 

I VARIABLE 

1 

I IN • 

Rx4 

1* 

Cl 

DXFPI 

I DXFP 

I VARIABLE 

1 

1 IN B 

RX4 

|X 

Cl 

XMFVC 

| XMFV 

| VARIABLE 

1 

| OUT 0 

RX4 

|X 

Cl 

XMOVC 

I XMOV 

| VARIABLE 

1 

| OUT • 

Rx4 

|X 

cl 

XCCVC 

I XCCV 

| VARIABLE 

1 

I OUT i 

R*4 

1* 

cl 

XFPVC 

| XFPV 

| VARIABLE 

1 

I OUT 0 

R*4 

1* 

Cl 

XOPVC 

| XOPV 

| VARIABLE 

1 

| OUT 0 

R*4 

1* 

Cl 

EOPI 

| EOPI 

| VARIABLE 

1 

I OUT 0 

RX4 

1* 

Cl 

EFPI 

1 EFPI 

| VARIABLE 

1 

| OUT • 

RX4 

1* 

Cl 

AAA 

I AAA 

| VARIABLE 

| ARRAY 21 

| OUT • 

RX4 

IX 

X 



51 

52 

53 

54 

55 
54 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
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SUBROUTINE CNTLOO 


C SUBROUTINES REQUIRED i PRPL01 * 

C * 

C XEND SUBROUTINES REQUIRED CNTLOO * 

C XBEGIN COWNMS REQUIRED CNTLOO * 

c ; 

C COMMONS REQUIRED « NONE * 

c * 

c XEND COMMONS REQUIRED CNTLOO * 

Cxx******************************************************************** 
CHARACTER * 4 MODN, QN, PCN, PFDN, REFN 


REAL MRREF, NRCONT , MRG , MRREFX , MRG6 5 , MRGI 00 

REAL HFVRL ,NOVRL , MFVRX, MOVRX 

DIMENSION AITAB(10),A1TAB(10) , EPLTAB(IO) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


c 

c CONTROL ER GAINS 

c 

C PCG50 LON CORNER POINT FOR CHAMBER PRESSURE PROPORTIONAL GAIN 

C P CGI 00 HIGH CORKER POINT FOR CHAMBER PRESSURE PROPORTIONAL GAIN 

C PCRL CHAMBER PRESSURE RATE LIMIT 

C 

DATA PCG50 /.6/» PCG100/1.0/, PCRL/300./ 


C 

C XFG50 LON CORNER POINT FOR MAINSTAGE CROSSFEED GAIN 

C XFG100 HI6M CORNER POINT FOR MAINSTAGE CROSSFEED GAIN 

C 

DATA XF650 /l. 15 /» XFG100/1.15/ 


C 

C MRG65 LOM CORNER POINT FOR MIXTURE RATIO PROPORTIONAL GAIN 

C MRGI 00 HIGH CORNER POONT FOR MIXTURE RATIO PROPORTIONAL GAIN 

C 

DATA MRG65 /.2/» MRG100/.5/ 


C 

C XOPPG OfOV PROPORTIONAL GAIN 

C XOPIG OPOV INTEGRAL GAIN 

C XOPVST START BIAS FOR OPOV 

C XOPDCO OfOV DELTA COMMAND OFFSET 

C 

DATA XOPPG /.BUS/, XOPIG/ .00068/, XOPVST/64.52/# XOPDCO/ 0./ 


C 

C XFPPG FPOV PROPORTIONAL GAIN 

C XFPIG FPOV INTEGRAL GAIN 

C XFPVST START HAS FOR FPOV 

C 

DATA XFPPG / 7./, XFPIG/. 40/, XFPVST/77 .22/ 


C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c RATE LIMITS 


151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
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SUBROUTINE CNTLOO 


C SUBROUTINES REQUIRED « PRPlOl * 

C K 

C XEND SUBROUTINES REQUIRED CNTLOO 

cxxxxxxxxxxxxxxjutxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C XBE6IN COfMONS REQUIRED CNTLOO * 

C x 

C COMMONS REQUIRED > NONE * 

C x 

C XEND COMMONS REQUIRED CNTLOO * 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

CHARACTER X A MODN, QN, PCN, PFDN, REFN 


REAL MRREF, NRCONT , MRC , MRREFX , MRG6 5 , MRG1 0 0 

REAL MFVRL»NOVRL »MFVRX,MOVRX 

DIMENSION AOTAB(10),A1TABC10),EPLTABC10) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CONTROLER GAINS 

PCG50 LON CORNER POINT FOR CHAMBER PRESSURE PROPORTIONAL GAIN 

PCGIOO HIGH CORMER POINT FOR CHAMBER PRESSURE PROPORTIONAL GAIN 

PCRL CHAMBER PRESSURE RATE LIMIT 

DATA PCG50 /.6/» PCG100/1.0/, PCRL/300./ 

XFG50 LON CORNER POINT FOR MAINSTAGE CROSSFEED GAIN 

XFGIOO HIGH CORNER POINT FOR MAINSTAGE CROSSFEED GAIN 

DATA XFG50 /l. 15 /. XFG100/1.15/ 


C 

C MRG65 LON CORNER POINT FOR MIXTURE RATIO PROPORTIONAL GAIN 

C MRG100 HI6H CORNER POONT FOR MIXTURE RATIO PROPORTIONAL GAIN 

C 

DATA MRG65 /. 2/, MRG100/.5/ 


C 

C XOPPG OPOV PROPORTIONAL GAIN 

C XOPIG OPOV INTEGRAL GAIN 

C XOPVST START BIAS FOR OPOV 
C XOPDCO OPOV DELTA COMMAND OFFSET 
C 

DATA XOPPG /. 0113/, XOPIG/ .00068/, XOPVST/64.52/, XOPDCO/ 0./ 

C 

C XFPPG FPQV PROPORTIONAL GAIN 

C XFPIG FPOV INTEGRAL GAIN 

C XFPVST START HAS FOR FPOV 
C 

DATA XFPPG / 7./, XFPIG/. 40/, XFPVST/77 .22/ 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c RATE LIMITS 


151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
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non 


SUBROUTINE CNT100 


u 201 

DATA CCVRL /200./,FPVRL/200./,HFVRL/200./,OPYRL/200./,MOVRL/200./ 202 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 203 
C CONSTANTS FOR FUEL DENSITY E9UATION 20A 

C 205 

C RHO = (AO + (A1+B1*PLPFD)XTLPFD ♦ (A2+B2XPLPFD)XTLPFDXTLPFD 206 

C 207 

DATA A0/ .38956E+01/, Al/ .6522E-01 /, A2/-.14013C-02/, 208 

* BO/ .42739E-02/, B1/-.21467E-03/, B2/ .30926E-05/ 209 

Cxxxxxxttxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxifxxxxxxxxinonnnaixxxxxxxxxxxxxxxxx 210 

C CONVERSION FACTOR. (448 CPU * 1 FTM3/SEC) 211 

DATA GTOC/448. 83303/ 212 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 213 

C 214 

c CONSTANTS FOR LOX FLOW CALCULATION 215 

C 216 

C C2 = C2A X (PC/RPL)XX2 + CZBx(PC/RPL) ♦ C2C 217 

C WO = (PC + 14.5) / C2 - WH 218 

J 219 

C NOTE. USE PCREF INSTEAD OF PC DURING THRUST LIMITING 220 

C 221 

DATA C2A/-. 030621/, C2B/. 016555 /, C2C/2. 92104/ 222 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXKXXXXXXXXXXXXXXXXXXXXXXX 223 

DATA AOTAB/ 42.75, 42.75, 39.00, 35.00, 35.00, 17.00, 9.40, 224 

0 9.40, 9.40,-32.40/ 225 

DATA A1TAB/ 0.20, 0.20, 0.25, 0.30, 0.30, 0.50, 0.58, 0.58, 226 

^ 0.58, 0.96/ 227 

DATA EPLTAB/ 0., 70., 75., 80., 85., 90., 95., 100., 105., 228 

0 110.0/ 229 

230 

231 

232 

DATA RPL / 3006. / , XOPVMX/ 64.43 / , IFIIST/ 1 / 233 

DIMENSION AAAC21) 234 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXX 235 

c INITIALIZATION 236 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXIXXXXXXXXXXXXXXXXXXXX 237 

IFdUPDAT .LT. 0)THEN 238 

C 239 

TNEXT = TIME-. 0025 240 

MNCYC = 3 241 

9FX * QFFM 242 

PCX = PC 243 

PFD1X = PFD1 244 

TFP1XX * TFP1 245 

TFP1X « TFP1XX 246 

PCRFXL ■ PCREF 247 

PCREFX « PCREF 248 

MRREFX * MRREF 249 

EOPIL * 0. 250 
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SUBROUTINE CNTLOO 


EFPIL =0. 251 

252 

DXFPIL = DXFPI 253 

DXOPIL = DXOPI 254 

255 

RPL50 = RPL * .5 256 

RPL65 * RPL x .65 257 

258 

CONVERT PC RATE LIMIT TO PSI PER 28 MSEC 259 

PCRLX * PCRL / 50. 260 

SCHEDULE DELTAS 261 

DPCG * PCG100 - PCG50 262 

DMRG * MRG100 - MRG65 263 

DXFG * XFG100 - XFG50 264 

DXCCV -100.-52. 265 

266 

COMPUTE OPOV DELTA POWER LEVEL 267 

OPOVDL - (1.5 * XOPVMX) - 97.5 268 

269 

CONVERT VALVE RATE LIMITS FROM X PER SECOND TO X PER 20 MSEC 270 

CCVRX = CCVRL / 50. 271 

FPVRX * FPVRL / 50. 272 

MFVRX = MFVRL / 50. 273 

MOVRX » MOVRL / 50. 274 

OPVRX = OPVRL / 50. 275 

IFIRST= 0 276 

ELSEIFC IFIRST . EQ . 1) THEN 277 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 278 

C UNLOAD THE ARRAY CONTAINING THE VARIABLES REQUIRED FOR RESTART X 279 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 280 

TNEXT » AAA( 1) 281 

MNCYC » IFIX(AAA( 2)+.l) 282 

QFX - AAA( 3) 283 

PCX * AAA( 4) 284 

PFD1X * AAA( 5) 285 

TFP1XX * AAA< 6) 286 

TFP1X » AAAC 7) 287 

PCRFXL * AAA( 8) 288 

PCREFX ■ AAAC 9) 289 

MRREFX = AAA (10) 290 

EOPIL » AAA(ll) 291 

EFPIL > AAA ( 12) 292 

DXFPIL = AAAC 13) 293 

DXOPIL = AAA(14) 294 

C 295 

XOVPS ■ AAAC 15) 296 

XFPVC » AAAC 16) 297 

XCCVXL - AAAC 17) 298 

XFPVXL » AAAC IS) 299 

. XMFVXL » AAA(19) 300 
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SUBROUTINE CNTLOO 


XMOVXt = AAA (20) 301 

XOPVXL = AAAC21) 302 

C 303 

RPL50 » RPL * .5 304 

RPL65 » RPl x .65 305 

PCRLX ■ PCRL / 50. 306 

DPCO » PCG100 - PCG50 307 

9MRG * NROIOO - MRG65 30* 

DXFG * XFG100 - XFG50 309 

DXCCV * 100. - 52. 310 

OFOVDL * (1.5 * XOPVMX) - 97.5 311 

CCVRX * CCVRL / 50. 312 

FPVRX = FPVRL / 50. 313 

MFVRX = NFVRL / 50. 314 

MOVRX * NOVRL / 50. 315 

OPVRX = OPVRL / 50. 316 

IFIRST » 0 317 

ENDIF 310 

IFdUPDAT .EO. 0)GO TO 50 319 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 320 

C TRANSIENT CONTROL SECTION 321 

Cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx XXXXXXXXXXXX 322 
IF(TIHE .LT. TNEXT) GO TO 50 323 

TNEXT « TINE ♦ 0.005 324 

MNCYC * HHCYC ♦ 1 325 

GO TO (100,200* 300,400)* MNCYC 326 

C 327 

C MINOR LOOP 1 328 

100 9FX » QFFM 329 

PCX = PC 330 

PFD1X * PFDI 331 

GO TO 50 332 

C 333 

C MINOR LOOP 2 334 

200 CONTINUE 335 

GO TO 50 336 

C 337 

C MINOR LOOP 3 338 

300 PCREFX * PCtEF 339 

MRREFX = MRREF 340 

TFP1X * TFP1XX 341 

TFP1XX * TFP1 342 

GO TO 50 343 

C 344 

C MINOR LOOP 4 345 

400 PCPCTX » PCX / RPL 346 

PCPCT * PCPCTX X 100. 347 

MNCYC * 0 348 

C 349 

C RATE LIMIT ON PCREF 350 
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SUBROUTINE CNTLOO 


410 D - PCREFX - PCRFXL 

IF(ABS(D) .GT.PCRLX) PCREFX = PCRFXL + SIGN(PCRLX,D) 

4X2 PCRFXL - PCREFX 
C 

C COMPUTE OPOV COMMAND LIMIT 

420 EPL ■ (PCREFX * 100.0 / RPL) + OPOVDL 
DO 421 1=10, 1,-1 

IFCEPL .GE. EPLTAB(I)) GO TO 422 

421 CONTINUE 
1*1 

422 OPOVCL = (AITAB(I) * EPL) + AOTAB(I) ♦ XOPDCO 
OPOVCL = MIN( OPOVCL. 100.0) 

C 

C COMPUTE GAIN SCHEDULES 
430 DRPL50 = (PCX - RPL 50) / (RPL - RPL50) 

PCG ■ PCG50 + (DRPL50 X DPCG) 

PCG = MAX( PCG50 »MIN(PCG1 00 , PCG) ) 

XFG * XFG50 + (DRPL50 X DXFG) 

XFG * MAX(XFG50,MI«XFG100,XFG)) 

DRPL65 * (PCX - RPL65) / (RPL - RPL65) 

NRG * MRG65 ♦ (DRPL65 X DMRG) 

MRG = MAX(MRG65»MIH(MRG100,MRG)) 

C 

C PC ERROR A PROPORTIONAL 

450 DPC ■ PCREFX - PCX 
EOPV * DPC x PCG 
DXOPP = EOPV X XOPPG 

C 

C OPV INTEGRATOR 
EOPI = EOPV 
C 

C -CHECK FOR THRUST LIMIT 

IFdUPDAT .LT. 0)X0VPS = XOPVST ♦ DXOPI 
IFdUPDAT .GT. 0) THEN 

IF((XOPVS.GE. OPOVCL). AND. (EOPI. GT. 0.0)) EOPI * 0.0 
ENDIF 
C 

451 IFdUPDAT .GT. 0) DXOPI = DXOPIL + (XOPIG*(EOPI ♦ EOPIL)) 

C 

DXOPIL = DXOPI 
EOPIL ■ EOPI 

C OPOV SUM A LIMIT CHECK 

DXOPV ■ DXOPP ♦ DXOPI 

XOPVS = DXOPV + XOPVST 

XOPVX * XOPVS 

C 

C COMPUTE FUEL DENSITY (RHOH) AND FLOWRATE (WH) 

470 RHOH = ((B2XPFD1X ♦ A2) * TFP1X ♦ (B1*PFD1X + AD) * TFP1X ♦ 
0B0XPFD1X + AO 
C 


351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 
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SUBROUTINE CNTLOO 


C CONVERT Q FROM GAL/MIN TO CU FT/SEC 
QFC ■ QFX / GTOC 
HH = QFC * RHOH 
C 

C CALCULATE OXIDIZER FLOWRATE (HO) AND MIXTURE RATIO 
C USE MEASURED PC IF IN NORMAL MODE 
C USE PCREF IF IN THRUST LIMITING MODE 

C 

T = PCX 
TT = T / RPL 

C2 = CCCZAXTT ♦ C2B) * TT) ♦ C2C 
HO = (CT + 14.51/C2) - HH 
MRCONT = HO / HH 
C 

C CROSSFEED 

DXFPX * DXOPV * XFG 
C 

C FPV CONTROL 
500 DMR ■ MRCONT - MRREFX 
EFPV ■ DMR x MRG 
EFPI =* EFPV 
C XFPV ERROR LIMIT 

IFCIUPDAT AT. 0)XFPVC * DXFPX + DXFPI + XFPVST 
C 

IFdUPDAT ,GT. 0)THEN 

IF((XFPVC.GE. 102.0). AMD. (EFPV. GT.O.)) EFPI ■ 0. 

IF( (XFPVC.LT. 0.00) .AND. (EFPV. LT.O.)) EFPI = 0. 

ENDIF 

IFdUPDAT .GT. 0) DXFPI * DXFPIL + (XFPIG X (EFPI ♦ EFPIL)) 
C 

DXFPIL * DXFPI 
EFPIL * EFPI 
DXFPP * XFPPG * EFPV 
DXFPV = DXFPI ♦ DXFPP 
C XFPV SUM 

XFPVX a DXFPV ♦ DXFPX + XFPVST 
XFPVX » MIN(XFPVX, 100.0) 

C 

C SCHEDULED VALVES 
80 XMFVX « 100. 

XMOVX ■ 100. 

XCCVX * 52. ♦ (DRPL50 X DXCCV) 

XCCVX * MAX(52.,MIN( 100., XCCVX)) 

C 

C VALVE RATE LIMITS 
C 


401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 
428 

421 

422 

423 

424 

425 

426 

427 

428 

429 
438 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 


IFdUPDAT .LT. 0)THEN 
XCCVXL = XCCVX 
XFPVXL = XFPVX 
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XMFVXL * XMFVX 
XMOVXL « XMOVX 
XOPVXL * XOPVX 
ENDIF 
C 

90 XCCVC = XCCVX 

D * XCCVX - XCCVXL 

IFCABSCD) .GT.CCVRX) XCCVC = XCCVX l + SIGHCCCVRX,D) 
92 XFPVC ■ XFPVX 

D = XFPVX - XFPVXL 

IFCABSCD) .6T.FPVRX) XFPVC « XFPVXL ♦ SIGNCFPVRX,D) 
96 XMFVC « XMFVX 

D = XMFVX - XMFVXL 

IFCABSCD) .GT.MFVRX) XMFVC = XMFVXL + SIGNCMFVRX,D) 
96 XMOVC » XMOVX 

D = XMOVX - XMOVXL 

IFCABSCD) .GT.MOVRX) XMOVC * XMOVXL + SIGNCMOVRX, D) 
98 XOPVC « XOPVX 

D = XOPVX - XOPVXL 

IFCABSCD) .GT.OPVRX) XOPVC » XOPVXL + SIGNCOPVRX,D) 
C 

XCCVXL = XCCVX 
XFPVXL » XFPVX 
XMFVXL * XMFVX 
XMOVXL * XMOVX 
XOPVXL « XOPVX 

Cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C LOAD STORE ARRAYS FOR RESTART CAPABILITY * 
cxxxxxxxxxxxxjraxx*x*xx**x**x*xxxxxxx***xxxxxxx 


AAA( 1) 

3 

TNEXT 

AAAC 2) 

3 

FLOATCMNCYC) 

AAA( 3) 

s 

QFX 

AAAC 6) 

x 

PCX 

AAAC 5) 

x 

PFD1X 

AAAC 6) 

X 

TFP1XX 

AAAC 7) 

X 

TFP1X 

AAAC 8) 

s 

PCRFXL 

AAAC 9) 

s 

PCREFX 

AAA(IO) 

3 

MRREFX 

AAAC11) 

3 

EOPIL 

AAAC12) 

S 

EFPIL 

AAAC13) 

X 

DXFPIL 

AAAC16) 

= 

DXOPIL 

AAAC 151 

3 

XOVPS 

AAAC 16) 

3 

XFPVC 

AAAC 17) 

3 

XCCVXL 

AAAC18) 

3 

XFPVXL 

AAAC19) 

3 

XMFVXL 

AAAC 20) 

3 

XMOVXL 


651 

652 

655 

656 

655 

656 

657 

658 

659 

660 
661 
662 
663 
666 

665 

666 

667 

668 

669 

670 

671 

672 

673 
676 

675 

676 

677 

678 

679 

680 
681 
682 
683 
686 

685 

686 

687 

688 

689 

690 

691 

692 

693 
696 

695 

696 

697 

698 

699 
500 
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SUBROUTINE CNTLOO 


AAAC21) = XOPVXL 
SO CONTINUE 
Cxxxxxxxxxxxxxx 
C PRINT CALLS X 

CXXXXXXXXXXXXXX 


IFCIPRPL . ER. 0)GQ TO 99 


CALL 

PRPL01C-9 

, • MODULE • , MODM// * CHTLOO ', DUMMY ) 

CALL 

PRPLOKl, 

•AO •,• 

* , AO 

) 

CALL 

PRPLOKl, 

•A1 •,• 

' » A1 

> 

CALL 

PRPLOKl, 

•A2 S' 

• ,A2 

) 

CALL 

PRPLOKl, 

•BO 

* ,10 

) 

CALL 

PRPLOKl, 

•Bl •,' 

',11 

) 

CALL 

PRPLOKl, 

•B2 S* 

',12 

) 

CALL 

PRPLOKl, 

•CCVRL •,• 

•.CCVRL 

1 

CALL 

PRPLOKl, 

•CCVRX •,' 

•.CCVRX 

) 

CALL 

PRPLOKl, 

■C2 

' ,C2 

) 

CALL 

PRPLOld, 

■C2A •,• 

*,C2A 

> 

CALL 

PRPL0K1, 

•C2B •,’ 

' ,C2B 

) 

CALL 

PRPLOKl, 

•C2C ■,* 

*,C2C 

1 

CALL 

PRPLOKl, 

iD •,• 

* ,D 

) 

CALL 

PRPLOKl, 

■DHR •,• 

• ,DMR 

) 

CALL 

PRPLOKl, 

•DMRG 

*,DMRG 

) 

CALL 

PRPLOKl, 

•DPC •»• 

*»DPC 

) 

CALL 

PRPL01C1, 

•DPCG 

1 , DPCG 

) 

CALL 

PRPL01 (1, 

■DRP150 V f 

* ,DRPL50 

} 

CALL 

PRPLOKl. 

•DRPL65 •,* 

'.DRPL65 

) 

CALL 

PRPLOKl, 

•DXCCV •,• 

', DXCCV 

) 

CALL 

PRPLOKl, 

•DXFG ■;» 

• , DXF6 

) 

CALL 

PRPLOKl, 

•DXFPI •,’ 

'.DXFPI 

) 

CALL 

PRPLOKl, 

•DXFPIl •,* 

'.DXFPIL 

> 

CALL 

PRPLOKl, 

•DXFPP •#• 

', DXFPP 

) 

CALL 

PRPLOKl, 

•DXFPV •,• 

', DXFPV 

) 

CALL 

PRPLOKl, 

• DXFPX •,• 

', DXFPX 

) 

CALL 

PRPLOKl, 

•DXOPI •„• 

•.DXOPI 

) 

CALL 

PRPLOKl, 

•DXOPIL •,• 

', DXOPIL 

> 

CALL 

PRPLOKl, 

• DXOPP •»• 

», DXOPP 

) 

CALL 

PRPLOKl, 

•DXOPV •,* 

'. DXOPV 

) 

CALL 

PRPLOKl, 

•EFPI •,• 

' , EFPI 

) 

CALL 

PRPLOKl, 

•EFPIL •,• 

'.EFPIL 

) 

CALL 

PRPLOKl, 

•EFPV •,* 

'.EFPV 

) 

CALL 

PRPLOKl, 

•EOPI •,• 

*,EOPI 

) 

CALL 

PRPLOKl, 

•EOPIL •,* 

', EOPIL 

) 

CALL 

PRPLOKl, 

•EOPV •,• 

'.EOPV 

) 

CALL 

PRPLOKl, 

•EPL 

'.EPL 

) 

CALL 

PRPLOKl, 

•FPVRL •,• 

• , FPVRL 

) 

CALL 

PRPLOKl, 

•FPVRX •,* 

', FPVRX 

) 

CALL 

PRPLOKl, 

•GTOC •,• 

'.GTOC 

) 

CALL 

PRPLOKl, 

•I ", f 

• , FLOATd 


CALL 

PRPLOKl, 

•IPRPi •,* 

'.FLOATdPRPL 

CALL 

PRPLOKl, 

•IUPDAT •.• 

', FLOAT (IUPDAT 


501 

502 

503 

504 

505 

506 

507 
SOS 

509 

510 

511 

512 

513 

514 

515 

516 

517 
51S 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 
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SUBROUTINE CNTLOO 
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SUBROUTINE CNTLOQ 


CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPL01C1 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPL0K1 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
CALL PRPLBKl 
RETURN 
END 


•XFG5B 

»,XFG50 

) 

601 

•XFPIG 

*, XFPIG 

) 

602 

•XFPPG •,' 

•.XFPPG 

) 

603 

•XFPVC ' . • 

*, XFPVC. 

) 

604 

'XFPVST *,* 

•.XFPVST 

) 

605 

* XFPVX 

•.XFPVX 

) 

606 

•XFPVXL 

•.XFPVXL 

) 

607 

•XMFVC * , ' 

•.XMFVC 

) 

60S 

•XMFVX »,» 

».XMFVX 

) 

609 

•XMFVXL • , » 

•.XMFVXL 

) 

610 

•XMOVC 

•.XMOVC 

) 

611 

•XMOVX • . * 

•.XMOVX 

) 

612 

•XMOVXL • . • 

•.XMOVXL 

) 

613 

•XOPDCO » , • 

'.XOPDCO 

} 

614 

•XOPIG 

•.XOPIG 

) 

615 

•XOPPG »,» 

•.XOPPG 

) 

616 

•XOPVC ',• 

•.XOPVC 

) 

617 

•XOPVMX 

'.XOPVMX 

) 

61S 

•XOPVS *,' 

•.XOPVS 

) 

619 

•XOPVST •,' 

XOPVST 

) 

620 

•XOPVX ' , ' 

'.XOPVX 

) 

621 

•XOPVXL •,' 

'.XOPVXL 

) 

622 

•XOVPS ',• 

'.XOVPS 

) 

623 


624 

625 


Appendix D 

TTBE Model Configration Input 

The detail TTBE model described in this report was delivered to NASA-MSFC. Presented in 
this appendix is the listing of the configuration input which the ROCETS system interprets to 
generate the TTBE simulation without the control model. 


177 



TTBB I/O CONTROL 


PACT 1 


CONFIG. INPOT 


************************************************************************ 1 

• ROCETS CONFIGURATION INPUT FILE * 2 

*- * 3 

• * 4 

* MODEL* ITBE * 5 

* * 6 

* CONFIGURATION FILE i CITBB002 * 7 

* * 8 

* PURPOSE i USED TO CENERATE TTBB WITHOUT CONTROL * 9 

* * 10 

* DATE: 08/15/90 * 11 

* * 12 

* ENGINEER: T.F. DENMAN, PRATT 1 WHITNEY * 13 

* *14 

************************************************************************ 15 

********************************************** 16 

* DEFINE CONFICRUATION PROCESSOR OPTIONS * 17 

********************************************** 18 

DEFINE OPTIONS 19 

UNITS i ENGLISH} 20 

TITLE i TTBB WITHOUT THE CONTROL - GENERATED FROM CTTBE002; 21 

CROSS i ON | 22 

PDS ! HDHS206. ROCETS. DATADICTj 23 

MAXIMUM ! SYMXSCHD - 20 24 

SYMXSCVL - 5000 t 25 

END OPTIONS 26 

************************************************************* 27 

* THE MAXIMUMS ON THE FOLLOWING VARIABLES CAN BE ALTERED * 28 

* BY USING THE MAXIMUM KEYWORD IN THE OPTIONS BLOCK * 29 

************************************************************* 30 

* MAXIMUM i SYMXNAMB - 10000 , 31 

* SYMXSCHD 20 32 

* SYMXSCVL - 5000 , 33 

* SYMXLIN - 100 , 34 

* SYMXLOUT 50 35 

* SYMXSTAT 200 36 

* SYMXBAL - 50 , 37 

* SYMXITMP - 20000 , 38 

* SYMXRTMP - 20000 , 39 

* SYMXCTMP - 20000 , 40 

* SYMXDTMP - 20000 , 41 

* SYMXPRNT - 10000 , 42 

* SYMXBLOK - 5000 , 43 

* SYMXCOL 10 44 

* SYMXDISC 500 j 45 
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TTBH W/0 CONTROL 


PAGE 2 


CONFIG. INPUT 


* DEFINE MODULES FROM OUTSIDE LIBRARY • 47 

********************************************** 48 

* 49 

* DEFINE INSTREAM SO 

* 51 

* FOR EXAMPLE i B2 

* S3 

* PIPE10 : HDHS206. PIPE. FORTRAN ; S4 

* 50 

* END INSTREAM 06 

* 57 

********************************************** 88 

* DEFINE ADDITIONAL EXTERNAL INPUTS • 59 

********************************************** 60 

DEFINE EXTERNALS 61 

PRIMMDI , 62 

PRIMFPB , 63 

PRIMOPB , 64 

PCREQ t 60 

END EXTERNALS 66 

********************************************** 67 

* CHANGE ITERATION VARIABLES FOR STATES * 68 

********************************************** 69 

DEFINE INTEGRATION 70 

ITERATE: HIPBSF FOR UIPBSF j 71 

ITERATE: BXPBSO FOR UTPBSO » 72 

ITERATE: HIYL1 FOR UTVL1 j 73 

ITERATE: ETVL10 FOR UTVL10 j 74 

ITERATE: HTVL11 FOR UTVL11 » 75 

ITERATE: ETVL12 FOR UTVL12 j 76 

ITERATE: BTVL13 FOR UTVL13 » 77 

ITERATE: HTVL14 FOR UTVL14 > 78 

ITERATE: KTVL1S FOR UTVL15 j 79 

ITERATE: BTVL16 FOR UTVL16 j 80 

ITERATE: ETVL17 FOR UTVL17 j 81 

IIERATE: HTVL18 FOR UTVL18 j 82 

ITERATE: HTVL19 FOR UTVL19 j 83 

ITERATE: BTVL2 FOR UTVL2 j 84 

ITERATE: HTVL20 FOR UTVL20 ; 85 

ITERATE: HTVL21 FOR UTVL21 j 86 

ITERATE: HTVL22 FOR UTVL22 j 87 

ITERATE: HTVL3 FOR UTVL3 j 88 

ITERATE: HTYL4 FOR UTVL4 ; 89 

ITERATE: HTVL5 FOR UTVL5 j 90 
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HBB W/0 CONTROL 


PACE 3 


CONFIG. INPUT 


ITERATE: 

HTVL6 

FOR UTVL6 i 


91 

ITERATE: 

BTVL7 

FOR UTVL7 | 


92 

ITERATE i 

HTVL8 

FOR UTVL8 | 


93 

ITERATE: 

HTVL9 

FOR UTVL0 | 


94 

ITERATE: 

PTPBSF 

FOR RHOPBSF i 


98 

ITERATE: 

PTPBSO 

FOR RHOPBSO j 


96 

ITERATE: 

PTVLl 

FOR RH0VL1 | 


97 

ITERATE: 

PTVL10 

FOR RH0VL10 i 


98 

ITERATE: 

PTVL11 

FOR RH0VL11 | 


99 

ITERATE: 

PTVL12 

FOR RH0VL12 ; 


100 

ITERATE: 

PTVL13 

FOR RB0VL13 i 


101 

ITERATE: 

PTVL14 

FOR RH0VL14 j 


in? 

ITERATE: 

PTVL18 

FOR RHDVL1B j 


103 

IIERAIE: 

PTVLl G 

FOR RHDVL1G i 


104 

ITERATE: 

PTVLl 7 

FOR RH0VL17 j 


108 

ITERATE: 

PTVLl 8 

FOR RH0VL18 j 


106 

ITERATE: 

PTVLl 9 

FOR RHOVL10 j 


107 

ITERATE: 

PTVL2 

FOR RR0VL2 j 


108 

ITERATE: 

PTVL20 

FOR RH0VL20 j 


109 

ITERATE: 

PTVL21 

FOR RH0VL21 j 


110 

ITERATE: 

PTVL22 

FOR RH0VL22 j 


111 

ITERATE: 

PTVL3 

FOR RR0VL3 j 


112 

ITERATE: 

PTVL4 

FOR RH0VL4 j 


113 

ITERATE: 

PTYLS 

FOR RR0VL5 | 


114 

ITERATE: 

PTVL8 

FOR RH0VL6 j 


118 

ITERATE: 

PTVL7 

FOR RH0VL7 j 


116 

ITERATE: 

PTVL8 

FOR RR0VL8 » 


117 

ITERATE: 

PTVL9 

FOR RH0VL9 ; 


118 

END INTEGRATION 



119 



:*** 


120 

* SET-UP BALANCES 

* 


121 

*********************±* 


122 

DEFINE BALANCES 



123 

BALANCE 

WCCVBAL 

: WCCV UNTIL PTVL9 

- PTVL9C j 

124 

BALANCE 

WF9BAL 

: WF9 UNTIL PTVL10 

- STBL0001 j 

128 

BALANCE 

WF10BAL 

: WF10 UNTIL PTVL11 

- STBL0002 j 

126 

BALANCE 

WHC2BAL 

: WHG2 UNTIL PTOPRB 

- STBL0003 | 

127 

BALANCE 

THTODBAL 

: TTHTOD UNTIL TTHTOD 

- SYBL0004 j 

128 

BALANCE 

WHCGBAL 

: WHGG UNTIL PTQSF 

- STBLOOOS f 

129 

BALANCE 

WHC1BAL 

: WHC1 UNTIL PTFPRB 

- STBL0006 f 

130 

BALANCE 

THTFDBAL 

: TTHTFD UNTIL TTHTFD 

- STBLQ007 } 

131 

BALANCE 

WHCGBAL 

: WHGG UNTIL PTFSF 

- STBLOOOS | 

132 

BALANCE 

WFINJBAL 

: WFINJ UNTIL PTMFI 

- STBLOOOS | 

133 

BALANCE 

WHPFPBAL 

: WHPFP UNTIL PTVL3 

- PTHPFD j 

134 

BALANCE 

WHPOPBAL 

: WHPOP UNTIL PTVL21 

- PTHPOD j 

138 
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TIBS V/O CONTROL 


PACE 4 


CONFIG. INPOT 


BALANCE ILPFPBAL : WLPFP UNTIL PTVL1 - PTLPFD l 136 

BALANCE WLPOPBAL » WLPOP UNTIL PIVL19 - PTLPOD j 137 

BALANCE VPRBPBAL t WPRBP UNTIL PTPBSO - PIPBPD j 138 

END BALANCES 139 

*********************************************************************** 140 

*********************************************************************** 141 

* 142 

DEFINE CONFIGURATION * 143 

* 144 


*********************************************************************** 146 

*********************************************************************** 146 


* 147 

• 148 

ABOVE THE ITERATION LOOP * 149 

• 160 

. * 151 

, 162 

DEFINE SYSTEM ABOVE 153 

> 184 

t ABOVE THE ITERATION LOOP MODULES AND EQUATIONS GO HERE 156 

. 186 

‘END SYSTEM ABOVE 167 

i 168 

* 150 

i * 160 

> INSIDE THE ITERATION LOOP * 161 

> * 162 

* 163 

t 164 

DEFINE SYSTEM INSIDE 168 

k 166 

********************************************************************** 167 

i HYDROGEN PROPERTIES * 168 

*********************************************************************** 169 

PROPERTY PACKAGES H2PR0P; 170 

LOCATION HTNK t RHO-FIPT.KT), TT-F(PT,HT)j 171 

LOCATION VL1 t RHO-F(PT,HT> , TT-F(PT,HT)| 172 

LOCATION VL2 i RHO-F(PT.HT) . TT-F(PT,HT)| 173 

LOCATION VL3 t RHO-F(PT,HT), TT-F(PT,HT>, CAMA-F(HT,PT) ; 174 

LOCATION VL4 i RHO-F(PT.HT) , TT-F(PT,HT>j 176 

LOCATION VL8 i RHO-F(PT,HT) , TT-F(PT,HT>j 176 

LOCATION VLB * RHO-F(PT.HT) , TT-F(PT,HT>, 177 

CP-F(HT,PT>, MU-F(HT,PT>, K-F(HT.PT); 178 

LOCAHON VL7 t RHO-F(PT.HI) , TT-F(PT,HT), 179 

CP-F(HT.PT), MU-F(HT,PT>, K-F(HT.PT) j 180 


TIRE V/O CONTROL 


PACE 8 


CONFIG. INPOT 


LOCATION VL 8 i RHO-F(PT.HT) , TT-F(PT,HT)j 181 

LOCAnON VL9 : RHO-F(PT.HT), TT-F(PT,HT>» 182 

LOCATION VLIO t RHO-F(PT.HT) , TT-F(PT,HT)j 183 

LOCATION VL11 i RHO-F(PT.HT) , TT-F(PT,BT>, 184 

CP-F(HT.PT), M0-F(HT,PT), K-FCHI,PT>i 185 

LOCATION VL12 : RHD-F(PT,HT>, TT-F(PT,HT), S-FIHT.PT>, 186 

CP-F(HT.PT), MD-F(HT.PT), K-F(HT.PT>| 187 

LOCATION VL13 i RHO-F(PT,HT) , TT-F(PT,BT)s 188 

LOCATION VL14 s RHO-F(PT,HT> , TT-F(PT,HT>| 189 

LOCATION VL18 : RHO-F(PI,OT>. IT-F(PT,HI>s 190 

LOCATION VL16 : RHO-F(PT,HT> . TT-F(PT.HT), CAMA-F(HT.PT) i 191 

LOCAnON PBSF t RHO-F(PT.HT) , TT-FtPT.HT), CAMA-F(HT.PT) 5 192 

END PROPERTY 193 

* OXYGEN PROPERTIES * 198 

*********************************************************************** 196 

PROPERTY PACKAGE: 02 PR 0 P; 197 

LOCATION OINK : RHO-F<PT,HT) , TT-F(PT,HT) ; 198 

LOCATION 7L17 t RHO-F(PT.ET) , TT-F(PT,HT)j 199 

LOCATION YL18 t RHO-F(PT.HT) , TT-F(PT,HT)j 200 

LOCATION VL19 « RHO-F(PT,KT) t TT-F(PT,HT)j 201 

LOCATION VL20 : RHO-F(PT.ET) , TI-F(PT,HT) ; 202 

LOCATION 7L21 i RHO-F(PT,ET> , TT-F(PT,HT)j 203 

LOCATION VL22 : RHO-F(PT.BT) , TT-F(PT.HT). S-F<ET,PT>j 204 

LOCATION PESO : RBO-F(PT.HT), TT-F(PT,HT)j 208 

END PROPERTY 206 

*********************************************************************** 207 

* HOT GAS PROPERTIES * 208 

*********************************************************************** 209 

EQUATION: PTHTFD - PTFSF t 210 

EQUATION: OFRHIFD - OFRFTBP t 211 

EQUATION: HFRHIFD - HFRFTBP j 212 

EQUATION: PIHTOD - PTOSF j 213 

EQUATION: OFRHTOD - OFROTBP j 214 

EQUAHON: HFRHTOD - HFROTBP ; 218 

EQUATION: 0FRVL16 - 0.0 * 216 

EQUATION: HFR7L16 - 0.0 t 217 

EQUATION: RCASAMB - 640.0 : 218 

EQUATION: GAMAAMB - 1.4 I 219 

EQUATION: OFRAMB - 0.0 j 220 

ELATION: HFRAMB - 0.0 ; 221 

EQUATION: OFRPBSF - 0.0 i 222 

HQUAHON: OFRPBSF - 0.0 » 223 

* 224 

PROPERTY PACKAGE: HCPROPj 228 
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roe f/0 CONTROL 


PACE 6 


CONFIC. INPOT 


LOCATION MCHB 

LOCATION FPRB 
LOCATION OPRB 
LOCATION FTBP 

LOCATION OTBP 

LOCATION HTFD 
LOCATION HTOD 
LOCATION FSF 
LOCATION OSF 
LOCATION MFI 
END PROPERTY 


i K-F(PT,TT) , MU-F{PT,IT) ,RHO-F(PT,TT),CAMA-F(PT,IT) , 


R-F(PT.TT), CP-F(PT,TI)j 
jCAMA-F(PT.TT), R-F(PT.TT), CP-F(PT,TT), 
jCAMA-F(PT.TT), R-F(PT.TT), CP-F(PT,TT>, 
«CAMA-F(PT,TT>, R-F(PT,TT>, CP-F(PT.TT), 


RHO-F(PT r TT) ! 
»CAMA-F(PT,TT), 
RHO-F(PT,TT) j 
* CP-F(PT,IT), 
: CP-F(PT.IT). 

jGAMA-F(PT.TT), 
»CAMA-F(PT,IT>, 
:GAMA-F(PT,TT), 


R-F(PT,TT), CP-F(PT,TT), 

CAMA-F(PT,TT)} 
GAMA-F(PT,TT)j 
R-F(PT,TT>, CP-F(PT.TT), 
R-F(PT.TT), CP-F(PT.TT), 
R-F(PT,TT) , CP-FIPT.IT>, 


RHD-F(PT.IT)t 
RHO-F(PT,TT) i 
Z-F(PT.TT) . 

Z-F(PT,TT) , 


RHO-F(PT,TT> •, 
RHO-F(PT,TT>: 
RHO-F(PT,TT)j 


ft EQUATIONS FOR PROPERTIES * 

*********************************************************************** 

* 


* INTERNAL ENERGIES - FUEL SIDE 


EQUATION: 

UTHTNK 

- 

UTHTNK 

- 

(l./RJ) 

* PTHTNK 

/ 

RHOHTNK 

EQUATION: 

UTYL1 

- 

HTYL1 

- 

( l./RJ) 

* PTVL1 

/ 

RH0VL1 

EQUATION: 

UTVL2 

- 

HTVL2 

- 

(l./RJ) 

* PTVL2 

/ 

RH0VL2 

EQUATION: 

UT7L3 

- 

BTVL3 

- 

(l./RJ) 

* PTVL3 

/ 

RH0VL3 

EQUATION: 

UTVL4 

- 

HTYL4 

m 

(l./RJ) 

* PTVL4 

/ 

RH0VL4 

EQUATION: 

UTVL8 

- 

BTVL8 

- 

(l./RJ) 

* PTVL8 

/ 

RH0VL8 

EQUATION: 

UTVL6 

- 

HTVL6 

- 

(l./RJ) 

* PTVL6 

/ 

RH0TL6 

EQUATION: 

UTVL7 

- 

HTVL7 

- 

(l./RJ) 

* PTVL7 

/ 

RH0VL7 

EQUATION: 

UTVL8 

- 

SIYL8 

- 

(l./RJ) 

* PTVL8 

/ 

RH0VL8 

EQUATION: 

UTVL9 

- 

HTVLfl 

- 

(l./RJ) 

ft PTVL9 

/ 

RH0VL9 

EQUATION: 

UTVL10 

- 

HTVL10 

- 

(l./RJ) 

* PTTL10 

/ 

RHOVLIO 

EQUATION: 

UTVL11 

- 

HTVL11 

- 

(l./RJ) 

* PTVL11 

/ 

RH0VL11 

EQUATION: 

UTVL12 

- 

HTVL12 

- 

(l./RJ) 

* PT7L12 

/ 

RH0VL12 

EQUATION: 

UTVL13 

- 

HTVL13 

- 

(l./RJ) 

* PTVL13 

/ 

RH0VL13 

EQUATION: 

CTVL14 

- 

HTVL14 

- 

(l./RJ) 

* PTVL14 

/ 

RH0VL14 

EQUATION: 

UTVL18 

- 

HTYL15 

- 

(l./RJ) 

* PTVL18 

/ 

RH0VL1B 

EQUATION: 

UTVL16 

- 

ETTL16 

- 

(l./RJ) 

* PTVL16 

/ 

RH0VL16 

ROTATION: 

UTPBSF 

- 

BIPBSF 

- 

(l./RJ) 

* PTPBSF 

/ 

RHOPBSF 

* INTERNAL ENERGIES 

- OXYGEN SIDE 




EQUATION: 

UTOTNK 

- 

HTDTNK 

- 

(l./RJ) 

* PTOINK 

/ 

RHOQTNK 

EQUATION: 

UTVL17 

- 

HTYL17 

- 

(l./RJ) 

* PTVL17 

/ 

RH0VL17 

EQUATION: 

UT7L18 

- 

BTVL18 

- 

(l./RJ) 

* PTVL18 

/ 

RH0VL18 

EQUATION: 

UT7L19 

- 

HTVL19 

- 

(l./RJ) 

* PTVL19 

/ 

RHOVLIO 

EQUATION: 

UT7L20 

- 

BTVL20 

- 

(l./RJ) 

* PTVL20 

/ 

RH0VL20 

EQUATION: 

UTVL21 

- 

HTVL21 

- 

(l./RJ) 

* PTVL21 

/ 

RH0VL21 

EQUATION: 

UTVL22 

- 

HTVL22 

- 

(l./RJ) 

* PTVL22 

/ 

RH0VL22 


226 

227 

228 

229 

230 

231 

232 

233 

234 
238 

236 

237 

238 

239 

240 

241 

242 

243 

244 
248 

246 

247 

248 

249 
280 
281 
282 

283 

284 
288 
286 

287 

288 
289 
260 
261 
262 

263 

264 
268 
266 

267 

268 

269 

270 
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TIDE V/0 CONTROL 


PACE 7 


CONFIG. INPUT 


EQUATION t CTPBSO - HTPBSO - (l./RJ) * PTPBSO / RHOPBSO j 271 

* ENTHALPIES - HOT GAS SIDE 272 

EQUATION: HTFPRB - CPFPRB * TTFPRB j 273 

EQUATION: HTOPRB - CPOPRB * TTOPRB ; 274 

EQUATION: HTFTBP - CPFTBP * IIFTBP l 275 

EQUATION: HTOIBP - CPOTBP * TTOTBP j 276 

EQUATION: HTFSF - CPFSF * ITFSF | 277 

EQUATION: HTOSF - CPOSF * TTOSF j 27S 

EQUATION: HIMFI - CPMFI * TTMFI j 279 

EQUATION: HIMCHB - CPMCHB * TIMCHB ( 280 

EQUATION: HIHTFD - CPHTFD * TTHIFD : 281 

EQUATION : KTHTOD - CPHTOD * TTHTOD j 282 

*********************************************************************** 283 

t FUEL SIDE NON-DERIVATIVE MODULES * 284 

*********************************************************************** 285 

* 286 

EQUATION : AREAFPOV - (.337998 / .32238 ) * AREAFPV j 287 

EQUATION : AREAOPOV - (.110886 / .11835 ) * AREAOPV j 288 

***************************** 289 

* LPFP EXIT DENSITY * 290 

***************************** 291 

EQUATION: RHOLPFD - RHOVLlj 292 

********************************** 293 

* LOW PRESSURE FUEL PUMP * 294 

********************************** 295 

MODULE: PUMPOlj 296 

NAME: LPFP j 297 

I/O LIST: INLET FLOW - LPFP, 298 

INLET PROPERTIES - HINK, 299 

EXIT PROPERTIES - LPFD, 300 

SHAFT - FL i 301 

DESIGN VALUES: SND - 15603.4, 302 

TRQD - 12924.8, 303 

WD - 148.7, 304 

HDD - 101574.6, 305 

GEAR - 1.0* 306 

MAP: PMAPOS* 307 

CMT: LOW PRESSURE FUEL PUMP* 308 

END MODULE 309 

***************************** 310 

* HPFP EXIT DENSITY * 311 

***************************** 312 

EQUATION: RHDHPFD - RH0VL3; 313 

*********************************** 314 

* HIGH PRESSURE FUEL PUMP • 315 


184 



ITBE I/O CONTROL 


PAGE 8 


CONFIG. INPUT 


*********************************** 316 

MODULE: PUMPOl; 317 

NAME: HPFP; 318 

I/O LIST: INLET FLOW - HPFP, 319 

INLET PROPERTIES - VL2 , 320 

EXIT PROPERTIES - HPFD, 321 

SEAFT - FH i 322 

DESIGN VALUES: SND - 34189.8, 323 

IRQD - 110141.9, 324 

WD - 148.7, 325 

HDD - 2229273.6, 326 

GEAR - 1.0: 327 

MAP: PMAP04 j 328 

CMT: HIGH PRESSURE FUEL PUMP: 329 

END MODULE 330 

************************************************** 331 

* N ON-INERTIAL FUEL TURBINE COOLING LINE * 332 

************************************************** 333 

MODULE: PIPE01; 334 

NAME: FTC j 338 

I/O LIST: INLET PROPERTIES - VL3 , 336 

EXIT PROPERTIES - FSF l 337 

DESIGN VALUES: CF - 2.025: 338 

CMT: FUEL TURBINE COOLING FLOW FROM VOLUME 3 TO VOLUME FSF; 339 

END MODULE 340 

*************************** 341 

* MAIN FUEL VALVE * 342 

*************************** 343 

MODULE: VALVOO; 344 

NAME: MFV i 348 

I/O LIST: UPSTREAM PROPERTY - VL3 , 346 

DOWNSTREAM PROPERTY - VL4 : 347 

DESIGN VALUES: AREA - 15.35313, 348 

RKLS - 2.289; * 349 

CMT: MAIN FUEL VALVE: 380 

END MODULE 381 

****************************************** 352 

* FUEL IGNITER NON-INERTTAL LINE * 383 

****************************************** 354 

MODULE: PIPE01: 358 

NAME: FIG ; 386 

I/O LIST: INLET PROPERTIES - VL4 . 387 

EXIT PROPERTIES - MCHBj 388 

DESIGN VALUES: CF - 0.44} 389 

CMT: FUEL IGNITER FLOW FROM VOLUME 4 TO THE MAIN CHAMBER; 360 
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END MODULE 361 

**************************************** 362 

* NON-INERTIAL FUEL LINE THREE * 363 

**************************************** 364 

MODULE* PIPEOlj 368 

NAMEt F3 i 366 

I/O LIST: INLET PROPERTIES - VLB , 367 

EXIT PROPERTIES - VL6 ? 368 

DESIGN VALUES! CF - 198. 6; 369 

CMTi PIPE FLOW FROM VOLUME 8 TO VOLUME 6} 370 

END MODULE 371 

*************************************** 372 

* NON-INERTIAL FUEL LINE FOUR * 373 

*************************************** 374 

MODULE: PIPE01; 378 

NAMEi F4 i 376 

I/O LIST! INLET PROPERTIES - VL6 , 377 

EXIT PROPERTIES - VL7 } 378 

DESIGN VALUES: CF - 133. 6| 379 

CMT: PIPE FLOW FROM VOLUME 6 TO VOLUME 7» 380 

END MODULE 381 

********************************* 382 

* COOLANT CONTROL VALVE * 383 

********************************* 384 

EQUATION i RH0VL9C - RH0VL9 } 388 

MODULE! PIPE04s 386 

NAME! CCV » 387 

I/O LIST! UPSTREAM PROPERTY - VL9C, 388 

DOWNSTREAM PROPERTY - VL8 | 389 

DESIGN VALUES! RKLS - 1.763 j 390 

CMTi COOLANT CONTROL VALVE} 391 

END MODULE 392 

*************************************** 393 

* NON-INERTIAL FUEL LINE NINE * 394 

*************************************** 398 

MODULE! PIPE06} 396 

NAMEt PS } 397 

I/O LIST! INLET PROPERTIES - VL10, 398 

EXIT PROPERTIES - VL11; 399 

nrcroi wir iipc. rv — 9. 400 

CMI: UPSTREAM PRESSURE CALC FROM VOLUME 11 TO VOLUME 10} 401 

END MODULE *02 

************************************** 403 

* NON-INERTIAL FUEL LINE TEN * 404 

************************************** 408 
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MODULE! PIPE06 j 408 

NAME! F10 j 407 

I/O LISTs INLET PROPERTIES - VL11, 408 

EXIT PROPERTIES - VL12; 409 

DESIGN VALUES! CF - 29.3? 410 

CMT: UPSTREAM PRESSURE CALC FROM VOLUME 12 TO VOLUME 11; 411 

END MODULE 412 

******************************************************** 413 

* LOW PRESSURE FUEL TURBINE DISCHARGE PRESSURE * 414 

******************************************************** 415 

EQUATION! PTLTFD - PTVL13; 416 

************************************* 417 

* LOW PRESSURE FUEL TURBINE * 418 

************************************* 419 

MODULE! IURB02; 420 

NAME! LPFT; 421 

I/O LIST! INLET PROPERTIES - VL12, 422 

EXIT PROPERTIES - LTFD, 423 

SHAFT WORK - FL ; 424 

DESICN VALUES: ETAD - 0.53. 428 

PRD - 1.33, 426 

SND - 1S603.4 . 427 

AREA - 1.0 , 428 

DIAM - 6.63 . 429 

DC1 0 , 430 

DC2 - 0 . 431 

GEAR - 1.0 ; 432 

MAPt TBMP04; 433 

CMT! LOW PRESSURE FUEL TURBINE; 434 

END MODULE 438 

***************************************** 438 

* NON -INERTIAL FUEL LINE ELEVEN * 437 

***************************************** 438 

MODULE! PIPE01; 439 

NAME! Fll ; 440 

I/O LIST! INLET PROPERTIES - VL13. 441 

EXIT PROPERTIES - VL14; 442 

DESICN VALUES! CF - 108.6; 443 

CMT! PIPE FLOW FROM VOLUME 13 TO VOLUME 14; 444 

END MODULE 448 

***************************************** 448 

* NON -INERTIAL FUEL LINE TWELVE * 447 

***************************************** 448 

MODULE! PIPE01; 449 

NAME; F12 t 480 


187 




TIBS I/O CONTROL 


PACE 11 


CONFIG. INPUT 


I/O LIST i INLET PROPERTIES - VL14. 4S1 

EXIT PROPERTIES - VL18; 402 

DESIGN VALUES) CF - 147.4; 403 

CMTt PIPE FLOW FROM VOLUME 14 TO VOLUME IS; 404 

END MODULE 400 

******************************************* 406 

* NON-INERTIAL FUEL LINE THIRTEEN * 407 

******************************************* 408 

MODULE) PIPE01; 409 

NAME) F13 ; 460 

I/O LIST) INLET PROPERTIES - VL13, 461 

EXIT PROPERTIES - VL10; 462 

DESIGN VALUES) CF - 143.0; 463 

CMTs PIPE FLOW FROM VOLUME 13 TO VOLUME 10; 464 

END MODULE 460 

466 

* NON-INERTIAL FUEL LINE FOURTEEN * 467 

******************************************* 468 

MODULE) PIPED 1; 469 

NAME) F14 ; 470 

I/O LIST: INLET PROPERTIES - VL1B, 471 

EXIT PROPERTIES - VL16; 472 

DESIGN VALUES) CF - 199.4; 473 

CMTt PIPE FLOW FROM VOLUME 10 TO VOLUME 16; 474 

END MODULE 470 

****************************************** 476 

* NON-INERTIAL FUEL LINE FIFTEEN * 477 

****************************************** 478 

MODULE) PIPE01; 479 

NAME) FIB ; 480 

I/O LIST) INLET PROPERTIES - VL16. 481 

EXIT PROPERTIES - MCHBj 482 

DESIGN VALUES) CF - 76.9; 483 

CMTt PIPE FLOW FROM VOLUME 16 TO THE MAIN CHAMBER; 484 

END MODULE 480 

******************************************************** 486 

* NON-INERTIAL FUEL LINE TO MAIN FUEL INJECTOR * 487 

******************************************************** 488 

MODULE) PIPE01; 489 

NAME) FSLV; 490 

I/O LIST: INLET PROPERTIES - VL16, 491 

EXIT PROPERTIES - MFI ; 492 

DESIGN VALUES: CF - 28.2; 493 

OCT: PIPE FLOW FROM VOLUME 16 TO THE MAIN FUEL INJECTOR; 494 

END MODULE 490 
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****************************************************** 496 

* NON -INERTIAL OXIDIZER TURBINE COOLING LINE * 497 

****************************************************** 498 

MODULE i PIPEOli 499 

NAME: QIC i BOO 

I/O LIST: INLET PROPERTIES - PBSF, 801 

EXIT PROPERTIES - OSF s 802 

DESICN VALDES: CF - 0.661: 803 

CUT: OXYGEN TURBINE COOLING FLOW FROM VOLUME PBSF TO VOLUME FSF; 804 

END MODULE 808 

* MANIFOLD COOLING HEAT IRANSFER 806 

EQUATION: TKMCO - .00167 : 807 

EQUATION: TKMCF - .00131 : 803 

EQUATION: AHIMCO - 728.0 : 809 

EQUATION: AHIMCF - 872.0 j 810 

EQUATION: QDOTFMCI -TKMCO* AHIMCO* (TTMFI-ITVL 13)* 811 

SQRKAMAXKO. , (1F13/U.36))): 812 

EQUATION: QDOTFMCO - - QDOTFMCI | 813 

EQUAnON: QDOTCMCI -IEMCF*AHTMCF*(TIMFI-TTVL13)* 814 

SQRT(AMAX1(0.,(WF11/18.80))): 818 

EQUAHON: QDOTOMCO - - QDOTCMCI j 816 

*********************************************************************** 817 

* FUEL SIDE DERIVATIVE MODULES • 818 

*********************************************************************** 819 

* 820 

********************************** 821 

• INERTIAL FUEL LINE ONE * 822 

********************************** 823 

MODULE: PIPEOO; 824 

NAME: FI | .828 

I/O LIST: INLET PROPERTIES - VL1 . 826 

EXIT PROPERTIES - VL2 i 827 

DESIGN VALUES: CF - 699.2, 828 

AREA - 28.26, 629 

XLHN - 20.00: 830 

CMT: FLOW DERIVATIVE FROM VOLUME 1 TO VOLUME 2: 831 

END MODULE 832 

********************************** 833 

* INERTIAL FUEL LINE TWO * 834 

********************************** 838 

MODULE: PIPEOO: 836 

NAME: F2 \ 837 

I/O LIST: INLET PROPERTIES - VL4 , 838 

EXIT PROPERTIES - VLB j 839 

DESIGN VALUES: CF - 190.1, 840 
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AREA - 10.60, 841 

XLEN - 17.80| 842 

CMTi FLOE DERIVATIVE FROM VOLUME 4 TO VOLUME 8; 843 

END MODULE 844 

********************************** 848 

* INERTIAL FUEL LINE SIX * 846 

********************************** 847 

MODULE i PIPEOOj 848 

NAME: F6 | 849 

I/O LIST: INLET PROPERTIES - VL4 , 880 

EXIT PROPERTIES - VL9 } 881 

DESIGN VALUES: CF - 79.6, 882 

AREA - 1.00, 883 

XLEN - 4.027{ 884 

CMT: FLOW DERIVATIVE FROM VOLUME 4 TO VOLUME 9: 888 

END MODULE 886 

************************************ 887 

* INERTIAL FUEL LINE BIGHT * 888 

************************************ 889 

MODULE: PIPEOOj 860 

NAME: F8 j 861 

I/O LIST: INLET PROPERTIES - VL4 , 862 

EXIT PROPERTIES - VLlOj 863 

DESIGN VALUES: CF - 68.2, 864 

AREA - 3.14, 868 

XLEN -134.00} 866 

CMT: FLOW DERIVATIVE FROM VOLUME 4 TO VOLUME 10} 867 

END MODULE 868 

*********************************** 569 

* INERTIAL FUEL LINE FIVE * 870 

*********************************** 871 

MODULE: PIPEOOi 872 

NAME: F8 } 873 

I/O LIST: INLET PROPERTIES - VL7 . 874 

EXIT PROPERTIES - VL8 } 878 

DESIGN VALUES: CF - 142.3, 876 

AREA - 17.86, 877 

XLEN - 92.00} 878 

CMT: FLOW DERIVATIVE FROM VOLUME 7 TO VOLUME 8} 879 

END MODULE 880 

************************************ 581 

* INERTIAL FUEL LINE SEVEN * 882 

************************************ 583 

MODULE: PIPEOOs 884 

NAME: F7 } 888 
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I/O LISTt INLET PROPERTIES - VL8 , 886 

EXIT PROPERTIES - PBSF; 887 

DESIGN VALUES: CF - 682.1, 888 

AREA - 4.73 , 889 

XLEN - 77.80| 890 

CMT: FLOW DERIVATIVE FROM VOLUME 8 TO VOLUME PBSFi 891 

END MODULE 892 

**************************************************** 893 

* INERTIAL FUEL LINE TO THE FUEL PREBURNER * 894 

**************************************************** 898 

MODULE: PIPEOO; 896 

NAME: FFPB: 897 

I/O LIST: INLET PROPERTIES - PBSF, 898 

EXIT PROPERTIES - FPRB; 899 

DESIGN VALUES: CF - 124.0, 600 

AREA - 1.00, 601 

XLEN - 80.00} 602 

CMT: FLOW DERIVATIVE FROM VOLUME PBSF TO VOLUME FPRB; 603 

END MODULE 604 

******************************************************** 608 

* - — INERTIAL FUEL LINE 10 THE OXIDIZER PREBURNER * 606 

******************************************************** 607 

MODULE: PIPEOO; 608 

NAME: FOPBs 609 

I/O LIST: INLET PROPERTIES - PBSF. 610 

EXIT PROPERTIES - OPRB; 611 

DESIGN VALUES: CF - 60.8, 612 

AREA - 1.00, 613 

XLEN - 80.00; 614 

CMT: FLOW DERIVATIVE FROM VOLUME PBSF TO VOLUME OPRB; 618 

END MODULE 616 

********************** 617 

* VOLUME ONE * 618 

********************** 619 

MODULE: VOLMOO; 620 

NAME: VL1 ; 621 

I/O LIST: UPSTREAM PROPERTIES - LPFD, 622 

INLET FLOW - LPFP, 623 

EXIT FLOW - FI . 624 

DOWNSTREAM PROPERTIES - VL2 , 628 

QDOT - VL1 ; 626 

DESIGN VALUES: VOL - 368.2, 627 

QDOT - 0.0 ; 628 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 1; 629 

END MODULE 630 
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********************** 631 

* VOLUME TWO * 632 

********************** 633 

MODULE: VOLMDD; 634 

NAME: VL2 j 638 

I/O LIST: UPSTREAM PROPERTIES - VL1 , 636 

INLET FLOW - FI , 637 

EXIT FLOW - HPFP, 638 

DOWNSTREAM PROPERTIES - HPFD, 639 

QDOT - VL2 •, 640 

DESIGN VALUES: VOL - 200.0, 641 

QDOT - 0.0 i 642 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 2; 643 

END MODULE . 644 

************************ 648 

* VOLUME THREE * 646 

************************ 647 

MODULE: V0LM01: 648 

NAME: VL3 j 649 

I/O LIST: UPSTREAM PROPERTIES - HPFD. 650 

INLET FLOW - HPFP, 651 

EXIT FLOW - MFV , FTC , 652 

DOWNSTREAM PROPERTIES - VL4 , FSF , 663 

QDOT - VL3 ; 654 

DESIGN VALUES: VOL - 347.9; 655 

Off: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 3} 656 

END MODULE 657 

*********************** 658 

* VOLUME FOUR * 669 

*********************** 660 

MODULE: VOLMOlj 661 

NAME: VL4 \ 662 

I/O LIST: UPSTREAM PROPERTIES - VL3 , 663 

INLET FLOW - MFV , 664 

EXIT FLOW - F2 . F6 , F8 , FIG . 665 

DOWNSTREAM PROPERTIES - VL8 . VL9 . VL10, MCHB, 666 

JJDOT - VL4 | 667 

DESICN VALUES: VOL - 186.0: 668 

Off: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 4; 669 

END MODULE 670 

*********************** 671 

* VOLUME FIVE * 672 

*********************** 673 

MODULE: VOLMOO: 674 

NAME: VL5 ; 675 
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I/O LIST: UPSTREAM PROPERTIES - VL4 , 676 

INLET FLOW - F2 , 677 

EXIT FLOW - F3 , 678 

DOWNSTREAM PROPERTIES - VL6 . 679 

QDOT - VLB t 680 

DESIGN VALUES: VOL - 866.3, 681 

(JDOT - 0.0 t 682 

CMTt DENSITY AND INTERNAL ENERCY DERIVATIVES FOR VOLUME 6: 683 

END MODULE 684 

************************************* 68B 

* NOZZLE COOLING VOLUME SIX * 688 

************************************* 687 

MODULE t NCLVOOt 688 

NAME: VL6 | 689 

I/O LISTs UPSTREAM PROP - VLB , 690 

INLET FLOW - F3 , 691 

EXIT FLOW - F4 , 692 

DOWNSTREAM PROP - VL7 , 693 

(JDOT - 6 HOT, 6AMB, 694 

METAL TEMPERATURES - MTL1, MTLSj 69S 

DESIGN VALUES: VOL -3160.0, 696 

AREA - 17.872 j 697 

CMTt DENSITY AND INTERNAL ENERCI DERIVATIVES FOR VOLUME 6s 698 

END MODULE 699 

*************************************** 700 

* NOZZLE COOLING VOLUME SEVEN * 701 

*************************************** 702 

MODULES NCLVOOt 703 

NAME: VL7 | 704 

I/O LISTs UPSTREAM PROP - VL8 , 70S 

INLET FLOW - F4 , 706 

EXIT FLOW - FB , 707 

DOWNSTREAM PROP - VL8 , 708 

(JDOT - 7H0T, 7AMB, 709 

METAL TEMPERATURES - MTL2, MILS; 710 

DESIGN VALUES: VX -1616.0, 711 

AREA - 40.376s 712 

CMTt DENSITY AND INTERNAL ENERCY DERIVATIVES FOR VOLUME 7: 713 

END MODULE 714 

************************ 71S 

* VXUMB EIGHT * 716 

************************ 717 

MODULE: VOLVO lj 718 

NAME: VL8 ; 719 

I/O LISTs UPSTREAM PROPERTIES - VL7 . VL9 , 720 
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INLET FLOW - FS , CCV , 721 

EXIT FLOW - F7 , 722 

DOWNSTREAM PROPERTIES - PBSF, 723 

QDOT - VLS i 724 

DESIGN VALUESt VOL - 1000. j 726 

Off t DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 8; 726 

END MODULE 727 

*********************** 728 

* VOLUME NINE * 720 

*********************** 730 

MODULE: VOLMOO} 731 

NAME: VLB * 732 

I/O LIST: UPSTREAM PROPERTIES - VL4 , 733 

INLET FLOW - F6 . 734 

EXIT FLOW - CCV , 736 

DOWNSTREAM PROPERTIES - VL8 , 736 

({DOT - VL9 j 737 

DESIGN VALUES: VOL - 500.0, 738 

({DOT - 0.0 j 739 

CUT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 9; 740 

END MODULE 741 

********************** 742 

* VOLUME TEN • 743 

********************** 744 

MODULE: VOLMOO} 746 

NAME: VL10} 746 

I/O LIST: UPSTREAM PROPERTIES - VL4 . 747 

INLET FLOW - F8 , 748 

EXIT FLOW - F9 , 749 

DOWNSTREAM PROPERTIES - VL11, 760 

ODDI - VLlDi 761 

DESIGN VALUES: VOL - 1000., 762 

({DOT - 0.0 } 763 

CMTj DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 10} 764 

END MODULE 768 

***************************************** 766 

* CHAMBER COOLING VOLUME ELEVEN * 787 

***************************************** 768 

MODULE: NCLVOO} 769 

NAME: VL11} 760 

I/O LIST: UPSTREAM PROP - VL10, 761 

INLET FLOW - F9 , 762 

EXIT FLOW - F10 , 763 

DOWNSTREAM PROP - VL12, 764 

({DOT - 11HT, HAM, 768 
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METAL TEMPERATURES - MTL3, MIL7; 766 

DESIGN VALUES: VOL - 144.2 , 767 

AREA - 10.927} 768 

CMT: DENSITY AND INTERNAL BNERCT DERIVATIVES FOR VOLUME 11: 769 

END MODULE 770 

***************************************** 771 

* CHAMBER COOLING VOLUME TWELVE * 772 

***************************************** .773 

MODULE: NCLVOO; 774 

NAME: VL12; 776 

I/O LIST: UPSTREAM PROP - VL11, 776 

INLET FLOW - F10 , 777 

EXIT FLOW - LPFT, 778 

DOWNSTREAM PROP - LTFD, 779 

QDOT - 12HT, 12AM, 780 

METAL TEMPERATURES - MTL4, MTL8; 781 

DESIGN VALUES: VOL - 144.2, 782 

AREA - 10.927: 783 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 12: 784 

END MODULE 788 

*************************** 786 

* VOLUME THIRTEEN * 787 

*************************** 788 

MODULE: V0LM01: 789 

NAME: VL13: 790 

I/O LIST: UPSTREAM PROPERTIES - LTFD, 791 

INLET FLOW - LPFT. 792 

EXIT FLOW - Fil , F13 , 793 

DOWNSTREAM PROPERTIES - VL14, VL18, 794 

QDOT - VL13; 796 

DESIGN VALUES: VOL - 600.0; 796 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 13; 797 

END MODULE 798 

*************************** 799 

* VOLUME FOURTEEN * 800 

*************************** 801 

MODULE: VOLMOO; 802 

NAME: VL14; 803 

I/O LIST: UPSTREAM PROPERTIES - VL13, 804 

INLET FLOW - Fll . 808 

EXIT FLOW - F12 , 808 

DOWNSTREAM PROPERTIES - VL16, 807 

QDOT - OMCI; 808 

DESIGN VALUES: VOL - 600.0; 809 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 14; 810 


195 



UBS I/O CONTROL 


PACE 19 


CONFIG. INPUT 


END MODULE 811 

************************** 812 

* VOLUME FIFTEEN * 813 

************************** 814 

MODULEt VOLMDOt 8 IS 

NAME: VLlSj 816 

I/O LIST« UPSTREAM PROPERTIES - VL13, 817 

INLET FLOW - F13 . 818 

EXIT FLOW - F14 , 819 

DOWNSTREAM PROPERTIES - VL16, 820 

QDOT - FMCI; 821 

DESIGN VALUES! VOL - 800. Oj 822 

CMTi DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 18; 823 

END MODULE 824 

************************** 828 

* VOLUME SIXTEEN * 826 

************************** 827 

MODULE: V0LM01; 828 

NAME: VLlBj 829 

I/O LIST: UPSTREAM PROPERTIES - VL14, VL1S. 830 

INLET FLOW - F12 . F14 , 831 

EXIT FLOW - FSLV. F18 . 832 

DOWNSTREAM PROPERTIES - MFI . MCHB, 833 

QDOT - VL16; 834 

DESIGN VALUES! VOL - 800.0; 838 

CMTj DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 16; 836 

END MODULE 837 

****************************************** 838 

* PREBURNER FUEL SPLITTER VOLUME * 839 

****************************************** 840 

MODULE! V0LM01; 841 

NAME! PBSF; 842 

I/O LIST! UPSTREAM PROPERTIES - VL8 , 843 

INLET FLOW - F7 , 844 

EXIT FLOW - FFPB. FOPB. OTC . 848 

DOWNSTREAM PROPERTIES - FPRB. OPRB. OSF . 846 

qDOI - PBSF; 847 

DESIGN VALUES! VOL - 800.0; 848 

CMI: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME PBSF; 849 

END MODULE 880 


i******************************************************************** 881 

OXYCEN SIDE NON -DERIVATIVE MODULES * 882 

******************************************************************** 883 

884 

LPOP EXIT DENSITY 888 
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EQUATION] RHOLPOD - RH0VL19; 886 

************************************** 887 

* LOW PRESSURE OXIDIZER PUMP * 888 

************************************** 880 

MODULE: PUMP01; 860 

NAME: LPOP) 861 

I/O LIST: INLET FLOW - LPOP, 862 

INLET PROPERTIES - VL18, 863 

EXIT PROPERTIES - LPOD, 864 

SHAFT - OL i 868 

DESIGN VALUES: SND - 8041.8, 866 

TRQD - 18818.1, 867 

WD - 896.2, 868 

HDD - 7873.6, 869 

GEAR - 1.0; 870 

MAP: PMAP07; 871 

CMT: LOW PRESSURE OXYGEN PUMP; 872 

END MODULE 873 

* BPOP EXIT DENSITY 874 

EQUATION: RHOHPOD - RH0VL21; 878 

*************************************** 876 

* HIGH PRESSURE OXIDIZER PUMP * 877 

*************************************** 878 

MODULE: PUMP01; 879 

NAME: HP OP; 880 

I/O LIST: INLET FLOW - HP OP, 881 

INLET PROPERTIES - VL20 , 882 

EXIT PROPERTIES - HPOD, 883 

SHAFT - OH » 884 

DESIGN VALUES: SND - 27240.8, 888 

TRQD - 80843.4, 886 

WD - 1078.4, 887 

HDD - 91778.2, 888 

GEAR - 1.0; 889 

MAP: PMAP06; 890 

CMT: HIGH PRESSURE OXYGEN PUMP; 891 

END MODULE 892 

******************************************* 893 

* NON- INERTIAL OXIDIZER LINE FOUR * 894 

******************************************* 898 

MODULE: PIPE01; 896 

NAME: 04 ; 897 

I/O LIST: INLET PROPERTIES - VL21, 898 

EXIT PROPERTIES - VL20; 899 

DESIGN VALUES: CF - 0.221; 900 
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CMTt HP OP RECIRC. FLOW FROM VOLUME 21 TO VOLUME 20s 901 

END MODULE 902 

* LINE 7 DOWNSTREAM DENSITY 903 

EQUATION: RHOPOCO - RH0VL21; 904 

******************************************** 90S 

* NON- INERTIAL OXIDIZER LINE SEVEN * 906 

******************************************** 907 

MODULE: PIPEOls 908 

NAME: 07 ; 909 

I/O LIST: INLET PROPERTIES - VL21, 910 

EXIT PROPERTIES - POGOj 911 

DESICN VALUES: CF - 0.050 j 912 

CMT: LIQUID POCO FLOW FROM VOLUME 21 TO POCO; 913 

END MODULE 914 

************************************************************ 918 

* LOW PRESSURE OXIDIZER TURBINE DISCHARGE PRESSURE * 916 

************************************************************ 917 

EQUATION: PILTOD - PTVL19; 918 

***************************************** 919 

* LOW PRESSURE OXIDIZER TURBINE * 920 

***************************************** 921 

MODULE: ITJRB02; 922 

NAME: LPOTj 923 

I/O LIST: INLET PROPERTIES - VL22, 924 

EXIT PROPERTIES - LTOD, 928 

SHAFT WORE NODE - OL j 926 

DESIGN VALUES: ETAD - 0.64, 927 

PRD - 9.73, 928 

SND - 8041.8, 929 

AREA - 1.0 930 

DIAM - 0.0 . 931 

CEAR - 1.0 932 

DC1 - 0 . 933 

DC2 Os 934 

MAP: TEMPOS j 838 

Off: LOW PRESSURE OXYGEN TURBINE s 936 

END MODULE 937 

******************************* 938 

* MAIN OXIDIZER VALVE * 939 

******************************* 940 

* LUMP LINE, VALVE, AND INJECTOR RESISTANCES 941 

EQUATION: RKLSMOV- 0.881; 942 

EQUATION: CFMOVX - 37.98 * AREAMOV / SQRT(RKLSMOV) | 943 

EQUATION: CF09 - 729.8 ; 944 

EQUATION: CF90V - SQRT( CFMOVX* *2 * CF09**2 / (CFM0VX**2 ♦ CF09**2>>; 948 
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* SCALE INJECTOR RESISTANCE WITH PRIMING FRACTION 946 

EQUATION: CFOIJ - CFOINJ / PRIMMOI ; 947 

EQUATION: CFMOV - SQRT(CF90V**2 * CFOIJ* *2 / <CF90V**2 * CF0IJ**2))j 948 

MODULE: PIPEOl; 949 

NAME: MOV ; 980 

I/O LIST: UPSTREAM PROPERTY - VL21, 981 

DOWNSTREAM PROPERTY - MCHB; 982 

CMT: MAIN OXIDIZER VALVE; 983 

END MODULE 984 

EQUATION: WOINJ - WMOV * PRIMMOI s 988 

EQUATION: TTMOI - TTYL21 j 986 

* 987 

* PREP EXIT DENSITY 988 

* 989 

EQUATION: RHOPBPD - RHOPBSOj 960 

************************** 961 

* PREBURNER PUMP * 962 

************************** 963 

MODULE: PCMPOlj 964 

NAME: PRBPi . 968 

I/O LIST: INLET FLOW - PRBP, 966 

INLET PROPERTIES - VL21, 967 

EXIT PROPERTIES - PBPD. 968 

SHAFT - OH ; 969 

DESIGN VALUES: SND - 27240.8, 970 

TRQD - 3127.3, 971 

WD - 98.73, 972 

HDD - 74264.7, 973 

GEAR - 1.0; 974 

MAP: PMAP08; 978 

CMT: PREBURNER OXYGEN PUMP; 976 

END MODULE 977 

****************************************** 978 

* N ON-INERTIAL OXIDIZER LINE SIX * 979 

****************************************** 980 

MODULE: PIPEOl; 981 

NAME: 06 ; 982 

I/O LIST: INLET PROPERTIES - PBSO, 983 

EXIT PROPERTIES - VL20; 984 

DESIGN VALUES: CF - 0.813; 988 

CMT: PRBP RE CIRC. FLOW FROM VOLUME PBSO TO VOLUME 20; 986 

END MODULE 987 

***************************************** 988 

* FUEL PREBURNER OXIDIZER VALVE * 989 

***************************************** 990 
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* LUMP LINE. VALVE . AND INJECTOR RESISTANCES 991 

ELATION t RKLSFPOV - 0.62881; 992 

EQUATION: CFFPVX - 37.98 * AREAFPOV / SQRT( RKLSFPOV) ; 993 

EQUATION: CF011 - 17.186; 994 

EQUATION: CF11FP - SQRI( CFFPVX* *2 * CF01 1**2/ (CFFPVX* *2 ♦ CF011**2)>; 998 

* SCALE INJECTOR RESISTANCE WITH PRIMING FRACTION 996 

EQUATION: CFFPB - CFOFPB / PRIMFPB ; 997 

EQUATION: CFFPOV - SQRT(CF11FP**2 * CFFPB**2/(CF11FP**2 ♦ CFFPB**2>>; 998 

MODULE: PIPE01; 999 

NAME: FPOV; 1000 

I/O LIST: UPSTREAM PROPERTY - PBSO, 1001 

DOWNSTREAM PROPERTY - FPRBj 1002 

CMI: FUEL PREBURNER OXIDIZER VALVE; 1003 

END MODULE 1004 

EQUATION: WOFPB - WFPOV * PRIMFPB ; 1008 

EQUATION: TTFPBI - TTPBSO ; 1006 

********************************************* 1007 

* OXIDIZER PREBURNER OXIDIZER VALVE * 1008 

********************************************* 1009 

* LUMP LINE. VALVE, AND INJECTOR RESISTANCES 1010 

EQUATION: RKLSOPOV - 0.84173; 1011 

EQUATION: CFOPVX - 37.98 * AREAOPOV / SQRI< RKLSOPOV) ; 1012 

EQUATION: CF010 - 6.7420; 1013 

EQUATION: CF100P - SQRTt CFOPVX* *2 * CF010**2/(CF0PVX**2 ♦ CF01Q**2)); 1014 

* SCALE INJECTOR RESISTANCE WITH PRIMING FRACTION 1018 

EQUATION: CFOPB - CFOOPB / PRIMOPB ; 1016 

EQUATION: CFOPOV - SQRT(CF100P**2 * CFOPB**2/(CF100P**2 * CF0PB**2)); 1017 
MODULE: PIPE01; 1018 

NAME: OPOV; 1019 

I/O LIST: UPSTREAM PROPERTY - PBSO. 1020 

DOWNSTREAM PROPERTY - OPRB; 1021 

CMT: OXIDIZER PREBURNER OXIDIZER VALVE; 1022 

END MODULE 1023 

EQUATION: WOOPB - WOPOV * PRIMOPB ; 1024 

EQUATION: TTOPBI - TTPBSO ; 1028 

********************************************************************** 1026 

* OXYGEN SIDE DERIVATIVE MODULES * 1027 

********************************************************************** 1028 

* 1029 

************************************** 1030 

* INERTIAL OXIDIZER LINE ONE * 1031 

************************************** 1032 

MODULE: PIPE03; 1033 

NAME: 01 ; 1034 

I/O LIST: INLET PROPERTIES - OINK. 1038 
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EXIT PROPERTIES - VL17: 1036 

DESIGN VALUES: CF - 1049.2, 1037 

AREA - 11.00, 1038 

XLEN - 1128., 1030 

DLTZ - 1128.: 1040 

CMT: FLOW DERIVATIVE FROM OINK TO VOLUME 17; 1041 

END MODULE 1042 

************************************** 1043 

* INERTIAL OXIDIZER LINE IWO * 1044 

************************************** 104S 

MODULE: PIPEOO: 1048 

NAME: 02 ; • 1047 

I/O LIST: INLET PROPERTIES - VL17, 1048 

EXIT PROPERTIES - VL18; 1049 

DESIGN VALUES: CF - S274.6, 1080 

AREA - 11.00, 1081 

XLEN - 300.0: 1082 

CMT: FLOW DERIVATIVE FROM VOLUME 17 TO VOLUME 18; 1083 

END MODULE 1084 

**************************************** 1088 

* INERTIAL OXIDIZER LINE THREE * 1086 

**************************************** 1087 

MODULE: PIPEOO: 1088 

NAME: 03 : 1089 

I/O LIST: INLET PROPERTIES - VL19, 1060 

EXIT PROPERTIES - VL20; 1061 

DESIGN VALUES: CF - 1262.6, 1062 

AREA - 1.00, 1063 

XLEN - 8.330: 1064 

CMT: FLOW DERIVATIVE FROM VOLUME 19 TO VOLUME 20: 1068 

END MODULE 1066 

*************************************** 1067 

* INERTIAL OXIDIZER LINE FIVE * 1068 

*************************************** 1069 

MODULE: PIPEOO: 1070 

NAME: 08 ; 1071 

I/O LIST: INLET PROPERTIES - VL21, 1072 

EXIT PROPERTIES - VL22: 1073 

DESIGN VALUES: CF - 147.4, 1074 

AREA - 1.00, 1078 

XLEN - 66.67: 1076 

CMT: FLOW DERIVATIVE FROM VOLUME 21 TO VOLUME 22: 1077 

END MODULE 1078 

**************************** 1079 

* VOLUME SEVENTEEN * 1080 
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**************************** 1081 

MODULE: VOLMOO; 1082 

NAME: VL17; 1083 

I/O LIST: UPSTREAM PROPERTIES - OINK, 1084 

INLET FLOW - 01 , 1085 

EXIT FLOW - 02 , 1086 

DOWNSTREAM PROPERTIES - VL18, 1087 

qP0T - VL17 j 1088 

DESIGN VALUES: VOL - 12408., 1080 

QDOT - 0.0 t 1090 

CMTi DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 17; 1091 

END MODULE 1092 

*************************** 1093 

* VOLUME EIGHTEEN * 1094 

*************************** 1096 

MODULE: VOLMOO; 1096 

NAME: VL18; 1097 

I/O LIST: UPSTREAM PROPERTIES - VL17, 1098 

INLET FLOW - 02 , 1099 

EXIT FLOW - LPOP, 1100 

DOWNSTREAM PROPERTIES - LPOD, 1101 

QDOT - VL18; 1102 

DESIGN VALUES: VOL - 3300., 1103 

QDOT - 0.0 ; 1104 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 18; 1105 

END MODULE 1106 

*************************** 1107 

* VOLUME NINETEEN * 1108 

*************************** 1109 

MODULE: V0LM01; 1110 

NAME: VL19; 1111 

I/O LIST: UPSTREAM PROPERTIES - LPOD, LTOD, 1112 

INLET FLOW - LPOP, LPOT, 1113 

EXIT FLOW - 03 , 1114 

DOWNSTREAM PROPERTIES - VL20, 1115 

QDOT - VL19; 1116 

DESIGN VALUES: VOL -1771.0; 1117 

CUT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 19; 1118 

END MODULE 1119 

************************* 1120 

* VOLUME TWENTY * 1121 

************************* 1122 

EQUATION : EIPOGO - HTVL20 ; 1123 

MODULE: V0LM01; 1124 

NAME: VL20; 1125 
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I/O LIST: UPSTREAM PROPERTIES - VL19, PBSO. VL21. 1126 

INLET FLOW - 03 ,06 ,04 , 1127 

EXIT FLOW - BP OP, POGO, 1128 

DOWNSTREAM PROPERTIES - HPOD, POGO, 1120 

QDOT - VL20; 1130 

DESIGN VALUES: VOL -4936.0; 1131 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 20; 1132 

END MODULE 1133 

**************************** 1134 

* VOLUME IWENTYONE * 1138 

**************************** 1136 

MODULE: V0LM01; 1137 

NAME: VL21; 1138 

I/O LIST: UPSTREAM PROPERTIES - EPOD, 1139 

* INLET FLOW - EPOP, 1140 

EXIT FLOW - 04 . 08 . 07 . MOV , PRBP, 1141 

DOWNSTREAM PROPERTIES - VL20, VL22, POGO. VL21, PBPD, 1142 

()DOT - VL21; 1143 

DESIGN VALUES: VOL -1260.0; 1144 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 21; 1148 

END MODULE 1146 

**************************** 1147 

* VOLUME IWENTYIWO * 1148 

**************************** 1149 

MODULE: VOLMDO; 1180 

NAME: VL22; 1181 

I/O LIST: UPSTREAM PROPERTIES - VL21, 1182 

INLET FLOW - 08 , 1183 

EXIT FLOW - LPOT, 1184 

DOWNSTREAM PROPERTIES - LTOD, 1188 

QDOT - VL22; 1186 

DESIGN VALUES: VOL - 998.0, 1187 

QDOT - 0.0 ; 1158 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME 22; 1189 

END MODULE 1160 

********************************************** 1161 

* PREBURNER OXIDIZER SPLITTER VOLUME * 1162 

********************************************** 1163 

MODULE: V0LM01; 1164 

NAME: PBSO; 1168 

I/O LIST: UPSTREAM PROPERTIES - PBPD, 1166 

INLET FLOW - PRBP, 1167 

EXIT FLOW - 06 , OPOV, FPOV, 1168 

DOWNSTREAM PROPERTIES - VL20, PBSO. PBSO. 1169 

{JDOT - PBSO; 1170 
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DESIGN VALUES: VX -347.0; 1171 

CMT: DENSITY AND INTERNAL ENERGY DERIVATIVES FOR VOLUME PESO; 1172 

END MODULE H 73 

*************************** 1174 

* POCO SUPRESSION * 1178 

*************************** 1176 

MODULE: POCOOQ; U77 

NAME t POCO i 1178 

I/O LIST: INLET HELIUM FLOW - HE4 . 1179 

LIQUID OXIDIZER PROP - VL20, 1180 

GASEOUS OXIDIZER FLOW - 07 , 1181 

EXIT FLOW - 08 . 1182 

EXIT PROP - VL17{ 1183 

DESIGN VALUES: VX - 2000. | 1184 

Off: POCO SUPPRESSOR: H8« 

END MODULE 1186 

*********************************************************************** 1187 

* HOT GAS SIDE NON-DERIVATIVH MODULES * 1188 

*********************************************************************** 1189 

* 1190 

***************************************** 1191 

* NON -INERTIAL HOT GAS LINE TWO * 1192 

***************************************** 1193 

MODULE: PIPEOB; HM 

NAME: HC2 } H9« 

I/O LIST: INLET PROPERTIES - OPRB, 1196 

EXIT PROPERTIES - OTBP; 1197 

DESIGN VALUES: AREA - 12.34, 1198 

PELS - 1.000; 1198 

CMT: CALCULATES OPRB PRESSURE FROM OTBP PRESSURE AND FLOW; 1200 

END MODULE l 201 

****************************************** 1202 

* NON-INERTIAL HOT CAS LINE FOUR * 1203 

****************************************** 1204 

MODULE: PIPE02: 1208 

NAME: HG4 ; l 206 

I/O LIST: INLET PROPERTIES - OTBP, 1207 

EXIT PROPERTIES - OSF ; 1208 

DESIGN VALUES: AREA - 0.0667, 1209 

RKLS - 1.0; 1210 

CMT: CALCULATES FLOW FROM VOLUME OTBP TO VOLUME OSF; 1211 

END MODULE 1212 

************************************************************* 1213 

* HIGH PRESSURE OXIDIZER TURBINE DISCHARGE PRESSURE * 1214 

************************************************************* 1216 
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EQUATION > ZHPOT - ZOIBP; 1216 

****************************************** 1217 

* HIGH PRESSURE OXIDIZER TURBINE * 1218 

****************************************** 1210 

MODULBi TURBO 1; 1220 

NAMEt HPOT: 1221 

I/O LISTt INLET PROPERTIES - OTBP, 1222 

BXIT PROPERTIES - HTOD, 1223 

SHAFT WORE - OH ; 1224 

DESIGN TALUESs ETAD - 0.7S2, 1228 

PRD - 1.81, 1226 

PSID - 1.00, 1227 

SND - 27240.8, 1228 

AREA - 1.0 1229 

DIAM - 10.19 , 1230 

GEAR - 1.0 i 1231 

MAPs TBMP08; 1232 

CM! i HIGH PRESSURE OXYGEN TURBINE i 1233 

END MODULE 1234 

******************************************************* 1238 

* N ON-INERTIAL HOT GAS OVERBOARD LEAKAGE LINE * 1236 

******************************************************* 1237 

MODULES PIPE02: 1238 

NAMES OLE ; 1230 

I/O. LISTs INLET PROPERTIES - OSF . 1240 

EXIT PROPERTIES - AMB j 1241 

DESIGN VALUESs AREA - 0.0, 1242 

RKLS - 1.0; 1243 

CMTt CALCULATES LEAKAGE FLOW FROM VOLUME OSF; 1244 

END MODULE 1248 

***************************************** 1246 

* NON- INERTIAL HOT CAS LINE SIX * 1247 

***************************************** 1248 

MODULES PIPE08; 1240 

NAMES HG6 ; 1280 

I/O LIST: INLET PROPERTIES - OSF , 1281 

EXIT PROPERTIES - MFI ; 1282 

DESIGN VALUES: AREA - 11.276, 1283 

RKLS - 1.000; 1284 

CMT: CALCULATES OSF PRESSURE FROM MFI PRESSURE AND FLOW; 1288 

END MODULE 1286 

***************************************** 1287 

* NON- INERTIAL ROT CAS LINE ONE * 1288 

***************************************** 1289 

MODULE: PTPE08; 1260 
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NAME: HC1 : 1261 

I/O LIST: INLET PROPERTIES - FPRB, 1262 

EXIT PROPERTIES - FTBP; 1263 

DESIGN VALUES: AREA - 31.86, 1264 

RKLS - l.OOO; 1265 

CMT: CALCULATES FPRB PRESSURE FROM FTBP PRESSURE AND FLOW; 1266 

END MODULE 1267 

******************************************* 1268 

* N ON-INERTIAL HOT GAS LINE THREE * 1269 

******************************************* 1270 

MODULE: PIPE02; 1271 

NAME: HC3 ; 1272 

I/O LIST: INLET PROPERTIES - FTBP, 1273 

EXIT PROPERTIES - FSF ; 1274 

DESIGN VALUES: AREA - 0.0, 1275 

RKLS - 1.0; 1276 

CMT: CALCULATES FLOW FROM VOLUME FTBP TO VOLUME FSF; 1277 

END MODULE 1278 

********************************************************* 1279 

* HIGH PRESSURE FUEL TURBINE DISCHARGE PRESSURE * 1280 

********************************************************* 1281 

EQUATION: ZHPFT - ZFTBP; 1282 

************************************** 1283 

* HIGH PRESSURE FUEL IURBINB • 1284 

************************************** 1285 

MODULE: TURBO 1; 1286 

NAME: HPFT; 1287 

I/O LIST: INLET PROPERTIES - FTBP, 1288 

EXIT PROPERTIES - HTFD, 1289 

SHAFT WORK - FH ; 1290 

DESIGN VALUES: ETAD - 0.804, 1291 

PRD - 1.45, 1292 

PSID - 1.00, 1293 

SND - 34189.8, 1294 

AREA - 1.0 , 1295 

DIAM - 10. 19 , 1296 

CEAR - 1.0 ; 1297 

MAP: IBMP03; 1298 

CMT: HIGH PRESSURE FUEL TURBINE; 1299 

END MODULE 1300 

****************************************** 1301 

* NON- INERTIAL HOT CAS LINE FIVE * 1302 

****************************************** 1303 

MODULE: PIPH08; 1304 

NAME: HC5 ; 1305 
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I/O LIST: INLET PROPERTIES - FSF . 1306 

EXIT PROPERTIES - MFI ; 1307 

DESIGN VALUES: AREA - 21.36, 1308 

RELS - 1.000} 1300 

CMTs CALCULATES FSF PRESSURE FROM MFI PRESSURE AND FLOW; 1310 

END MODULE 1311 

*************************************************** 1312 

* NON-INERTIAL HOT CAS FUEL INJECTOR LINE * 1313 

*************************************************** 1314 

MODULE: PIPE08; 1316 

NAME: FINJ; 1316 

I/O LIST: INLET PROPERTIES - MFI , . 131" 

EXIT PROPERTIES - MCHB; 1318 

DESIGN VALUES: AREA - 27.90, 1319 

RELS - 1.000; * 1320 

CMT: CALCULATES MFI PRESSURE FROM MCHB PRESSURE AND FLOW; 1321 

END MODULE 1322 

****************************** 1323 

* NOZZLE PERFORMANCE * 1324 

****************************** 1326 

MODULE: NOZLOO; 1326 

NAME: NQZL; 1327 

I/O LIST: INLET PROPERTIES - MCHB, 1328 

EXIT PROPERTIES - AMB ; 1329 

DESIGN VALUES: AREA - 82.06, 1330 

AR - 77.5, 1331 

CS - 0.98; 1332 

CMT: NOZZLE PERFORMANCE CALCUALTION; 1333 

END MODULE 1334 

********************************************* 1336 

* NOZZLE HOT CAS SIDE HEAT TRANSFER * 1336 

********************************************* 1337 

EQUATION : AREANZLC - AREANOZL ; 1338 

MODULE: QN0Z01; 1339 

NAME: QNOZ; 1340 

I/O LIST: HOT GAS PROPERTIES - MCHB, 1341 

QDOT - N0Z1, N0Z2, 1342 

IM - Mill, MTL2, 1343 

NOZZLE AREA - NZLG; 1344 

DESIGN VALUES: RCRV - 6.1; 1346 

CMT: NOZZLE HEAT TRANSFER RATES; 1346 

END MODULE 1347 

********************************************** 1348 

* CHAMBER HOT GAS SIDE HEAT TRANSFER * 1349 

********************************************** 1380 
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MODULE : QCHMOl; 1381 

NAME: QCHM; 1382 

I/O LIST: HOT GAS PROPERTIES - MCHB, 1383 

qDOT - CHMlt CHM2, 1384 

IM - MTL3, MTL4, 1388 

NOZZLE AREA - NZLG; 1386 

DESIGN VALUES: RCRV - 6.1{ 1387 

CMT: CHAMBER HEAT TRANSFER RATES ; 1388 

END MODULE 1360 

*********************************************************************** 1360 

* HOT CAS SIDE DERIVATIVE MODULES * 1361 

*********************************************************************** 1362 

* 1363 

****************************** 1364 

* OXIDIZER PREBURNER * 1368 

****************************** 1366 

MODULE: PBRNOlj 1367 

NAME: OPRBj 1368 

I/O LIST: FUEL FLOW - FOPB, 1369 

FUEL PROPERTIES - PBSF. 1370 

OXIDIZER FLOW - OOPB, 1371 

OXIDIZER PROPERTIES - OPBI, 1372 

HELIUM FLOW - HE2 . 1373 

HELIUM PROPERTIES - HEIR, 1374 

EXIT FLOW - HC2 ; 1378 

DESIGN VALUES: VOL - 347. , 1376 

OFBL - 0.08 , 1377 

OFLT - 0.4 , 1378 

ILIT - 1 ; 1379 

CMT: PRESS., TEMP. , OXID. FRAC., AND BE. FRAC. DERIV. FOR OPRB: 1380 

END MODULE 1381 

******************************************* 1382 

* OXIDIZER TURBINE BT-PASS VOLUME * 1383 

******************************************* 1384 

MODULE: V0LM02: 1388 

NAME: OTBPj 1386 

I/O LIST: UPSTREAM PROPERTIES - OPRB, 1387 

INLET FLOW - HG2 . 1388 

EXIT FLOW - HPOT, HG4 . 1389 

DOWNSTREAM PROPERTIES - OSF , OSF , 1390 

JJDOT - OTBPj 1391 

DESIGN VALUES: VOL - 600.0; 1392 

CMT: PRESS., TEMP., OXID. FRAC., AND HE. FRAC. DERIV. FOR OTBP; 1393 

END MODULE 1394 

**************************************************** 1398 
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* OXIDIZER PRHBORNER SECONDARY FLOW VOLUME * 1396 

**************************************************** 1397 

MODULE: V0LM02| 1398 

NAME: OSF '; 1338 

I/O LIST: UPSIREAM PROPERTIES - HTOD, OTBP, PBSF, 1400 

INLET FLOW - SPOT, HG4 , OTC , 1401 

EXIT FLOW - HG0 , OLK , 1402 

DOWNSTREAM PROPERTIES - MFI , AMB . 1403 

qDOT - OSF ; 1404 

DESIGN VALOESi VOL - 800.0: 1408 

CMTt PRESS., IEMP., OXID. FRAC., AND BE. FRAC. DERIV. FOR OSF; 1408 

END MODULE 1407 

************************** 1408 

* FUEL PREBURNER * 1409 

************************** • 1410 

MODULE! PBRNOls 1411 

NAME! FPRB; 1412 

I/O LIST: FUEL FLOW - FFPB, 1413 

FUEL PROPERTIES - PBSF, 1414 

OXIDIZER FLOW - OFPB, 1415 

OXIDIZER PROPERTIES - FPBI, 1416 

HELIUM FLOW - HE1 , 1417 

HELIUM PROPERTIES - HEIR, 1418 

EXIT FLOW - HG1 i 1419 

DESIGN VALUES! VOL - 347.0, 1420 

OFBL - 0.08 , 1421 

OFLT - 0.4 , 1422 

ILIT - 1 : 1423 

CMTt PRESS., TEMP., OXID. FRAC., AND HE. FRAC. DERIV. FOR FPRB: 1424 

END MODULE 1425 

*************************************** 1426 

* FUEL TURBINE BY-PASS VOLUME * 1427 

*************************************** 1428 

MODULH: V0LMD2; 1429 

NAME t FTBP i 143Q 

I/O LIST! UPSTREAM PROPERTIES - FPRB, 1431 

INLET FLOW - HG1 , 1432 

EXIT FLOW - HP FT, HG3 , 1433 

DOWNSTREAM PROPERTIES - FSF , FSF , 1434 

QDOT - FTBP* 1435 

DESIGN VALUES: VOL - 500.0; 1436 

CMT: PRESS., TEMP., OXID. FRAC., AND HE. FRAC. DERIV. FOR FTBP; 1437 

END MODULE 1438 

************************************************ 1439 

* FUEL PREBURNER SECONDARY FLOW VOLUME * 1440 
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************************************************ 1441 

EQUATION t 0FRVL3 - 0.0 ; 1442 

EQUATION t HFRVL3 - 0.0 ; 1443 

MODULE i V0LM02; 1444 

NAME i FSF ; 1445 

I/O LISTi UPSTREAM PROPERTIES - HTFD, FTBP, VL3 , 1446 

INLET FLOW - HP FT, HG3 , FTC , 1447 

EXIT FLOW - ECS . 1448 

DOWNSTREAM PROPERTIES - MFI . 1449 

QDOT - FSF ; 1450 

DESIGN VALUES: VOL - 500.0; 1451 

CMT: PRESS.. TEMP., OXID. FRAC. , AND HE. FRAC. DERIV. FOR FSF; 1462 

END MODULE 1463 

************************************* 1454 

* MAIN FUEL INJECTOR VOLUME * 1455 

************************************* i486 

MODULE: V0LM02; 1457 

NAME: MFI | 1458 

I/O LIST: UPSTREAM PROPERTIES - FSF , OSF , VL16, 1459 

INLET FLOW - EGA . HC6 , FSLV. 1460 

EXIT FLOW - FINJ, 1461 

DOWNSTREAM PROPERTIES - MCHB. 1462 

QDOT - FMCO, OMCOj 1463 

DESIGN VALUES: VOL -4210.0; 1464 

CMT: PRESS., TEMP., OXID. FRAC., AND HE. FRAC. DERIV. FOR MFI; 1465 

END MODULE 1466 

*********************************** 1467 

* MAIN CHAMBER COMBUSTION * 1468 

*********************************** 1469 

MODULE: MCBB01; 1470 

NAME: MCHB; 1471 

I/O LIST: FUEL FLOW - F15 , 1472 

FUEL PROPERTIES - VL16, 1473 

* FUEL ICNITER FLOW - FIG , 1474 

OXIDIZER FLOW - OINJ, 1475 

OXIDIZER PROPERTIES - MDI , 1476 

OXIDIZER ICNITER FLOW- OIG , 1477 

MFI INLET FLOW - FINJ. 1478 

MFI INLET PROPERTIES - MFI , 1479 

HELIUM FLOW - HE3 , 1480 

HELIUM PROPERTIES - HETK, 1481 

EXIT FLOW - NZLG, 1482 

QDOT - CHM1, COM2, 1483 

ILITPB1 - FPRB, 1484 

ILITPB2 - OPRB; 1485 
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TTBE I/O CONTROL 


PACE 34 


CONFIG. INPUT 


DESIGN VALUES i VOL -2228.0, 1486 

OFBL - 0.08 , 1487 

OFLT - 0.4 , 1488 

ILIT - 1 t 1489 

CMTt PRESS., IEMP. , OXID. FRAC., AND BE. FRAC. DERIV. FOR MCHBj 1490 
END MODULE 1491 

1492 

* NOZZLE NODE ONE METAL TEMPERATURE DYNAMICS * 1493 

****************************************************** 1494 

MODULE t MHTLOO; 1495 

NAME! MTLli 1498 

I/O LIST: QDOT - GROT, N0Z1; 1497 

DESIGN VALUES: IMTL - 1 , 1498 

RM - 18.2; 1499 

CMT: METAL TEMPERATURE DERIVATIVE FOR NOZl/MTLlj 1800 

END MODULE 1601 

****************************************************** 1602 

* NOZZLE NODE TWO METAL TEMPERATURE DYNAMICS * 1603 

****************************************************** 1604 

MODULE: METLOO; 1606 

NAME: MTL2; 1606 

I/O LIST: tJDOT - 7H0T, N0Z2; 1607 

DESIGN VALUES: IMTL - 1 . 1508 

RM - 48.8; 1609 

CMT: METAL TEMPERATURE DERIVATIVE FOR N0Z2/MTL2; 1610 

END MODULE 1611 

******************************************************* 1612 

* CHAMBER NODE ONE METAL TEMPERATURE DYNAMICS * 1613 

******************************************************* 1814 

MODULE: METLOO; 1618 

NAME: MTL3; 1616 

I/O LIST: QDOT - 11HT, CHM1; 1517 

DESIGN VALUES: IMTL - 1 , 1618 

RM - 31.0 ; 1619 

CMT: METAL TEMPERATURE DERIVATIVE FOR CHM1/MTL3; 1620 

END MODULE 1621 

******************************************************* 1822 

* CHAMBER NODE TWO METAL TEMPERATURE DYNAMICS * 1623 

******************************************************* 1624 

MODULE: METLOO; 1628 

NAME: MIL4; 1626 

I/O LIST: qDOT - 12BT, CHM2; 1627 

DESIGN VALUES: IMIL - 1 , 1628 

RM - 31.0 ; 1629 

CMT: METAL TEMPERATURE DERIVATIVE FOR CHM2/MTL4; 1630 
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TIBS W/0 CONTROL 


PACE 38 


CONFIG. INPOT 


END MODULE 1831 

****************************************************** 1832 

* LOW PRESSURE FUEL TURBOPUMP ROTOR DYNAMICS * 1833 

****************************************************** 1834 

MODULE: ROIROO: 1838 

NAME: FL j 1836 

I/O LIST: TORQUE - LPFT, LPFP; 1837 

DESICN VALUES: PMOM - 1.2863 t 1838 

CMT: ROTOR SPEED DERIVATIVE FOR LPFT/LPFP; 1839 

END MODULE 1840 

******************************************************* 1841 

* HIGH PRESSURE FUEL TURBOPUMP ROTOR DYNAMICS * 1842 

******************************************************* 1843 

MODULE: ROIROO; 1844 

NAME: FH » 1848 

I/O LIST: TORQUE - HP FT, HPFP; 1846 

DESIGN VALUES : PMOM - 2.8608 ; 1847 

CMT: ROTOR SPEED DERIVATIVE FOR HP FT/HP FP ; 1848 

END MODULE 1849 

********************************************************** 1860 

* LOW PRESSURE OXIDIZER TURBOPUMP ROTOR DYNAMICS * 1851 

********************************************************** 1682 

MODULE: ROIROO; 1563 

NAME: OL ; 1884 

I/O LIST: TORQUE - LPOT, LPOP; 1868 

DESIGN VALUES : PMOM - 2.3900 ; 1666 

CMT: ROTOR SPEED DERIVATIVE FOR LPOT/LPOP; 1867 

END MODULE 1868 

*********************************************************** 1869 

* high PRESSURE OXIDIZER TURBOPUMP ROTOR DYNAMICS * 1560 

*********************************************************** 1661 

MODULE: ROIROO; 1862 

NAME: OB ; 1863 

I/O LIST: TORQUE - BPOT, BPOP, PRBP; 1864 

DESICN VALUES : PMOM - 1.4496 ; 1868 

CMT: ROTOR SPEED DERIVATIVE FOR HPOT/HPOP; 1866 

END MODULE 1867 

* 1668 

END SYSTEM INSIDE 1869 

* 1870 

— * 1871 

* * 1872 

* BELOW THE ITERATION LOOP * 1873 

« * 1874 

* 1678 
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UBS I/O CONTROL 


PACE 36 


CONFIG. INPUT 


* 1876 

DEFINE SYSTEM BELOV 1877 

* 1878 

*************************************************************** 1879 

* EULER INTEGRATE THE DENSITY OF THE INJECTOR VOLUMES * 1880 

*************************************************************** 1881 

EQUATION t VOLMOI -1188.9 ; 1882 

EQUATION t VOLOPBI - 41. B» 1883 

EQUATION i VOLFP.BI - 81.8; 1884 

EQUATION * RHOMOI - RHOMOI ♦ DT * (1M0V - VOINJ) / VOLMOI ; 1888 

EQUATION * RHOOPBI - RHOOPBI * DT * (WOPOV - WOOPB) / VOLOPBI; 1886 

EQUATION i RHOFPBI - RHOFPBI * DT * (VFPOV - VDFPB) / VOLFPBI; 1S87 

****************************************************** 1888 

* LOW PRESSURE FUEL TURBOPUMP ROTOR BREAK AWAY - * 1889 

****************************************************** 1890 

MODULE; R0IRO1; 1891 

NAME; FLBR; 1892 

I/O LIST: TORQUE - LPFT, LPFP, 1893 

SHAFT - FL ; 1894 

DESIGN VALUES: BTQ - 40.0; 1896 

CMT: ROTOR BREAK AWAY FOR LPFT/LPFP; 1896 

END MODULE 1897 

******************************************************* 1898 

* HIGH PRESSURE FUEL TURBOPUMP ROTOR BREAK AWAY - * 1899 

******************************************************* 1600 

MODULE: ROTROl; 1601 

NAME: FHBR; 1602 

I/O LIST: TORQUE - HPPT, BPFP, 1603 

SHAFT - FH i 1604 

DESIGN VALUES : BTQ - 160.0; 1608 

CMT: ROTOR BREAK AWAY FOR HPFT/HPFP; 1606 

END MODULE 1607 

********************************************************** 1608 

* LOW PRESSURE OXIDIZER TURBOPUMP ROTOR BREAK AWAY - * 1609 

********************************************************** 1610 

MODULE: ROTROl; 1611 

NAME: OLBR; 1612 

I/O LIST; TORQUE - LPOT, LPOP, 1613 

SHAFT - X ; 1614 

DESIGN VALUES: BTQ - 80.0; 1616 

CMT: ROTOR BREAK AWAY FOR LPOT/LPOP; 1616 

END MODULE 1617 

*********************************************************** 1618 

* high PRESSURE OXIDIZER TURBOPUMP ROTOR BREAK AWAY - * 1619 

*********************************************************** 1620 
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TTBH W/0 CONTROL 


PAGE 37 


CONFIG. INPUT 


MODULE: ROTROl; 1621 

NAME: OHBR; 1622 

I/O LIST: TORQUE - HPOT, HP OP, 1623 

SHAFT - OH ; 1624 

DESIGN VALUES : BTQ - 90.0; 1628 

CMT: ROTOR BREAK AWAY FOR HPOT/HPOP; 1626 

END MODULE 1627 

* 1628 

END SYSTEM BELOV 1629 
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